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FOREWORD 


Foreword 


Ralph  R.  Goodman 
Director 

SAC L ANT  ASW  Research  Centre 
La  Spezia,  Italy 


Geologists  and  geophysicists  have  used  both  acoustic  and  seismic  methods  to  study 
the  seabed  by  considering  the  propagation  of  signals  arising  from  both  natural 
seismic  events  and  man-made  impulsive  sources.  Although  significant  advances 
have  been  made  in  instrumentation,  such  as  long  towed  geophysical  arrays,  air 
guns  and  ocean  bottom  seismometers,  the  picture  of  the  seafloor  is  s till  far  from 
complete. 

Underwater  acoustics  concerns  itself  today  with  the  phenomena  of  propagation 
and  noise  at  frequencies  and  ranges  that  require  an  understanding  of  acoustic 
interaction  at  both  of  its  boundaries,  the  sea  surface  and  seafloor,  over  depths 
ranging  from  tens  to  thousands  of  meters.  Much  of  the  earlier  higher  frequency 
(>  1  kHz)  work  included  the  characterization  of  the  seafloor  in  regimes  of 
reflection  coefficients  which  were  empirically  derived  from  surveys.  The  results 
of  these  studies  met  with  only  limited  success,  confined  as  they  were  to  those 
areas  where  survey  data  existed  and  lacking  a  physical  understanding  of  the 
processes  of  reflection  and  scattering.  The  development  in  recent  years  of 
systematic  methods  of  measuring  the  physical  and  acoustical  characteristics  of  the 
seabed  have  led  to  a  much  better  understanding  of  propagation.  However  these 
techniques  are  exceedingly  slow  and  expensive  and  do  not  fulfill  the  needs  of 
present  day  acoustics.  The  relatively  recent  recognition  of  the  importance  of  the 
role  that  shear  and  evanescent  waves  in  the  seabed  have  on  waterborne  sound 
has  been  a  major  factor  in  the  recognition  of  the  need  for  acousticians  and 
seismologists  to  exchange  and  discuss  their  common  goals  and  developments. 

The  realization  of  the  benefits  that  could  be  achieved  from  such  an  exchange  led 
our  Centre  to  sponsor  the  Symposium  on  Ocean  Seismo-Acoustics,  with  partial 
support  from  the  Office  of  Naval  Research,  London.  It  was  held  at  La  Spezia, 
Italy,  June  10-14,  1985  and  its  purpose  was  to  enhance  cross-disciplinary  dialogue 
and  generate  ideas  for  new  research  directions. 

The  collection  of  papers  originating  from  the  Symposium  have  been  published  as  a 
volume*  of  the  NATO  Conference  Series,  with  the  papers  organized  into  eiph* 
topics: 

-  Propagation:  Theoretical  Developments 

-  Propagation:  Modelling  and  Experiments 

-  Propagation:  Fluctuations 

-  Boundary  Scattering 

-  Sediments:  Seismo-Acoustic  Waves  and  Geoacoustic  Characteristics 

-  Seismo-Acoustic  Waves:  Propagation,  Attenuation,  and  Inverse 

Techniques 

-  Seismo-Acoustic  Noise 

-  Technology  and  New  Approaches  in  Seismo-Acoustic  Measurements 

This  present  document  is  a  compilation  for  normal  NATO  purposes  of  the  papers 
presented  at  the  Symposium  by  Centre  staff. 


*  NATO  Conference  Series.  IV,  Marine  sciences;  v.16:  "OCEAN  SEISMO-ACOUSTICS: 
Low-frequency  underwater  acoustics" ,  edited  by  Tuncay  Akal  and  Jonathan  M.  Berkscr., 
(including  bibliographies  and  indexes).  Published  by  Plenum  Press,  New  York. 
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AKAL  et  al:  Ambient  seafloor  noise 


TIME  DEPENDENCE  OF  INFRASONIC  AMBIENT  SEAFLOOR  NOISE 
ON  A  CONTINENTAL  SHELF 

T.  Akal,  A.  Barbagelata,  G.  Guidl,  and  M.  Snoek 

SACLANT  ASW  Research  Centre 
Viale  San  Bartolomeo  400 
1-19026  La  Spezla,  Italy 


ABSTRACT 

SACLANTCEN  is  conducting  a  long  term  programme  of  ambient  noise 
measurements  that  cover  the  frequency  band  from  1  to  100  Hz.  A  sea-bed 
sensor  consisting  of  three  seismometers  and  a  hydrophone  has  been  deployed 
close  to  the  edge  of  the  continental  break  In  water  of  100  m  depth.  The 
data  are  being  systematically  collected  and  telemetered  via  radio  link  to 
SACLANTCEN.  The  initial  data  (105  days)  have  been  used  to  Investigate  the 
noise  characteristics  and  its  statistics. 


INTRODUCTION 

Ambient  noise  in  the  ocean  has  many  sources.  As  a  result,  long  term 
measurements  are  required  to  identify  these  sources  and  define  the  noise 
characteristics  of  an  ocean  area.  Studies  of  ambient  noise  have  shown 
that  seismic  acoustic  energy  in  certain  parts  of  the  noise  spectrum  can  be 
related  to  various  known  sources  [1  to  13).  Figure  1  identifies  several 
well  known  sources  that  transfer  energy  basically  from  biological, 
meteorological,  geological  or  cultural  origin  into  the  underwater  seismo- 
acoustic  noise  field. 

SACLANTCEN  has  begun  a  programme  to  measure  the  time  dependence  of 
ambient  noise  along  the  Ligurian  shelf  edge  (Fig.  2a).  The  frequency  band 
of  interest  is  1  to  100  Hz.  Initially  both  short  (33  hr)  and  long 
(105  days)  term  measurements  were  made. 

Two  major  factors  affect  ambient  noise  characteristics,  the  meteorolo- 
gic  conditions  of  the  area  and  the  distribution  of  the  noise  sources 
around  the  measurement  site.  In  support • of  the  interpretations  of  noise 
data,  environmental  and  shipping  data  that  could  affect  noise  charac¬ 
teristics  have  been  Incorporated  Into  the  study. 
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INSTRUMENTATION 

The  instrumentation  system  used  for  these  measurements  consisted  of 
the  following  components  (Fig.  2b): 

-  the  ocean  bottcm  seismometer 

-  the  radio  transceiver  (mounted  on  a  surface  buoy) 
the  wave-rider  buoy  (attached  to  the  surface  buoy) 
the  radio  repeater  (mounted  on  a  mountain  top) 

the  receiver  an..  control  unit 

-  the  Interface  to  the  digital  computer 
the  computer  system. 

Signals  from  the  sensors  were  immediately  converted  into  digital  form. 
They  suffered  no  significant  degradation  when  transmitted  through  the 
various  components  >.o  SACLANTCEN. 

Ocean  Bottom  Seismometer  (OBS) 

The  OBS  consisted  of  geoacoustical  and  auxiliary  sensors,  preampli¬ 
fiers  and  A/D  electronics  that  are  contained  in  a  disc-shaped  housing. 
Power  is  provided  by  a  battery  installed  in  the  surface  buoy  and 
transmitted  via  a  coatial  cable. 

The  geoacoustic  sensors  comprised  three  active  1  Hz  geophones 
(Teledyne  S-500)  mounted  orthogonally  and  an  omnidirectional  hydrophone. 
These  sensors  can  be  operated  at  any  angle  and  need  no  locking  for 
transportation  and  handling.  They  have  a  velocity  sensitivity  of 
450  V-sec/m  and  a  -5  dB  passband  of  1  to  100  Hz.  In  addition  to  these 
sensors,  the  system  includes  auxiliary  sensors:  two  tilt  meters  to  moni¬ 
tor  the  inclination,  a  compass  to  determine  orientation,  and  a  thermistor 
to  measure  temperature. 

The  four  geaacoustic  sensors  are  sampled  at  600  Hz  per  channel.  The 
data  is  converted  into  a  floating  point  format  and  combined  with  low  speed 
data  from  positional  sensors.  No  gain  adjustment  is  necessary  during  the 
experiment  because  the  instantaneous  floating  point  format  allows  the 
achievement  of  a  very  high  dynamic  range  (120  dB).  Calibration  can  be 
carried  out  by  injecting  stabilized  AC  current  into  the  calibration  coil 
of  the  geophones.  The  data  is  transmitted  in  biphase  code  to  the  surface 
buoy  via  a  coaxial  cable  at  38400  bits/sec. 


The  Radio  Buoy 

The  radio  buoy,  mounted  or.  a  4  m  diameter  circular  buoy  moored  in 
100  m  depth,  about  8  km  from  the  coast,  contained  recharge-ble  batteries 
(i.e.  power  for  32  hours  of  continuous  data  acquisition),  a  FM  modulated 
VHP  radio  transmitter  (220  MKz/4  W)  and  a  standby  radio  receiver  for  power 
ON/OFF  commands . 


The  Repeater 


The  unattended  radio  repeater  was  installed  on  510  m  high 
Mt.  Castellano,  to  overcome  topographic  Interference  between  the  OBS  radio 
buoy  and  SACLANTCEN.  The  repeater  received  and  retransmitted  OBS  radio 
buoy  signals  ar  different  frequencies.  A  dedicated  telephone  line  between 
the  repeater  and  SACLANTCEN  permitted  the  relay  of  orders  to  the  OBS  and 
the  monitoring  of  the  performance  of  the  OBS  instruments. 


1 


SACLANTCEN  CP-37 


1-3 


AKAL  et  al:  Ambient  seafloor  noise 


The  Recclvor  and  ConLrol  Unit 


This  untt,  Installed  at  SACLANTCKN,  allowed  the  performance  of  the 
following  functions :  reception,  demodulation  and  forwarding  of  the  com¬ 
puter  Interface  signals,  monitoring  of  the  performance  of  the  geoacouatlc 
and  environmental  sensors  mounted  on  the  0B8 ,  through  dlgltal-to-analogue 
conversion  and  display  of  the  transmitted  signals,  monitoring  and  manual 
control  of  the  performance  of  the  OBS  and  repeater  station. 


The  Computer  Interfac. c 

This  unit,  Installed  In  the  computer  area,  received  the  baseband 
signal  from  the  receiver  unit.  Bits  and  words  were  synchronized  and 
transferred  to  the  digital  computer  In  a  parallel  format.  Commands  were 
received  from  the  computer  which  switched  the  full  chain  on  and  off 
according  to  a  pre-set  time  schedule. 


The  Computer  System 

Acquisition  and  preprocessing  were  performed  on  a  Hewlett  Packard 
21  MX  computer  with  a  disc  unit  and  other  peripherals.  Data  were  imme¬ 
diately  stored  on  discs  for  display,  plot  and  preprocessing  options  and 
for  transfer  to  digital  tape. 


MEASUREMENT  AND  ANALYSIS  TECHNIQUES 

This  system  was  used  to  measure  ambient  sea  noise  for  105  days 
(December  1984  to  March  1985).  Automatic  acquisitions  of  ambient  noise  of 
5  minutes  duration  were  made  dally;  wave-rider  data  and  meteorological 
data  were  obtained  over  longer  periods.  Figure  3,  a  sample  record, 
illustrates  the  variation  of  the  ambient  noise  levels  at  different  sensors 
for  different  meteorological  conditions,  together  with  the  wave-rider 
recordings.  Spectral  estimates  of  these  signals  were  made  by  applying 
standard  techniques  [14]  to  25-sec  time-windows  over  a  1  to  100  Hz  fre¬ 
quency  band,  and  averaging  the  power  spectral  density  levels  over  a  1  Hz 
frequency  bandwidth.  Figure  4  shows  the  power  spectral  density  of  the 
sample  signals  of  Fig.  3. 


Noise  Fluctuations  in  33-hour  Period 


Measurements  of  5,  10,  30,  60  minutes  and  24  hours  duration  were  con¬ 
ducted  to  study  both  the  short  and  long-time  variations  of  ambient  noise. 
The  power  spectral  density  contours  in  a  frequsney/rangs  plans  provide  the 
eimplaet  display  of  the  behaviour  of  the  ambient  noise  field  fluctustions, 
which  are  a  function  of  both  frequency  and  time. 

Figure  5  shows  an  example  of  the  changes  of  the  ambient  noise  spectrum 
levels  ovsr  a  33-hour  period  at  different  sensors.  On  21  March  85,  the 
sea  surface  wave  field  at  the  measurement  site  wee  flat,  under  e  8WW  wind 
of  3  n. rai/hr.  During  the  day  the  wind  increased  to  15  n. ml/hr,  causing 
waves  as  high  as  80  cm.  Twenty-seven  hours  after  the  beginning  of  the 
data  acquisition,  wind  direction  changed  to  North  at  2  n. mi/hr.  The 
characteristics  of  the  ambient  noise  spectrum  levels  evidsnt  in  Fig.  5 
provide  a  natural  division  of  the  noise  field  into  two  distinct  regions: 
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Fig.  3  Staple  recordings:  signals  at  various  OBS  sensors,  nave-rider 
and  nctcoro logic  information  under  different  meteorological 
condi c ions . 
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In  the  low  frequency  (1  to  30  Hz)  region,  horizontal  sensors  generally 
have  a  higher  energy  content  than  the  vertical  sensor  and  the  omni¬ 
directional  hydrophone  has  a  much  steeper  fall-off  slope  than  those  of 
the  seisalc  sensors-  Increasing  noise  levels  between  the  7th  and  27th 
hours,  particularly  evident  in  the  vertical  coaponenc,  correspond  to 
an  increase  in  wind  speed  and  sea  state;  conditions  then  reverted  to 
those  of  the  Initial  part  of  the  measurecent  period. 

In  the  higher  frequency  region  (30  to  100  Hz),  shipping  noise  is 
demonstrated  by  sporadic  peaks  in  the  contours.  The  peaks  normally 
correspond  to  the  activities  of  aedlua-size  trawlers  and  fishing  boats 
passing  through  the  measurement  site. 
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Fig.  7  a 

Particle  nations  under  sea-state  0  and  6. 
The  hodographs  are  displayed  for  difterent 
frequency  kinds.  Soto,  for  sea-state  6,  x , 
y  and  z  scales  are  exaggerated  6  times  as 
compared  to  those  for  sea -state  0. 


Fig.  7b 

Azimuthal  power  pattern 
foi  two  different  fre¬ 
quency  bands. 


For  all  sensors  we  observed  a  frequency-dependent  increase  of  spectrum 
levels  with  increasing  sea  state  (see  selection  of  different  sea  states  in 
Fig.  4).  At  low  frequencies,  the  increase  is  greater  than  at  the  high 
frequencies.  Sea-state  effects  normally  disappear  above  the  60  to  70  Hz 
band.  Within  the  5  to  15  Hz  baud  peaks  correspond  to  the  sea  surface 
waves.  These  peaks  increase  with  increasing  sea  state  while  their  fre¬ 
quency  range  narrows.  Figure  7a  displays  two  sets  of  hodographs  (for  the 
xy  and  xz  planes)  obtained  at  sea-states  0  and  6.  The  first  column 
illustrates  tho  particle  motions  within  the  1  to  100  Hz  frequency  band 
while  the  remainder  of  the  columns  are  devoted  to  20  Hz  bandwidths 
covering  the  100  Hz  measurement  band.  The  energy  is  evidently  con¬ 
centrated  in  the  1  to  20  Hz  bandwidth,  for  which  particle  motion  is  evi¬ 
dent  in  both  vertical  and  horizontal  planes,  characteristic  of  interface 
waves  [15,16].  Figure  7b  shows  the  azimuthal  power  patterns  for  sea  state 
0  and  6  of  the  same  event.  The  patterns  of  sea-state  were  analysed  with 
both  large  (Af  ■  5)  and  narrow  (A f  ■  0.5)  frequency  bands.  The  direc¬ 
tional  patterns  indicate  the  wind  direction  for  which  the  sea  surface 
waves  had  the  same  orientation. 


Effects  of  Shipping 

The  high  frequency  reg'on  (30  to  100  Kz)  is  predominantly  influenced 
by  ship  traffic  noise.  lie  dominant  shipping  of  the  area  was  coastal 
trawlers,  fishing  boats  and  pleasure  craft  (nearby  shipping  activity) 
and  distant  commercial  shipping  passing  through  deep  water  to  or  from 
Genova  .at  distances  between  20  and  130  n.mi  from  the  station. 

Figure  8  illustrates,  for  different  sensors,  the  superposition  of  105 
power  spectra,  that  were  recorded  daily  for  the  study.  To  simplify  this 
large  amount  of  data  we  applied  basic  statistical  analysis  techniques  [17] 
to  obtain  average  power  and  mean  power  levels  for  each  sensor. 
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This  difference  (low  levels  on  the  hydrophone  spectra)  is  probably  due 
to  the  cut-off  of  the  waterborne  propagating  noise.  For  the  water  depth 
and  the  sediment  type  of  the  area,  the  waterborne  acoustic  propagation 
cut-off  frequency  corresponds  to  10-15  Hz  [18].  Apparently  the  energy  is 
coupled  into  the  noise  field  from  distant  shipping  (deep  water),  and  con¬ 
fined  to  the  continental  shelf  because  of  the  waterborne  acoustic  fre¬ 
quency  cut-off  process.  In  this  case,  the  energy  lost  within  the  water 
column  couples  into  the  seismic  waves  which  are  detected  by  the  geophones. 
This  phenomenon  appears  as  high  energy  levels  within  the  5  to  25  Hz  fre¬ 
quency  band  on  the  geophone  spectra.  This  observation  and  the  relatively 
lower  levels  observed  on  the  vertical  geophones  as  compared  with  the  hori¬ 
zontal  ones,  indicate  that  the  noise  field  is  mainly  related  to  interface 
and  Love  waves. 


COMPARISON  WITH  OTHER  DATA 

Experimental  infrasonic  ambient  noise  data  have  been  reported  by 
others.  The  data  set  for  sea  states  0  and  6  presented  in  this  paper  are 
superimposed  (Fig.  10)  on  Kibblewhite  and  Ewans  [13]  data  which  summarize 
various  other  ambient  noise  data  sets.  The  Mediterranean  continental 
shelf,  lacking  the  long  swells  of  the  open  oceans,  yields  lower  pressure 
amplitudes  at  low  frequencies.  The  same  data  set  compared  to  Wenz  curves 
[1]  within  the  10  to  100  Hz  band  fits  into  the  shallow  water  region  of 
usual  traffic  noise  (Fig.  11).  In  the  1  to  10  Hz  frequency  band,  spectrum 
levels  fit  nicely  into  the  seismic  background  and  sea  surface  wave  noise. 
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Fig.  11  Comparison  at  data  set  for  sea-states  0  and  6  with  I Yens  curves  111. 


CONCLUSION 

Systematic  noise  data  collected  over  the  1  to  100  Hz  frequency  hand 
for  a  long  period  revealed  that  noise  levels  correlated  well  with  wind-sea 
conditions  and  nearby  shipping  activity. 

Meteorological  noise  overlaps  with  shipping  noise  in  the  transient 
region  of  5  to  25  Hz  band.  In  the  same  band,  the  hydrophone  (water-borne 
noise)  revealed  very  low  energy  levels  as  compared  to  the  geophones.  This 
effect  is  most  probably  caused  by  the  water-borne  acoustic  frequency  cut¬ 
off  due  to  the  decreasing  water  depth. 

The  collection  of  data  for  a  longer  period  has  produced  better  noise 
statistics  for  the  area.  The  data  fit  well  into  previous  ambient  noise 
data  sets.  There  is  strong  evidence  that  at  infrasonic  frequencies  noise 
energy  propagates  essentially  as  interface  waves  and  ducted  Love  waves. 
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ABSTRACT 

The  Importance  of  the  geoacoustic  properties  of  marine  sediments  for 
low-frequency  shallow  water  propagation  is  well  established.  Reliable 
compressional  properties  can  be  obtained  from  laboratory  analysis  of 
collected  samples.  The  shear  properties,  however,  are  highly  affected  by 
the  deterioration  of  the  chemical  and  mechanical  bindings  caused  by 
de-pressurization  and  change  of  temperature  in  the  core  sampling  process. 
These  properties  therefore  have  to  be  determined  in  situ.  Isolated  shear 
waves  are  not  easily  generated  due  to  partial  coupling  into  compressional 
waves  at  all  interfaces.  In  the  case  of  stratified,  isotropic  media, 
however,  coupling  can  be  avoided  if  the  polarization  of  the  shear  waves  is 
parallel  to  the  interfaces.  Here  an  experimental  technique  is  described 
for  generation  and  detection  of  horizontally  polarized  shear  waves  -  or 
Love  waves  -  in  a  horizontally  stratified  sea-bed.  Shear  velocity  profiles 
are  determined  indirectly  by  numerical  modelling  of  the  experimental  data. 
Results  are  presented  for  different  types  of  marine  sediments. 


INTRODUCTION 

The  importance  of  the  geoacoustic  properties  of  marine  sediments  for 
low-frequency  acoustic  propagation  is  well  established  [1].  Geoacoustic 
properties  also  provide  a  major  source  of  information  for  the  determination 
of  the  static  and  dynamic  characteristics  of  the  sea-bed. 

When  seismic  waves  propagate  through  the  marine  sediments,  the  dynamic 
properties  of  the  sediments  determine  the  speed  and  the  attenuation.  An 
investigation  of  the  seismic  propagation  characteristics  can  therefore  pro¬ 
vide  accurate  information  on  the  geotechnical  characteristics  of  the  sea¬ 
bed. 


Assuming  that  the  sea-floor  sediments  behave  like  elastic  media,  the 
seismic  propagation  is  controlled  by  the  bulk  modulus  K,  the  shear  modulus 
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G,  and  the  density  p.  These  parameters  are  related  to  the  compressional 
wave  velocity,  Cp,  and  the  shear  wave  velocity,  Cs,  by 

Cp  =  [(K  +  4G/3)/p)1/2  and 

Cs  -  (G/p)1/2 

If  the  compressional  wave  velocity  (Cp)  and  the  shear  wave  velocity  (Cs) 
of  an  idealized  marine  sediment  can  be  measured  concurrently  with  a 
measurement  of  bulk  density,  the  elastic  parameters  of  the  sediments  can  be 
completely  characterized  (see  Tables  1  and  2). 


tabic  1.  Table  2. 


BASIC  WAVE  TYPES  AND  VELOCITIES 

1 

ELASTIC  PARAMETERS  IN  TERMS  OF  WAVE 
VELOCITIES  AND  BULK  DENSITY  (p) 

a)  BODY  WAVES 

K:  BULK  MODULUS  (incompressibility) 

K=p(Cj-4C]/3) 

-COMPRESSIONAL 

(longitudinal) 

E:  YOUNGs  MODULUS 

C,=[(K+4G/3)/p]'/1 

E=  2C[p(1+o)  or 

-SHEAR 

(transverse) 

E=pCi(l+o)(1-2o)/(l-o) 

C.=(G //>)" 

17. POISSON'S  RATIO  (transv./long.  strain) 

b)  DUCTED  WAVES 

or  1/2(Ci-2C])/(CJ-C|) 

G:SHEAR  MODULUS  (rigidity) 

~  L0VE  Ct=(G,/p)vl 

G=pC| 

c)  SURFACE  WAVES 

A= LAME'S  CONSTANT 

-  SCHOLTE  ,  , 

C„~(G/p)w 

X=p(Cj-2C5) 

Several  laboratory  and  field  methods  are  available  to  measure  density 
and  wave  velocities  in  bottom  sediment  cores.  However,  the  reliability  of 
such  measurements  is  degraded  by  sample  disturbance  and  temperature  and 
pressure  changes.  In  particular,  the  shear  properties  are  highly  effected 
’•y  the  deterioration  of  the  chemical  and  mechanical  bindings  caused  by  the 
differences  in  temperature  and  pressure  between  the  sampling  and  the 
laboratory  measurements.  Also,  in  the  laboratory,  the  small  size  of  the 
sample  does  not  allow  direct  information  to  be  obtained  about  the  response 
of  these  materials  to  low-frequency  seismic/acoustic  waves.  Therefore  the 
most  reliable  geoacoustic  measurement  techniques  for  ocean-bottom  sediments 
are  therefore  in  situ  techniques. 

The  choice  of  an  actual  experimental  technique  is  etrongly  influenced 
by  the  structure  of  the  sea  bed.  If  the  bottom  was  homogeneous  to  a 
significant  depth,  the  wave  speeds  and  attenuation  could  be  measured 
directly  because  different  seismic  wave  types  could  be  separated.  In 
reality,  however,  recent  sediments  consist  of  a  relatively  thin  top  layer 
of  unconsolidated  clay  or  silty  clay  overlying  semi -consolidated  silts  and 
sands  [2,3].  This  structure  results  from  the  lowering  of  sea  levels 
during  the  glaciation  of  the  Pleistocene  epoch  in  which  sand  was  deposited 
over  wide  areas  of  the  continental  shelves.  The  shelves  were  subsequently 
covered  by  unconsolidated  sediments  during  sea  level  rising  during  the 
post-glacial  times.  Sedimentary  processes  such  as  sorting,  compaction  and 
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consolidation,  also  have  contributed  to  the  highly  stratified  structure  of 
the  sea  bed. 

In  such  a  layered  sea  bed  the  compressional  and  shear  waves  couple, 
yielding  an  often  complex  wave  field.  The  geoacoustic  properties  therefore 
must  be  determined  indirectly  (inversely)  from  the  propagation  charac¬ 
teristics.  The  most  important  attenuation  mechanism  of  low  frequency  pro¬ 
pagation  in  shallow  water  is  the  conversion  of  acoustic  energy  in  the  water 
column  into  shear  waves  In  the  bottom  [1].  Because  of  this  importance  this 
mechanism  has  received  special  attention  during  the  past  decade.  Because 
the  propagation  of  seismic  interface  waves  (also  known  as  Scholte  waves)  is 
almost  entirely  controlled  by  their  shear  properties,  significant  effort 
has  been  applied  to  the  analysis  of  these  properties  [4,5].  However,  the 
inversion  involves  significant  numerical  effort  not  easily  automated. 

Here  it  will  be  demonstrated  that  horizontally  polarized  shear  waves 
(SH),  also  called  Love  waves,  can  be  generated  by  simple  means  in  the  low- 
velocity  duct  formed  by  the  upper  sediments.  Because  the  shear  waves  do 
not  couple  into  compressional  waves  in  a  horizontal  stratification  of 
isotropic  layers,  their  propagation  characteristics  directly  reflect  the 
shear  properties.  The  Love  waves  therefore  form  an  alternative  to  the 
seismic  interface  waves  as  the  basis  of  the  shear  property  inversion. 
Further,  the  scalar  nature  of  the  wave  equation  for  Love  waves  signifi¬ 
cantly  decreases  the  calculations  involved  in  the  inversion  compared  to  the 
requirements  for  the  interface  waves. 


EXPERIMENTAL  PROGRAMME 

The  Love  wave  experiments  were  caried  out  as  part  of  a  multinational 
research  programme  to  measure  the  velocity  and  attenuation  of  seismic  waves 
and  supporting  geotecnnical  parameters  in  marine  sediments.  The  programme 
was  initiated  by  SACLANT  to  bring  together  experts  in  geophysics,  geology 
and  acoustics  to  compare  various  measurement  and  analysis  techniques  for 
data  gathered  under  Identical  environmental  conditions. 

For  this  program  four  institutes  played  a  leading  role  in  the  specific 
problems  related  to  sampling,  analysis  and  interpretation  of  ocean  bottom 
sediments.  The  institutes  and  their  tasks  were: 

-  Institute  of  Marine  Geology  (IGM-CNR),  Bologna,  Italy;  sampling  the 
sea-bed  using  different  techniques,  geotechnical  analysis  and  interpreta¬ 
tion,  and  geological  description; 

-  Centro  de  Flslca  Applicada  (CENFA),  Madrid,  Spain;  measuring  reflec¬ 
tion  loss  under  controlled  conditions  in  a  layered  sea-bed  with  variable 
sensitivity  gradients,  and  testing  a  model  developed  by  the  CENFA  s^'an- 
tists; 


US  Naval  Ocean  Research  and  Development  Activity  (NORDA),  Bay  St. 
Louis,  MS,  USA;  sediment  sampling,  laboratory  analysis  cf  the  samples  for 
compressional  and  shear  wave  velocities,  sediment  mass  property  analysis, 
measurements  with  the  ELAC  Sea-Bed  Classifier  to  define  high-resolution 
layering,  and  sediments  classification; 

SACLANT  ASW  Research  Centre  (SACLANTCEN),  La  Spezla,  Italy; 
programme  coordination,  measurement  of  seismic  wave  velocities  and  atte¬ 
nuation  using  ocean-bottom  seismometers,  sampling  and  geoacoustic  analysis 
of  sediments  and  site  surveying  (using  side-scan  sonar  and  seismic 
profilers);  SACLANTCEN  theoretical  propagation  models  were  used  with 
seismic  data  to  determine  the  sea-bed  properties. 
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MEASUREMENT  TECHNIQUES 

The  techniques  developed  by  the  different  Institutions  were  applied  to 
four  different  sedimentary  areas  In  the  vicinity  of  La  Spezia,  Italy 
(Fig.  1) 

Shear-wave  velocity  measurements 

A  directional  energy  source  was  developed  to  excite  both  interface  and 
ducted  Love  waves.  Figure  2  illustrates  the  techniques  used  by  SACLANTCEN. 
Five  ocean  bottom  seismometers  (OBS),  each  with  tri-axlal  geophones  and  a 
hydrophone,  were  deployed  along  the  measurement  range  in  each  area.  For 
the  short  ranges  applied  (<25  m) ,  a  high  energy  source  level  was  not 
needed.  The  signals  were  generated  by  shooting  a  spear  from  a  hydraulic 
gun  mounted  on  a  frame  into  a  mass  resting  on  the  sea-bed  (Fig.  3).  This 
type  of  source  generates  a  signal  below  the  elastic  limits  of  the  sedi¬ 
ments,  thus  avoiding  changes  in  sediment  properties.  The  direction  of  the 
gun-mount  could  be  changed  as  appropriate  to  generate  different  wave  types. 
The  experiment  was  conducted  by  diving  scientists  [6]. 

Figure  4  shows  examples  of  the  signals  received  at  a  distance  of  25  m 
from  the  source.  For  example,  in  Fig.  4a,  the  geophone  received  maximum 
energy  from  a  transverse  (shear)  pulse  traveling  in  the  x-dlrection.  The 
alignment  of  the  corresponding  particle-motion  diagram  (hodograph)  also 


Fig.  1  Expcritacnt  areas,  La  Spozia ,  Italy. 
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Fig.  2  Techniques  used  to  measure  Love  wave  velocity 


Fig.  3 

SACLANTCEN  technique  to  generate 
seismic  signals,-  a  gun  shoots  a 
spear  against  a  mass  coupled  to 
the  sea-bed  with  a  heavy  metal 
plate. 


Fig.  4  Acoustic  signals  received  at  a  distance  of  2S  m  from  the  source, 
and  corresponding  particle  motion  diagrams  (hodographs) . 
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confirms  that  the  signal  is  confined  to  the  horizontal  plane  (i.—axis,. 
This  dispersed  form  of  the  transverse  signal  indicates  that  this  is  no  a 
pure  body  wave  but  a  Love  wave  propagating  within  the  upper  sedln.eni.8ry 
layer.  This  assumption  was  confirmed  by  synthetic  seismograms  produce--’,  by 
the  numerical  model  SAFARI  (7],  The  radial  (Fig.  Ab)  and  vertical 
(Fig.  Ac)  signals  were  much  more  complex;  the  particle  motions  were  con¬ 
fined  to  either  the  radial  or  vertical  plane  with  very  small  transverse 
deflections,  a  characteristic  behaviour  of  interface  waves  [A]. 

Bottom  reflectivity  measurements 

Figure  5  illustrates  the  system  used  by  NORDA  and  CENFA  to  measure  the 
acoustic  reflection  char.-  teristics  of  the  bottom.  NORDA's  ELAC  Sea-Bed 
Classifier  (30  kHz)  was  used  both  to  obtain  high  resolution  sub-bottom 
layering  information  and  high-resolution  bottom  reflection  loss  measure¬ 
ments.  SACLANTCEN  also  used  an  EG&G  UNIBOOM  seismic  source  for  bottom 
reflection  measurements. 

Sediment  sampling 

Figure  6  shows  the  techniques  used  to  sample  sea  bed  sediments.  Samples 
were  collected  by  the  diver-scientists  who  used  a  grab  sampler,  a  box 
corer,  a  gravity  corer,  and  IGM's  newly  developed  compresslonal  air-powered 
corer . 

Side-scan  sonar 


The  site  surveys,  using  side-scan  sonar  (EG&G  MARK1B)  provided  knowledge 
of  both  the  surface  sediment  distribution  and  the  sea  bed  topography. 

Continuous  seismic  profiling 

An  EG&G  UNIBOOM  source/receiver  and  a  streamer  were  towed  behind  the 
survey  ship.  Resolution  was  poor  in  water  depths  less  than  10  m. 


Fig.  5  JtOPDA  a nd  CENFA  sgstcns  to  treasure  the  acoustic  reflection 
characteristics  of  the  bottom. 
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Fig.  6  Techniques  used  to  sample  sea-bed  sediments 

a)  hand  cotcr 

b)  gravi tg  corer 

c)  ZCH's  compressional  air-powered  corer 


Laboratory  analysis 

The  cores  were  analyzed  every  2  cm  along  their  length  for  the  following 
properties  and  characteristics; 

-  compressional  wave  velocity 

-  shear  wave  velocity 

-  porosity 

-  water  content 

-  wet  density 

-  dry  density 

-  torvane  and  shear  strength 

-  Atter.berg  limits 

-  grain  size 

-  mineralogic  composition 


EXPERIMENTAL  RESULTS 

Data  from  two  areas,  1  and  2,  with  different  bottom  conditions  are  pre¬ 
sented  here  (see  Fig.  1). 

Area  A  is  situated  along  the  Cinque  Terre  coast  of  the  Ligurian  Sea; 
Area  B  is  situated  at  the  eastern  entrance  to  the  Portovenere  channel  of 
the  Gulf  of  La  Spezla. 


SACIWCEN  CP-37 


2-6 


AKAi.  et  al :  Love  waves 


Area  1 

This  study  area  had  a  water  depth  of  IS  m.  The  sea-bed  consisted  of 
coarse  sand  and  fine  gravel-  Slde-3can  sonar  and  observations  by  diving 
scientists  Indicated  that  the  area  was  covered  with  ripples  (approximately 
X  =  2  m)  created  by  wave  action.  To  avoid  distorting  the  propagation  path 
the  axis  of  the  measurement  range  was  laid  out  parallel  to  the  axes  of  the 
ripples.  Seismic  profiles  (Fig.  7)  show  a  sub-bottom  reflector  approxima¬ 
tely  2  ms  from  the  bottom.  The  bottom  sediments  apparently  were  derived 
from  land  slides  and  the  reworking  of  sediments  along  the  steep  coastline. 

Figure  8  summarizes  the  measured  parameters  along  a  core  sample  taken  in 
this  area.  The  high  percentage  of  gravel  in  the  sea  bed  limits  the  core 
penetration  to  less  than  20  cm.  This  gravel  consists  mainly  of  limestone, 
slate,  and  sandstone,  apparently  derived  from  the  adjacent  mountains  which 
have  a  similar  lithology. 


PUNTA  MONASTEROU  AREA  1 
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Fig.  8  Area  1.  Measuzod  parameters  along  a  core  sample. 
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Figure  9  is  a  time  vs  distance  plot  of  the  output  of  the  transverse 
(shear)  wave  sensitive  geophone.  Distinct  wavelets  can  be  observed  tra- 
veling  v'th  speeds  up  to  90  m/s.  The  inversion  was  performed  by  means  of 
the  SAFARI  model  [7],  using  trial  and  error.  The  very  simple  3hear  velo¬ 
city  profile  shown  in  Fig.  10a  yielded  the  synthetic  seismograms  shown  in 
Fig.  10b,  which  agree  well  with  Fig.  9. 


M 


Fig.  9  Area  1.  Time  vs  distance  plot  of  the  shear  wave 
sensitive  geophone. 


Area  2 


Side-scan  sonar  records  of  the  eastern  entrance  fo  the  Portovenere 
Channel  indicate  a  flat  bottom.  Sub-bottom  profiles  taken  along  the 
measuremen:  base  indicate  four  distinct  reflectors  (Fig.  11). 

Figure  12  illustrates  the  physical  parameters  measured  along  one  of  the 
cores  taken  in  Area  2.  The  material  consists  of  silty  clays  and  clayey 
silts  with  very  small  variations.  The  mineralogy  of  these  sediments  con¬ 
sists  mostly  of  quartz  and  calcite  with  some  magnesium.  Of  the  clay 
minerals,  60%  is  illite,  10%  chlorite  and  30%  smectite  and  kaolinite.  The 
percentage  of  illite  increases  towards  the  bottom  of  the  core  sample  with 
the  depth  in  the  core  (and  hence  the  sea  bottom).  Attenberg  limits 
measured  were  51  to  55%  for  the  liquid  limit  (W^)  and  about  25%  for  the 
plastic  limit  (Wp).  The  liquid  limit  tends  to  increase  with  depth, 
verifying  that  the  sediments  become  finer  towards  the  bottom  of  the  core. 
Figure  13  shows  the  shear  waves  detected  by  transverse  wave  sensitive 
geophones.  A  distinct  wavelet,  traveling  at  30  m/s,  indicates  hove  wave 
speeds  in  the  upper  sedimentary  layers  about  1/3  of  that  in  Area  l. 
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Fig.  13 


Area  2. 
geophones 


Tine  vs  distance  plots  of  the 
•  and  shear-wave  speed. 


shear  wave  sensitive 


CONCLUSIONS 

™<  its:  srss,  ss  £rs*sL that  —  —  -  «* 

sediment  layers.  The  Love  waves  are  horizon^*  i*?  L’PPerinoat  unconsolidated 
horizontally  atr.tiflod  JZ!,  “rl*1!’ "»'i  1» 
types.  This  property  greatly  slmnHf<»«  . .  convert  Into  other  wave 
behaviour.  l.o"e  L  J,  Sy  »L.for.  „r.  »J«Uln«  of  their 

•fnatio.  of  the  eedi.ent  rteer  pro^m"  l°* 

based  on  the  properties  of  the  more  complex  seiamir  <  ye*rs  been 
experiment  demonstrated  that  two  wave  types  fstho^  in*e»;f«ce  wave,.  The 
be  generated  and  detected  in  a  single  exp«imfnJ«l  ?"d  Lov“>  °*n  ln  fact 
therefore  is  a  realistic  possibility  JK  ?«  1  8^up'  The  combination 
prediction  accuracy  of  deduced  shear  sneer/  sPPHed,  can  improve  the 

blnation  oa.  olnru.tn  th.  dojroo  *f  ."l.otro.v  if.V  F“"h=-  tte 
layers.  6  1  anisotropy  in  the  upper  sedimentary 
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MEDIUM-INDUCED  LOW-FREQUENCY  FLUCTUATIONS  IN  ACOUSTIC  TRANSMISSION  LOSS: 
EXAMPLES  FROM  MEASUREMENTS  IN  SELECTED  GEOGRAPHICAL  AREAS 
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ABSTRACT 

The  ocean  is  a  complex,  highly  variable  acoustic  medium.  A  propa¬ 
gating  acoustic  signal  is  affected  by  a  host  of  phenomena,  Including  the 
sea-surface  and  bottom,  volume  inhomogeneities,  internal  waves  and  tides, 
and  non-stationary  water  masses.  These  effects  cause  fluctuations  in  the 
amplitude  and  phase  of  an  acoustic  signal  and  an  accompanying  loss  in  its 
coherence  properties.  The  responsible  mechanisms,  and  hence  the  acoustic 
effects,  cover  a  wide  range  of  temporal  and  spatial  scales  and,  in 
general,  can  be  understood  only  in  terms  of  deterministic  and  random  for¬ 
ces  acting  in  concert.  Although  the  listing  of  these  mechanisms  is 
generally  easier  than  their  isolation  in  realistic  situations,  it  has  been 
possible  to  correlate  fluctuations  in  acoustic  transmission  loss  with 
environmental  variability,  especially  for  very  low  acoustic  frequencies. 
In  particular,  measurements  conducted  by  SACLANTCEN  in  diverse  geographi¬ 
cal  areas  have  identified  semi-diurnal  effects  (tidal  as  well  as  heating), 
internal  waves,  inertial  oscillations,  and  moving  water  masses  as  signifi¬ 
cant  contiibutors  to  low-frequency  acoustic  fluctuations. 


INTRODUCTION 

The  basic  problem  in  using  acoustics  in  the  ocean  Is  the  complexity  of 
the  medium.  The  parameters  controlling  the  propagation  vary,  usually 
unpredictably,  both  spatially,  and,  more  significantly,  temporally. 
Consequently,  an  acoustic  signal  propagating  in  such  a  medium  is  affected 
not  only  by  interactions  with  the  bottom  and  surface  boundaries,  but  also 
by  volume  inhomogeneities  caused  by  non-uniformities  in  temperature,  den¬ 
sity,  and  salinity  distributions.  The  degradation  of  the  acoustical 
signal  is  manifested  by  fluctuations  in  its  amplitude  and  phase  and  by  an 
.ccorapanylng  loss  in  its  coherence  properties. 

Generally,  this  degradation  results  in  a  decrease  in  performance  of 
underwater  acoustic  systems,  including  sonar.  In  particular,  the  estima¬ 
tions  of  range  and  angular  bearing  of  a  sonar  system  can  be  adversely 
affected  by  temporal  and  spatial  coherence  losses.  For  this,  and  other 
reasons,  an  understanding  of  the  causes  of  signal  degradation  is  essential. 


SACLANTCEN  CP-37 


3-1 


ALI  et  el :  Acoustic  transmission  loss  fluctuations 

This  paper  will  consider  environmentally-induced  signal  fluctuations  as 
they  affect  transmission  loss. 

Although  the  mechanisms  leading  to  fluctuations  in  acoustic  propaga¬ 
tion  are  diverse,  an  essential  common  feature  is  an  associated  non¬ 
uniformity  in  the  medium,  either  temporal  or  spatial  or  both.  Depending 
on  the  temporal  and  spatial  scales  involved,  the  mechanisms  can  be  con¬ 
sidered  either  deterministic  or  random  [1].  The  general  circulation  of 
the  ocean  ("ocean  climate")  and  its  associated  current  systems  (Gulf 
Stream,  Kuroshlo,  etc.)  are  characterized  by  horizontal  scales  of  variabi¬ 
lity  limited  only  by  the  size  of  the  basin,  vertical  scales  of  a  few 
100  m,  and  temporal  scales  from  a  few  days  to  seasonal.  These  are  deter¬ 
ministic  structures.  The  intermediate  scales  of  variability,  including 
ocean  motions  such  as  fronts  and  eddies,  can  also  be  considered  to  be 
deterministic  perturbations  from  the  mean  structure.  The  associated  sca¬ 
les  of  variabilities  are  of  the  order  of  100  to  1000  km  in  the  horizontal, 
to  ocean  depth  in  the  vertical,  and  days  to  months  in  time. 

Internal  waves,  fine-structure,  and  microstructure  comprise  the 
smaller  scales  and  are  considered  random  features.  The  internal  waves  are 
characterized  by  scales  from  100  m  to  10  km  or  more  in  the  horizontal,  1 
to  100  m  in  the  vertical,  and  from  about  10  min  to  1  day  in  time.  Because 
they  owe  their  existence  to  the  restoring  forces  due  to  the  density  gra¬ 
dient  and  the  Coriolis  force,  the  frequency  spectra  of  internal  waves  are 
bounded  by  .  the  inertial  frequency  at  the  low  end  and  by  the  buoyancy  fre¬ 
quency  (Brunt-V  aisala)  at  the  high  end.  Variability  induced  by  internal 
waves  has  been  found  to  be  a  very  significant  source  of  sound  scattering 
and  has  received  considerable  attention  in  recent  years  (2,1,4].  In 
essence,  this  increased  activity  arises  from  two  complementary  fa-ts:  on 
the  one  hand  the  spatial  scales  of  internal  waves  match  the  acousti.  wave¬ 
length  over  a  broad  frequency  range,  and  thereby  affect  the  acoustic 
field;'  on  the  other  hand  the  internal  wave  is  one  of  the  few  ocean 
features  for  which  a  reasonably  effective  statistical  model  (Garrett-Munk) 
is  available.  Fine-  and  microstructure  variability  involve  horizontal 
scales  from  several  metres  to  hundreds  of  metres,  vertical  scales  from 
centimetres  to  about  10  m,  and  temporal  scales  of  the  order  of  millise¬ 
conds.  Such  variability  generally  affects  sound  propagation  in  the  fre¬ 
quency  range  from  approximately  one  to  tens  of  kilohertz. 

Table  1  [1]  summarizes  the  temporal  scales  of  environmentally-induced 
fluctuations  often  observed  in  acoustic  propagation  experiments.  Of  par¬ 
ticular  relevance  to  low  frequency  acoustic  propagation  Is  the  information 
contained  within  the  dotted  area  of  Table  1.  The  remainder  of  the  paper 
is  devoted  to  a  discussion  of  three  SACLANTCEN  measurements,  the  results 
of  which  show  a  clear  correlation  between  transmission  loss  fluctuations 
and  the  environmental  phenomena  highlighted  here. 

In  particular,  long-period  internal  waves  (inertial  oscillations)  and 
semi-diurnal,  tide-related  effects  figure  prominently  in  the  measurements 
conducted  in  two  areas  of  the  North  Atlantic  Ocean  and  one  in  the 
Mediterranean.  Moreover,  although  tidal  forces  are  important  at  both 
sites,  differences  among  the  local  oceanographic  conditions  lead  to 
interesting  contrasts  between  the  nature  and  the  degree  to  which  these 
tidal  forces  affect  the  acoustic  propagation.  In  particular,  at  the 
Mediterranean  test  site  volume  effects  in  the  water  column  seemed  to  be  of 
greater  significance  than  in  the  North  Atlantic  test  site,  where  acoustic 
interaction  with  the  bottom  was  the  predominant  source  of  the  fluc¬ 
tuations.  Therefore,  discussions  of  the  results  will  be  grouped  under 
physical  mechanisms,  viz.,  under  "volume  effects"  and  "bottom  effects", 
rather  than  under  the  less  relevant  geographical  indicators. 
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Table  1  Temporal  variation  observed  in  acoustic  data 
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24  h 
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IS  day  -«*- 
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MESOSCALE  PHENOMENA 
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OCEAN  CLIMATE 
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PHASE  AND  AMPLITUDE  FLUCTUATIONS  OFTEN 
SIMPLY  CORRELATED  WITH  TIDAL  VARIATIONS. 


GREATEST  EFFECTS  IN  SHALLOW  WATER 

_ I 


FLUCTUATIONS  IN  PROPAGATION  PATH,  ACOUSTIC 
INTENSITY;  DISTORTION  OF  SIGNAL  SHAPE; 

SLOW  MODULATION  OF  ACOUSTIC  RAYS 


FLUCTUATION  IN  AMPLITUDE,  PHASE  AND 
PROPAGATION  PATH 


1  MEASUREMENTS  SHOWING  VOLUME  EFFECTS 

1.1  Propagation  Characteristics  of  the  Test  Environment 

The  acoustic  propagation  experiment  was  conducted  in  a  shallow  water 
region  of  the  Mediterranean  where  the  water  depth  varied  from  about  40  m 
to  85  m,  as  seen  in  Fig.  1.  Although  the  bathymetry  of  the  area  is  fairly 
complex,  that  along  the  propagation  run  is  relatively  simple.  The  water 
circulation  in  the  region  consists  of  a  three-layer  system:  an  eastward¬ 
flowing  surface  layer  of  Atlantic  water,  an  intermediate  layer  in  which 
turbulent  mixing  occurs,  and  a  westward-flowing  bottom  layer  of  more 
saline  (Levantine)  water.  A  temperature/salinity  profile  of  the  measured 
data  confirms  this  general  behaviour. 

For  the  experimental  situation  depicted  in  Fig.  1,  broad-band 
(explosive)  sources  were  dropped  at  approximately  hourly  intervals.  A 
vertical  array  of  hydrophones,  35  km  distant,  received  the  signals. 
Simultaneously,  samplings  were  taken  of  the  pertinent  oceanographic  para¬ 
meters:  sound  speed,  temperature,  salinity  and  density  (STEV  casts). 
The  test  was  conducted  during  summer  conditions.  Figure  2  presents  typi¬ 
cal  depth  profiles  of  the  environmental  parameters  in  this  area.  This 
type  of  profile  tends  to  favour  downward  acoustic  refraction  and,  there¬ 
fore,  results  in  greater  bottom  interaction  than  »ould  occur  in  winter. 
Vertical  stratification  and  some  fine-structure  are  evident  In  the  profi¬ 
les.  A  closer  analysis  of  the  sound  speed  profile  over  shorter  intervals 
in  depth  and  sound  speed  reveals  more  clearly  the  presence  of  fine- 
structure,  characterized  by  vertical  dimensions  of  the  order  of  cen¬ 
timetres  to  one  or  two  metres.  As  already  noted,  such  fine-structure  will 
tend  to  scatter  sound  of  frequency  ranges  higher  than  of  interest  herein. 
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Fig.  2 

Depth  profiles  of  temperature,  sound  speed, 
and  salinity  (Mediterranean  site). 
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The  effect  on  acoustic  propagation  is  shown  in  Fig.  3  <5>,  which  pre¬ 
sents  contours  of  measured  transmission  loss,  in  1/3  octave  bands,  in  the 
f requency/range  plane.  The  existence  of  an  optimum  frequency  range  (i.e. 
minimal  transmission  loss)  for  acoustic  propagation  is  clearly  evident. 
In  this  case,  it  lies  between  approximately  100  and  400  Hz.  The  explana¬ 
tion  for  this  is  that  acoustic  energy  at  very  low  frequencies  suffers 
large  attenuation  as  a  result  of  bottom  penetration  which  increases 
directly  with  wavelength.  On  the  other  hand  the  energy  at  the  very  high 
frequencies  is  greatly  attenuated  by  absorption  in  the  water  column  and, 
possibly,  by  scattering  from  fine-structure.  Hence  an  optimum  frequency 
range  exists  somewhere  between  the  two  extremes.  In  other  words,  for  the 
conditions  typefled  by  Fig.  3,  the  ocean  behaves  like  a  band-pass  filter 
for  propagating  broadband  acoustic  signals. 


1.2  Temporal  Variability  in  the  Environmental  Data 

Of  primary  interest  for  our  purpose  is  the  temporal  variation  of  the 
relevant  parameters  at  a  fixed  range,  here  35  km.  Figure  4  shows  an 
example  of  the  variation  of  sound  speed  with  depth  during  a  period  of  25 
hours. 

The  depth  range  from  approximately  25  to  35  m  comprises  the  steepest 
portion  of  the  thermocline.  The  region  down  to  20  m  or  so  Is  essentially 
Isovelocity  with  a  sound  speed  of  approximately  1538  m/s.  The  fluc¬ 
tuations  in  sound  speed  are  quite  evident,  particularly  at  a  depth  of  25  m 
or  so  within  the  thermocline.  The  plots  clearly  indicate  an  oscillation 
in  the  width  of  the  mixed  layer  (surface  duct).  Figure  4b,  In  particular, 
indicates  that  typical  features,  e.g.,  "kinks"  in  the  curve,  appear  to  be 
migrating  vertically  with  time  -  a  behaviour  often  suggesting  the  presence 
of  internal  waves  [Gregg,  6J.  The  frequency  content  of  these  oscillations 
is  of  particular  interest  because  it  provides  clues  to  the  responsible 
mechanisms.  Examples  of  frequency  spectra  of  the  relevant  environmental 
parameters,  obtained  from  Fast  Fourier  Transforms  of  the  corresponding 
normalized  time  series  at  25  m  depth,  are  shown  In  Fig.  5. 
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Fig.  4  Temporal  variation  of  sound  speed  over  approximately  25  hours. 
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Fig •  5  Spcct c<t  of  sound  speed  and  temperature  at  25  n. 

The  dominant  fluctuations  occur  In  the  frequency  range  from  0.08  to  0.06 
cycle/h,  or  for  periods  from  12  to  17  hours.  This  is  a  broad  response  and 
includes  both  the  semi-diurnal  period  and  Inertial  period  (approximately 
17  hours  here).  The  time  series  was  too  short  to  allow  for  the  spectral 
separation  of  the  two  responses.  From  the  results  obtained  at  other 
depths  it  does  appear,  however,  that  the  relative  significance  of  the 
semi-diurnal  tidal  and  inertial  effects  depends  on  depth.  In  particular, 
from  the  analysis  of  the  depth-dependent  frequency  spectra  it  appears  that' 
there  are  fluctuations  of  inertial  period  in  the  surface  layers,  and  that 
as  the  depth  Increases  the  dominant  fluctuations  are  semi-diurnal. 

Inertial  oscillations  may  occur  In  connection  with  a  sudden  change  of 
the  wind  (i.e.  a  wind  Impulse  of  short  duration)  and  changes  in  barometric 
pressure  (71.  The  whole  of  the  Mediterranean  is  dominated  by  meteorologi¬ 
cal  forcing;  therefore  it  is  likely  that  inertial  oscillations  are  present 
in  the  study  area.  That  they  dominate  the  surface  waters  is  probably 
related  to  their  meteorological  origin.  The  local  water  stratification 
retains  most  of  this  energy  in  the  surface  layers,  so  that  at  depth  the 
semi-diurnal  components  are  discernible. 

In  the  spectra  of  Fig.  5  there  also  appear  to  be  higher-frequency 
phenomena,  with  periodicities  of  about  2,4,5,  and  7  h.  These  may  be  asso¬ 
ciated  with  internal  wave  activity.  Note  that  the  drop-off  rate  (sloping 
line)  of  the  sound  speed  spectrum  Is  approximately  ,  3  slope  consistent 
with  that  predicted  by  the  Garretl-Munk  model  of  Internal  oaves. 


1.3  Fluctuations  in  Acoustic  Transmission  Loss 

Figure  6  shows  the  contours  of  measured  transmission  loss  (1/3  octave 
bands)  in  the  frequency/time  plane  for  source  and  receiver  depths  of  50  m 
and  40  m,  respectively.  The  higher  frequencies  (  >1.6  kHz)  exhibit  fat 
more  pronounced  fluctuations  than  the  lower  frequencies.  This  may  indi¬ 
cate  that  the  responsible  environmental  phenomena  have  physical  dimensions 
that  are  comparable  to  the  acoustic  wavelengths  of  the  higher  frequencies. 
In  order  to  demonstrate  this  selective  frequency  effect  more  clearly,  the 
spectra  of  transmission  loss  at  200  Hz  and  1600  Hz  were  compared 
(Fig.  7). 
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fig.  6  Contours  of  neasurctl  transnission  loss  fluctuations  at  a 
fixed  range  (35  ko). 
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Fig.  7  Spec era  of  trensoissicn  loss  foe  Cvo  selected  freg-encics. 


The  spectra,  obtained  for  the  same  source/receivor  as  those  In 
Fig.  6,  emphasize  the  difference  in  the  effect  of  variability.  These 
results  suggest  that  the  optimum  frequency  range  is  less  susceptible  to 
environmental  variability  than  other  frequency  ranges.  A  comparison  of 
these  spectra  with  those  shorn  in  Fig.  5  reveals  a  good  correlation  bet¬ 
ween  the  spectra  of  environmental  variability  and  the  higher  frequency 
transmission  loss  spectrum. 
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Fig.  8  Ce oaetry  of  the  north  Atlantic  aeasureaents. 


2  MEASUREMENTS  SHOWING  BOTTOM  EFFECTS 

2.1  Propagation  Characteristics  of  the  Test  Environment 

The  North  Atlantic  test  site  Is  a  very  complex  geological, 

oceanographic,  and  biological  province  and  has  been  the  object  of  many 

investigations  [8J.  The  test  site  Is  located  In  an  area  affected  by  a 

large-scale  permanent  front  resulting  from  the  convergence  of  two  water 
masses.  In  particular,  at  this  location,  cold,  less  saline  Arctic  water 
meets  warmer  Atlantic  water  to  form  the  Polar  oceanic  front.  The  substan¬ 
tial  differences  la  salinity  and  temperature  of  the  two  water-types 
(35%„  and  6*  to  7°C  for  the  Atlantic  vs  34, 6^,  and  about  0*C  for  the 

Arctic)  cause  the  Polar  front  to  be  characterized  by  steep  gradients  In 
these  parameters,  with  significant  Implications  for  sound  propagation.  An 
added  complication  arises  from  the  changing  position  of  the  front,  which 
osclllat*.-  with  semi-diurnal  tidal  periodicity. 

TI»e  geometry  and  relevant  details  of  the  measurements  are  provided  In 
Fig.  8.  The  two  tracks,  A-B  and  D-C,  represent  the  propagation  paths  for 
measurements  conducted  one  year  apart.  The  environmental  conditions  alorg 
the  two  paths  were  somewhat  different.  Although  both  tracks  A-B  and  IhC 
cross  the  Polar  front,  in  witer  changing  from  deep  to  shallow,  track  A-B 
is  over  a  generally  hard  bo. '.ora  (sand)  while  track  D-C  Is  over  a  generally 
soft  bottom  (sand,  silt,  clay).  From  the  point  of  view  of  temporal 
variability,  however,  a  more  fundamental  difference  between  the  t*o 
receiver  locations  arises  from  their  relative  proximity  to  the  Poltr 
front.  In  particular,  station  C  is  directly  affected  by  movements  of  tie 
front,  the  water  above  20  m  changing  from  Polar  to  Atlantic  with  nn 
approximately  12-h  periodicity,  while  the  waters  below  20  m  remain  of  tlie 
Polar  type.  In  comparison,  station  B  remains  In  virtually  isothermal 
Polar  water,  north  of  the  front,  as  shown  schematically  in  Fig.  8. 
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In  effect,  the  Polar  front  serves  as  a  demarcation  line  between  two 
different  propagation  areas.  To  the  north  of  the  front,  in  isothermal 
water,  an  important  part  of  the  propagation  will  be  in  shallow  water  under 
upward-refracting  conditions,  resulting  in  small  total  transmission  losses 
as  a  result  of  less  bottom-interaction.  On  the  other  hand,  in  waters 
affected  by  the  front,  propagation  will  be  under  strong  downward 
refracting  conditions,  as  a  result  of  the  steep  gradients,  and  thereby 
subject  to  higher  losses,  particularly  at  low  frequencies.  As  noted, 
measurement  site  B  was  always  in  isothermal  waters,  whereas  C  was  periodi¬ 
cally  subject  to  the  effects  of  the  front.  From  this  one  would  expect  the 
transmission  losses  measured  at  site  B  to  be  both  lower  in  magnitude  and 
subject  to  less  temporal  fluctuation  than  those  at  site  C.  Though 
generally  valid,  this  conclusion  is  nevertheless  based  on  an  over¬ 
simplification  of  the  propagation  conditions.  In  particular,  the  propaga¬ 
tion  depends  not  only  on  the  front,  but  also  on  frequency,  source  and 
receiver  depths,  bottom  conditions  (which  are  different  for  the  two 
tracks),  and,  possibly,  other  features  such  as  currents  and  inertial 
oscillations. 


2.2  Selected  Examples  of  Temporal  Variability 

In  the  following,  typical  results  of  measurements  made  over  a  48-h 
period  will  be  presented,  emphasizing  the  data  obtained  at  position  C. 
The  experimental  procedure  was  essentially  the  same  as  that  for  the 
Mediterranean  measurements,  described  previously. 

Figure  9  provides  a  typical  example  of  the  fluctuations  in  measured 
transmission  loss  at  site  C  for  several  frequencies,  along  with  the  pre¬ 
dicted  tidal  curve  for  the  area.  It  is  clear  that  the  fluctuations  in 
transmission  loss  correlate  well  with  tidal  periodicity,  particularly  for 
the  lower  frequencies.  This  is  consistent  with  the  assumption  that  the 
lower  frequencies  are  more  affected  by  the  tidally-induced  change  in 
downward-refracting  conditions.  Both  diurnal  and  semi-diurnal  periodici¬ 
ties  are  evident,  the  former  seemingly  the  more  significant.  Further, 


16  JULY|  _ 17  JULY _ 18  JULY 
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TIME  (  h  ) 

Fig. ,9  Tempoial  variation  of  transmission  losses  and  predicted 
tidal  curve  at  site  C. 
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Spectral  composition  of  trans¬ 
mission  loss  at  site  C  for  two 
selected  frequencies  (SD  240  m  / 
PD  50  m). 
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Fig.  11 

Comparison  of  spectra  of  transioissv 
loss  at  25  Hz  and  sound  speed • 


Frequency  (cy/h) 


there  Is  a  correspondence,  albeit  inexact,  between  high  tide  and  high 
transmission  loss,  ;onsistent  with  the  fact  that  the  front  reaches  site  C 
during  periods  of  high  tide.  Although  the  curves  shown  are  for  a  par¬ 
ticular  sou:  -e  'receiver  combination  (240  m  /  50  m)  th  observed  trends  are 
fairly  representative  of  all  the  data  obtained  at  this  site.  The  spectral 
composition  of  the  50  Hz  and  1600  Hz  signals  is  shown  in  Fig.  10,  again 
demons  tr'.  ting  the  more  marked  fluctuations  at  the  lower  frequency,  as  well 
aa  the  c.  .ispiouous  presence  of  diurnal  and  semi-diurnal  periodicities. 

Figure  11  provides  a  comparison  between  the  25  Hz  transmission  loss 
spectrum  at  a  depth  of  20  m  (source  depth  ■  18  m)  and  the  sound  speed 
spectrum  at  a  depth  of  18  m.  As  seen  the  correlation  between  the  two 
spectra  is  quite  good. 

It  is  instructive  to  compare  the  variability  for  the  two  positions  B 
and  C.  An  example  of  the  fluctuations  in  temperature  is  shown  in  Fig.  12. 
Several  features  are  Immediately  evident  from  the  figure.  First,  it  is 
clear  that  the  amplitudes  at  C  are  greater  than  thosa  at  B  ind,  in  addi¬ 
tion,  show  a  definite  dependence  on  depth.  Further,  although  both  data 
sets  reveal  tiocl  effects,  they  are  more  conspicuous  at  position  C,  par¬ 
ticularly  at  18  m  depth.  This  behaviour  no  doubt  ref’ects  the  positif  i  of 
the  stations  with  respect  to  the  Polar  front,  station  C  being  dir.ut’y 
affected  by  the  movements  of  the  fro-f,  while  station  B  is  constantly  in 
isothermal  waters,  as  discussed  earlier.  The  significant  differences  in 
results  between  the  two  locations  suggest  that  tidally  advected  changes  in 
water  masses,  as  at  C,  are  more  important  than  the  indirect  tidal  effects 
(changes  in  water  depth,  currents,  etc.)  that  are  evident  at  B.  The 
greater  energy  measured  at  18  =  (for  site  C)  seems  to  confirm  this  conclu¬ 
sion  since  the  transition  zone  between  the  Atlantic  and  Polar  waters 
occurs  at  this  depth. 
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transmission  loss  are  greater  at  site  C  th in  1ow  tre^cy 
difference  between  the  sites  vis-a-vis  thf  m  k  \  bUt  there  is  Uttle 
loss.  This  Is  consistent  with  th®  atrnna  J'igtler  frequency  transmission 
at  C.  Wth  th«  stronger  downward  refracting  conditions 


CONCLUSIONS 


thCahB  that  fIuctuati°ns  in  scoSli^  tr/ns^"6  i"  ,  the  North  Atla"«c 

w/th  the  diurnal,  semi-diurnal,  sod  ine  tial  ^  1088  correlate  well 

environmental  parameters.  However  the  nit  .r  a  C  VarlaWllty  of  the 
Phenomena  differ  for  the  two  areas.’  At  thT  *  dir  .  °f  the  causal 

effects  appear  to  be  the  principal  cause  X  hf  X  T**?  te8t  slte  voluBe 
N.  Atlantic  test  sites  bottom  effects  nraa  ^  fluctuations,  while  at  the 
are  described  below:  effects  predominate.  Both  types  of  effects 


SACLANTCEN  CP-37 


3-11 


ALI  et  al:  Acoustic  transmission  loss  fluctuations 


Volume  Effects 


•  Inertial  effects  dominate  In  the  surface  layers,  whereas  the  semi¬ 
diurnal  effects  are  of  greater  importance  at  greater  depths.  This 
suggests  a  meteorological  forcing  function. 

•  The  fluctuations  in  the  magnitude  of  acoustic  transmission  loss  are 
greatest  for  the  higher  frequencies  (1.6  kHz  and  above)  and  least 
for  an  optimum  frequency  range  from  approximately  100  to  400  Hz. 

Bottom  Effects 

•  At  the  N.  Atlantic  site  the  inertial  period  coincides  with  the 
semi-diurnal  period  and  cannot  be  resolved  by  spectral  analysis. 

•  The  transmission  loss  oscillations  correlate  with  tidal  movements 
of  the  Polar  front  and  are  more  marked  at  the  lower  frequencies. 
The  effects  on  acoustic  propagation  are  associated  more  closely 
with  the  tldally-advected  changes  in  water  masses  than  with  ».'• 
more  subtle,  indirect  tidal  effects. 
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ABSTRACT 

Long-range,  low-frequency  directional  measurements  of  acoustic 
backscattering  from  the  seafloor  were  made  in  the  Tyrrhenian  Sea  using 
explosive  sound  sources.  The  signals  were  received  by  a  horizontal  towed 
array  of  hydrophones  and  were  processed  by  a  beamforming  procedure  to 
obtain  the  directional  distribution  of  the  scattered  field  as  a  function 
of  time.  These  data  are  used  to  form  images  of  scatterers  and  to  estimate 
the  backscattering  strength  of  specific  physiographic  features.  From  the 
data  obtained  the  scattering  strengths  were  estimated  to  range  from  -25  dB 
to  -35  dB  and  did  not  exhibit  strong  dependence  on  frequency. 


INTRODUCTION 

Low-frequency  acoustic  backscattering  from  the  seafloor  has  often 
been  studied  using  omnidirectional  measurement  techniques  [1-4]. 
Generally,  omnidirectional  sources  and  receivers  are  used.  The  calculation 
of  scattering  strength  assumes  uniform  scattering  in  a  ring-shaped  area  on 
the  seafloor  for  a  given  element  of  travel  time.  Such  measurements  can  be 
contaminated  by  non-bottom  returns  having  the  same  travel  time  as  the 
scattering  ring  [5]  or  by  non-uniform  roughness  within  the  scattering  ring 
[6].  Althoug..  careful  measurements  using  omnidirectional  geometries  will 
yield  proper  estimates  of  scattering  strength,  these  techniques  are  more 
limited  than  directional  methods  in  application.  Directional  methods  make 
it  possible  to  obtain  images  of  the  scattering  features  and  to  estimate 
the  scattering  strength  of  specific  physiographic  features  such  as 
seamounts  and  continental  slopes.  Images  of  seafloor  topography  from  long- 
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range  directional  measurements  of  acoustic  backscattering  have  been 
obtained  by  various  methods:  broad-side  beam,  side-scan  mode  with  towed 
array  [7],  full  beamforming  with  a  towed  array  [8,9],  and  full  beamforming 
with  two-dimensional  arctic  arrays  [10].  Backscattering  values  have  been 
reported  for  some  of  these  measurements  [9,10]. 

Previously  the  authors  have  reported  beamforming  methods  applied  to 
sensors  of  a  towed  array  receiving  Backscattered  sound  [11]  and  aspects  of 
resolution  of  seafloor  images  that  were  obtained  [12].  In  this  paper  we 
further  describe  the  technique  used  for  processing  broadband  data  received 
from  an  explosive  source  and  report  additional  measurements  of  bottom 
backscattering  strength  for  the  Tyrrhenian  Sea. 


EXPERIMENTAL  METHODS 

An  omnidirectional  explosive  source  (broadband)  is  deployed  near  a 
hydrophone  array.  The  sounds  reflected  and  scattered  back  to  elements  of 
the  array  are  processed  to  form  beams,  as  in  Fig.  1.  For  processing,  the 
hydrophone  signals  are  split  into  time  increments.  For  each  time  incre¬ 
ment,  a  beamforming  procedure  is  used  to  obtain  the  angular  distribution 
of  energy  over  a  given  frequency  band  of  the  source. 


Fig.  1 

Experimental  geometry  for  backscattering  measurement  (not  to  scale): 
returns  from  an  omnidirectional  sound  source  (not  shown)  are  received 
on  directional  receiving  beams.  The  beams,  which  actually  scan  a 
vertical  section  of  the  volume,  are  indicated  by  their  wedge-shaped 
areas  on  the  seafloor.  For  simplicity  only  three  beams  on  the  right- 
hand  side  of  the  array  are  shown  in  their  far  field  and  the  comple¬ 
mentary  (ambiguity)  beams  on  the  left-hand  side  of  the  array  are 
omitted. 
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and  is  based  on^he^assumption* deSCribed  itl  Refs-  [H]  and  [12] 
velocity  and  is  an  application  of  freaw^Ia  sound  signals  of  uniform 
the  Fast  Fourier  Transform  [131  Th  q  Cy”domaln  beamfarming  baaed  on 
line  array  results  in  the  estimation  r°cedure  applied  to  sensors  of  a 

“t  fZZTi z  - 

kx  «  (t>/c)  sin  a. 

w  =  ^  ^he  SthoTof  processing "broad^b  t],6.an8ular  acoustic  frequency. 
Fig.  2.  For  each  k*  value  at‘ tbe celt £  d“ta  ls  shown  «r-PMc2lly  in 

im~  *-  <•  *—STS  ST32  jyiS:  \s  - 


°  Data  points 
*  interpolated  points 


Fifif.  2 

“oooooq^ooo  Beam-power  interpolation  for 

oj-  §  §  S  o  o  °  °  °  wavenumber  k„:  beam-power  levels 

«  °  o  °o  lx*  l  l  l  l  ^e  interpolated  for  the  points 

°  °  °  °  O  o  o  indicated  by  stars  from  the  data 

s'  at  the  points  indicated  by 

/  circles.  The  interpolated 

s'  points  fall  along  the  azimuth 

_  line  (dotted)  of  the  central 

- — - r— - - ►  k*  frequency  of  the  band  m  . 

*n  C 

idn”e,„Th  Urj£  'hTP° «“lS.l*Sehthe,  «.!«. 

cor.esponds  to  the  central  frequency  of  the^nd  ^4  156818  azlnmth  that 
value  over  a  given  frequency  band  U  nhr  ^ a  avera*e  beam-power 
averaging  the  interpolated  power?,  mis  tlrtT  8t  “ch  **  value  by 
space  to  averaging  data  point?  along  Inline  corresponds  in  « 

sea t tering3  ring1  twundaries  fcf/tfi  Ih,  1?teraecti?tt  of  the  beam  and  the 
time  window,  as  in  Fig.  3!  WiLT?'  t?Vel  tlme  and  Ptocessing- 

t ering,  the  scattering  strength  f  for  TV*  ‘T™*  *ncoherent  scat- 

8  "  5,  for  a  given  frequency  hand  may  be 


cT 

2  cos  i3 

Fisr.  2  Plan  view  of  scatterinc  area  at 

of  scattering  area  for  small  grazing  angles^0^^  dimensions 
c  is  the  sound  speed  in  seawater; 

T  the  processing  time; 

<>  the  3-dB  beamwidth  (rad); 

e  the  grazing  angle  at  the  bottom ;  and 
R  the  range. 
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expressed  In  terms  of  the  sonar  equation  as 

S  «  RL  -  SL  +  TL1  +  TL2  -  10  log  A  ,  (2) 

2 

where  S  is  the  scattering  strength  in  dB/m  ;  RL  and  SL  are  respectively 
the  equivalent  reverberation  level  and  source  level  in  dB  re  1  uPa;  TL1 
and  TL2  are  the  one-way  transmission  losses  in  dB;  and  A  is  the  effective 
scattering  area  of  the  bottom  in  m2 .  If  the  process  involves  a  com¬ 
bination  of  scattering,  reflection,  and  diffraction,  then  the  estimated 
scattering  strength  may  have  errors  or  may  be  '.iappropriate.  However,  if 
the  process  is  entirely  that  of  specular  reflection  for  an  area  of  the 
seafloor  greater  than  a  Fresnel  zone,  a  reflectioi  oss  parameter  would  be 
more  appropriate. 


MEASUREMENTS 

Measurements  were  made  in  the  Tyrrhenian  Sea  at  the  locations  shown  in 
Fig.  4;  the  same  figure  also  indicates  the  major  scattering  features. 
Explosive  SUS  sources  (0.8  kg  TNT)  were  set  at  a  depth  of  245  m  near  a 
receiving  array  towed  at  a  depth  of  100  m.  The  hydrophone  signals  were 
split  into  increments  that  were  processed  by  the  beamforming  technique 
described  above  and  implemented  in  the  processing  system  presented  schema¬ 
tically  in  Fig.  5.  Repetition  of  the  beamforming  procedure  for  adjacent 


Fig.  4  Location  of  experiment:  estimated  positions  of  major  scattering 
features  are  shown  in  shading. 


Fig.  5  Acquisition  and  processing  system.. 


SACLANTCEN  CP-37 


4-4 


1 


! 


BERKSON  et  al:  Directional  measurements  of  backscattering 


time  segments  gave  a  display  of  the  beam-time  history  of  features  of  the 
seafloor  that  scattered  sound  back  to  the  array.  Examples  of  beam-time 
histories  are  shown  in  Figs.  6  and  7.  It  can  be  seen  that  the  direct  arrival  Is 
followed  by  arrivals  due  to  multiple  reflection  from  the  seafloor  at  high 
grazing  angles.  These  appear  strong  on  all  beams  because  of  sidelobe 
reception  of  the  overloaded  signals.  Later  arrivals  are  primarily  due  to 
backscattering  from  the  seafloor  at  low  grazing  angles.  The  records  also 


Fig.  6 

Beam-time  history  at  600  Hz:  SO-Hz  bandwidth  and  array  ientation  45° ; 
positive  azimuth  toward  ship;  explosive  source  occurs  at  t  -  6  s. 


:■  •  Fig.  7 

Beam-time  history  at  700  Hz:  50-Hz  bandwidth  and  array  orientation  270’; 
positive  azimuth  toward  ship;  explosive  source  occurs  at  t  =  6  s. 
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show  continuous  arrivals  (artifacts)  that  are  unrelated  to  returns  from  the 
explosive  source  and  are  comprised  of  towship  and  shipping  noise.  The  major 
features  are  caused  by  backscattering  from  major  physiographic  features, 
e.g.  island  slope  of  Sardina  and  the  Baconl  Seamounts. 

Scattering  strengths  were  calculated  using  Eq.  (1)  and  are  shown  in 
Fig.  8  as  a  function  of  frequency.  Also  shown  are  other  scattering  strength 
values  calculated  from  directional  measurements  at  low  frequency  [9,10]. 
The  frequency  dependence  of  the  scattering  strengths  of  these  data  is 
apparently  weak.  This  is  in  agreement  with  most  measurements  of  scattering 
from  the  seafloor  at  frequencies  of  2  kHz  to  100  kHz  [14].  More  oata  are 
required  before  the  relationship  at  low  frequencies  can  be  established. 
However,  at  this  point  these  limited  data  certainly  do  not  suggest  the  very 
strong  frequency  dependence  (about  1.6  power  dependence,  or  4.8  dB  per 
octave)  fcund  in  some  of  the  high-frequency  (13  to  290  kHz)  data  of 
McKinney  and  Anderson  [15]. 


Fig.  8 

Scattering  strength  versus  frequency:  data  for  coasts  of  South  Italy 
and  Alaska  taken  from  Refs.  9  and  10  respectively . 


CONCLUSION 

Directional  measurements  of  backscattering  were  made  with  an  oanl- 
directional  explosive  source  and  a  linear  towed  receiving  array.  The 
hydrophone  signals  were  processed  to  form  receiving  beams 
and  the  data  are  displayed  as  images  of  the  scattering  features  of 
the  seafloor.  The  scattering  strengths  obtained  exhibit  weak  frequency 
dependence,  but  more  data  are  required  before  the  frequency  dependence  of 
the  seafloor  features  can  be  established. 
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ABSTRACT 

Recently  It  has  become  apparent  that  shear  properties  of  ocean  sedi¬ 
ments  are  most  easily  determined  from  measured  propagation  characteristics 
of  seismic  Interface  waves.  We  use  a  newly  developed  numerical  model  of 
wave  propagation  in  horizontally  stratified  viscoelastic  media  to 
reconstruct  from  Scholte-wave  records  the  shear-speed  and  shear- 
attenuation  profiles  In  the  upper  sediment  layers  of  a  shallow-water  test 
area.  For  an  unconsolidated  sand-silt  bottom  we  find  the  shear  speed  to 
increase  with  depth  from  85  m/s  at  the  sea  floor  to  300  m/s  at  60  m  depth. 
On  the  other  hand,  the  shear  attenuation  at  3  Hz  decreases  with  depth  from 
0.45  dB/A  at  the  sea  floor  to  0.15  dB/A  at  60  m  depth,  with  an  approxima¬ 
tely  quadratic  frequency  dependence. 


INTRODUCTION 

It  is  well  established  that  shear  rigidity  of  the  ocean  bottom 
affects  propagation  of  waterborne  sound  through  the  conversion  of  acoustic 
energy  Into  shear  waves.  This  energy  conversion  is  of  particular  impor¬ 
tance  in  low-frequency  shallow-water  acoustics,  in  which  the  excitation  of 
shear  waves  in  the  bottom  often  becomes  the  dominant  loss  mechanism  for 
waterborne  sound  [1].  Under  these  circumstances  a  realistic  physical 
model  of  the  ocean  bottom  is  a  viscoelastic  solid  described  by 
corapressional  and  shear-wave  velocities,  by  the  attenuation  factors  asso¬ 
ciated  with  these  waves,  and  by  the  material  density. 

While  density  and  the  corapressional-wave  properties  (speed  and  atte¬ 
nuation)  in  f.ediments  can  be  determined  by  direct  methods,  the  shear-wave 
properties  are  difficult  to  measure.  This  is  because  of  the  usually  high 
attenuation  of  these  waves  and  because  it  is  difficult  to  generate  a  wave 
that  consists  of  predominantly  transverse  particle  motion.  The  shear 
speed  and  attenuation  can  be  indirectly  determined,  however,  through  the 
measured  propagation  characteristics  of  the  ocean-bottom  interface  wave. 
The  existence  of  this  wave  is  Intrinsically  related  to  the  shear  proper¬ 
ties  of  the  sediments. 

We  shall  review  the  characteristics  of  seismic  interface  waves  and 
give  some  simple  formulas  that  relate  their  speed  and  attenuation  to  those 
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of  shear  bulk  waves.  We  shall  also  provide  examples  from  the  literature 
of  Indirectly  determii.ed  shear  speeds  in  sediments.  Finally,  we  will 
demonstrate  how  the  use  of  a  sophisticated  numerical  model  of  seismic-wave 
propagation  allows  us  to  reconstruct  from  seismic  records  the  shear  speed 
jnd  shear  attenuation  profiles  in  the  upper  sediment  layers. 


SEISMIC  INTERFACE  WAVES 

This  particular  wave  type  -  sometimes  named  a  boundary,  surface,  or 
interface  wave  -  is  a  guided  wave  propagating  along  the  interface  between 
two  media  with  different  shear  speeds  (2,31-  The  wave  is  generally  given 
different  names  according  to  the  media  of  propagation  involved  { A  J . 
Hence,  if  propagating  on  a  free  surface  of  a  solid,  it  is  called  a 
Rayleigh  wave,  if  propagating  along  a  liquid/solid  interface,  it  is  called 
a  Scholte  wave  and  if  associated  with  a  solid/solid  boundary,  it  is  called 
a  Stoneley  wave.  Note  that  at  least  one  of  the  media  must  be  a  solid  for 
the  interface  wave  to  exist.  In  the  case  of  a  water/sediment  interface, 
the  pertinent  wave  type  is  a  Scholte  wave  with  the  following  charac¬ 
teristics: 

•  The  wave  propagates  along  the  sea  floor  with  exponentially 
decaying  amplitude  away  from  the  guiding  interface  (the  wave  is 
evanescent  in  both  media). 

■»  Its  particle  motion  is  elliptical  in  the  depth/range  plane. 

•  There  is  no  low-frequency  cutoff. 

•  Its  propagation  speed  and  attenuation  are  closely  related  to 
the  shear  properties  of  the  sediment. 

For  a  simplified  environment  consisting  of  two  homogeneous  media  in 
contact,  the  Scholte  wave  is  non-dispersive.  Its  propagation  speed  (Vscj,) 
and  attenuation  (Agch)  then  relate  in  an  extremely  simple  manner  to  the 
sediment  shear  properties  as:  , 

VSch  =  0.9  Vs  (1) 

and 

^  sch  =  ^  ^  s »  (  2  ) 

where  the  attenuation  is  in  dec  bels  per  unit  length.  Using  numerical 
simulation  these  formulas  were  determined  to  be  accurate  to  within  a  few 
percent  for  unconsolidated  sediments  (clay,  silt,  sand)  with  low  shear 
speeds.  The  formulas  indicate  that  the  Scholte  wave  properties  are  vir¬ 
tually  independent  of  densities  and  cospressional  wave  properties  in  the 
two  media  [5j.  Consequently  the  Scholte  Interface  wave  is  ideal  for 
determining  the  shear  properties  of  ocean  sediments. 

In  a  realistic  environment  with  sediment  layering,  the  Scholte  wave 
becomes  dispersive.  In  addition,  a  finite  number  of  ducted  shear  modes 
will  be  present.  Figure  1  shows  dispersion  curves  for  a  solid  layer  of 
low  shear  speed  on  a  high-speed  solid  substratum.  Th:  fundament.’’  mode 
(M0)  is  an  interface  wave  that  in  the  low-frequency,  limit,  kH  0,  propa¬ 
gates  as  a  Rayleigh  wave  on  the  substratum,  and  in  the  high  frequency 
limit,  kH  +  propagates  as  a  Rayleigh  wave  on  the  surface  of  the  upper 
layer.  The  higher-order  inodes  are  ducted  shear  modes  in  the  upper  layer. 
They  all  have  a  cutoff  frequency  at  which  they  propagate  with  the  shear 
speed  of  the  substratum.  In  the  high-frequency  limit  these  modes  propa¬ 
gate  with  the  shear  speed  of  the  upper  layer.  In  this  case  there  is  no 
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Fig.  1  Dispersion  carves  for  elastic  layer 

on  a  semi-infinite  elastic  substratum 
(  Tolstoy  and  Usdm  121  ). 

Stoneley  wave  due  to  the  stringent  conditions  necessary  for  the  existence 
of  this  type  of  wave.  A  Stonely  wave  exists  only  when  the  change  In  shear 
speed  across  a  solid/solJd  boundary  is  less  than  10%. 

In  the  following  section  we  show  how  the  complex  information  con¬ 
tained  in  experimentally  determined  dispersion  curves  similar  to  those  of 
Fig.  1  can  be  used  to  extract  shear  information  about  the  ocean  bottom. 
Such  an  analysis  generally  requires  the  use  of  a  numerical  model. 
However,  if  experimental  results  show  negligible  dispersion,  we  can  use 
simple  formulas  (Eqs.  1  and  2)  to  determine  shear  properties  for  the  upper 
sediment  layer.  No  dispersion  means  that  the  sediment  is  homogeneous  to  a 
depth  of  1  to  2  wavelengths,  the  approximate  penetration  depth  of  the 
Scholte  wave. 


SEDIMENT  SHEAR  PROPERTIES 

Over  the  past  20  years  various  investigators  [6-17]  have  used  Scholte 
waves  to  determine,  the  shear  properties  of  unconsolidated  sediments. 
Their  experimental  results  are  summarized  in  Table  1.  Note  that  experi¬ 
ments  were  done  in  water  depths  ranging  from  1  to  5260  m,  over  sandy  or 
silty  bottoms,  and  for  frequencies  between  3  and  35  Hz.  The  measured 
Scholte  wave  attenuations  (with  geometrical  spreading  removed)  vary  more 
than  three  orders  of  magnitude,  with  the  highest  attenuation  being  asso¬ 
ciated  with  the  highest  centre  frequency.  The  inferred  shear  speeds  are 
between  25  and  260  m/s  and  in  excellent  agreement  with  the  values  quoted 
by  Hamilton  [18], 

Assuming  the  shear  attenuation  to  be  linear  with  frequency  [18],  we 
can  conveniently  express  it  in  decibels  per  wavelength.  This  has  been 
done  in  the  last  column  of  Table  1,  using  the  simple  formulas  of  Eqs.  (1) 
and  (2)  to  relate  measured  attenuation  of  the  Scholte  wave  with  the  bulk 
shear-wave  attenuation.  We  obtain  values  that  still  vary  more  than  two 
orders  of  magnitude  (0.02  -  2.3  dB As),  whereas  Hamilton  [18],  based  on 
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Table  1  Interface  wave  experiments 


Investigators 

Year 

Water 

depth 

(»> 

Bottoa 

type 

Centre 

freq. 

(Hr) 

Measured 

art. 

(dB/k*) 

Inferred 
shear  speed 
(■/a) 

Inferred 
shear  att. 
(dBAg) 

sand 

20 

300 

100 

1.4 

Bucker,  Whitney,  Keir 

1964 

70 

sand 

75 

200 

195 

1.4 

Davie*7 

1965 

4410 

_ 

6 

_ 

50-190 

- 

Herron,  Doroan,  Drake* 

1968 

5 

blit 

5 

- 

40-115 

390 

slit 

_ 

_ 

100 

- 

Hamilton  et  al. 

1970 

983 

silt 

“ 

90 

“ 

.0 

Schlraer 

1980 

130 

o-nd 

4.5 

7 

120 

0.2 

McDaniel,  Beeb*11 

1980 

32 

sand 

10 

- 

200 

- 

,  12 

Essen  et  al. 

i68! 

l 

silt 

4 

- 

75-250 

- 

Tutblll  et  al,** 

1981 

7 

nud 

4.5 

- 

25-50 

- 

Whltraarsh,  LiWall*1* 

1982 

5260 

- 

4.5 

- 

25-170 

- 

1 « 

Holt,  Hovee,  Syrstad 
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sparse  In  situ  measurements  cor  silt  and  sand  sediments,  quotes  atte¬ 
nuation  values  around  0.5  -  2  .1B/AS.  Note,  however,  that  the  very  low 
attenuation  values  are  all  associated  with  low-frequency  experiments  (3  to 
5  Hz). 


There  are  two  possible  explanations  for  the  observed  large  variation 
in  shear  attenuation  with  frequency.  First,  that  the  shear  attenuation  in 
sediments  decreases  with  depth  as  reported  by  Hamilton  [181.  This  leads 
to  lower  attenuation  values  for  the  low-frequency  experiments,  since  the 
Scbolte  wave  here  penetrates  deeper  into  the  bottom.  The  second  explana¬ 
tion  is  that  the  attenuation  does  not  depend  on  the  first  power  of  fre- 
qi incy  as  assumed  above,  but  rather  follows  a  power  law  (f  )  where  the 
exponent  is  between  1  and  2.  Some  experimental  evidence  for  this  latter 
assumption  does  exist  [18].  In  particular,  Stoll  in  a  recent  paper  [19] 
has  strongly  advocated  the  hypothesis  of  a  marked  frequency  dependence  of 
attenuation  in  fluid-s'- curated  sediments,  based  both  on  recent  fi**ld  and 
laboratory  data,  and  on  theoretical  results  from  the  Biot  model.  However, 
more  experimental  work  is  needed  in  order  to  definitely  resolve  the  depth 
and  frequency  dependence  of  shear-wave  attenuation  in  marine  sediments. 


NUMERICAL  MODELLING  OF  SCHOLTE  WAVES 

Figure  2  displays  the  Scholte  wave  data  from  which  shear  speed  and 
shear  attenuation  profiles  can  be  determined  using  a  sophisticated  propa¬ 
gation  model.  Stacked  time  signals  for  both  the  vertical  and  horizontal 
(radial)  particle  velocities  are  shown  as  recorded  by  a  geophone  on  the 
ocean  bottom  [17].  The  source  was  an  explosi'  e  charge  detonated  near  the 
sea  floor  in  20  m  of  water.  The  charge  size  was  increased  with  range  as 
indicated  by  the  black  dots  in  Fig.  2.  We  clearly  see  rhe  dispersed  low- 
frequency  Scholte-wave  arrivals  with  group  velocities  between  78  and 
235  m/s. 
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Vertical  component 


b) 


Horizontal  componont 


Fig.  2  Stacked  time  signals  for  v  ■ tical  (a) 
and  horizontal  (b)  particle  velocities 
as  recorded  by  a  geophone  on  the  sea  floor. 


Before  performing  detailed  modelling  of  this  propagation  situation, 
it  is  convenient  to  carry  out  a  dispersion  analysis  of  the  experimental 
data.  By  applying  the  multiple  filter  technique  [17)  to  the  recorded  time 
series  at  range  1.7  km,  we  obtain  the  dispersion  diagrams  shown  in  Pig.  3, 
oo  j°"tours  indlcate  energy  levels  in  arbitrary  decibels  (max  level 
■  99  dB)«  It  is  evident  that  the  radial  component  of  the  particle  velo¬ 
city  gives  more  information  about  the  propagation  situation  than  does  the 
vertical  component.  Thus  Fig.  3b  shows  that  energy  is  arriving  in  three 
discrete  modes,  of  which  the  slowest  arrival  is  the  Scholte  mode  (M  )  with 
its  energy  centred  around  2  Hz.  The  first  shear  mode  (Mi)  is  strongly 
excited,  with  maximum  energy  around  2.8  Hz,  while  the  second  shear  mode 
(M2)  is  only  weakly  excited. 
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Frequency  (Hz) 


Fig.  3  Dispersion  curves  obtained  by  applying  a  multiple  filter 
technique  to  the  experimental  records  at  range  1.7  km. 
Energy  is  seen  to  be  propagating  in  three  discrete  modes: 
Mq  ,  M,  ,  and  M j  . 


The  modelling  Is  done  with  a  newly  developed  numerical  model  of 
seismic  wave  propagation  in  horizontally  stratified  media  [20,21].  The 
aim  is  to  construct  a  model  environment  leading  to  computed  dispersion 
characteristics  that  agree  with  the  experimental  results  of  Fig.  3.  The 
modelling  is  done  in  a  trial-and-error  fashion,  in  which  environmental 
parameters  are  changed  in  a  systematic  way  until  acceptable  agreement  is 
obtained  between  theory  and  experiment.  Since  we  are  interested  in  deter¬ 
mining  the  shear  properties  for  an  unconsolidated  sediment,  we  can  fix, 
a  priori,  the  compressional-wave  properties  and  densities,  which  are  known 
to  have  negligible  effect  on  the  propagation  characteristics  of  bottom 
interface  waves  [5].  The  compressional-wave  properties  (Cc,Pc)  and  den¬ 
sities  (p)  used  as  input  to  the  seismic  model  are  shown  in  Fig.  4.  Also 
shown  are  the  final  choices  of  shear-speed  (C8)  and  shear-attenuation  (S8) 
profiles  for  the  bottom. 
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Fig.  4  Model  environment  used  for  generating  synthetic  seismograms. 

The  computed  low-frequency  dispersion  curves  for  the  model  environ¬ 
ment  are  shown  in  Fig.  5.  In  Fig.  5a  we  have  superimposed  the  theoretical 
dispersion  curves  (dashed  lines)  on  the  experimental  data*  Clearly  there 
is  good  agreement  between  theory  and  experiment  for  the  modal  arrival 
structure,  particularly  for  the  two  lower-order  modes.  This,  in  turn, 
means  that  we  have  chosen  an  appropriate  shear-speed  profile.  Figure  5b 
shows  the  relative  energy  distribution  in  the  first  three  modes  as  deter¬ 
mined  from  the  numerical  model.  Again  there  is  good  agreement  with  the 
experimental  results  in  Fig.  5a,  indicating  that  the  choice  of  shear- 
attenuation  profile  is  also  appropriate.  It  should  be  pointed  out, 
however,  that  it  was  necessary  to  assume  a  quadratic  frequency  dependence 
of  shear  attenuation  in  order  to  get  agreement  on  energy  levels  over  the 
entire  frequency  band. 

The  inferred  shear  speeds  (85  to  300  m/s)  and  shear-speed  gradients 
(  <  6  m/s/m)  agree  well  with  values  given  by  Hamilton  [18]  for  sand-silt 
bottoms.  Reported  data  concerning  the  shear-attenuation  profiles  are  so 
sparse  that  no  comparison  with  those  in  the  literature  can  be  made. 

As  a  final  confirmation  of  the  validity  of  this  modelling  exercise, 
we  have  created  synthetic  seismograms  for  both  the  vertical  and  horizontal 
(radial)  particle  velocities  (Fig.  6).  We  notice  a  good  overall  agreement 
with  the  experimental  results  in  Fig.  2,  the  main  difference  being  fast, 
low-frequency  arrivals  present  in  the  synthetic  seismograms  but  filtered 
out  in  the  experimental  data. 


CONCLUSIONS 

Direct  measurements  of  shear  properties  of  ocean  sediments  are  dif¬ 
ficult  to  perform.  However,  both  shear-speed  and  shear-attenuation  profi¬ 
les  can  be  inferred  from  measured  dispersion  characteristics  of  bottom 
interface  waves  through  the  use  of  a  sophisticated  propagation  model.  The 
modelling  is  considerably  simplified  because  the  propagation  charac¬ 
teristics  of  Scholte  waves  depend  entirely  on  the  sediment  shear  proper¬ 
ties  and  have  only  negligible  dependence  on  the  compressional-wave 
properties  and  the  density.  For  a  shallow  water  test  area  we  find  sedi¬ 
ment  shear  speeds  (85  to  300  m/s)  and  shear  attenuations  (0.15  to 
0.45  dB/A  at  3  Hz)  that  are  in  good  agreement  with  values  reported  in  the 
literature.  Moreover,  in  accordance  with  recent  sediment  studies,  we  find 
the  shear  attenuation  to  have  an  approximately  quadratic  frequency 
dependence. 
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3)  Vertical  component 


b) 


Horizontal  component 


Fig.  6  Synthetic  seismograms  for  vertical  (a) 
and  horizontal  (b)  particle  velocities 
at  the  sea  floor. 
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EVALUATION  OF  LOW-FREQUENCY  BOTTOM  BACKSCATTERING  STRENGTH 
VS  GRAZING  ANGLE  BY  MEANS  OF  MULTIPLE  BEAMFORMING 


D.  Marandino  and  T.G.  Goldsberry 
SACLANT  ASW  Research  Centre 
Viale  San  Bartolomeo  400 
1-19026  La  Spezia,  Italy 


ABSTRACT 

Bottom  backscattering  strength  was  measured  in  deep  water  as  a  func¬ 
tion  of  grazing  angle  In  the  Balearic  Abyssal  Plain  of  the  Mediterranean. 
The  measurements  were  made  with  a  towed,  narrowband  low-frequency  omni¬ 
directional  source  and  a  towed,  horizontal  linear  array.  The  method  takes 
advantage  of  the  multiple  beamforming  capability  of  the  receiving  array 
and  processor  to  accurately  discriminate  returns  at  given  grazing  angles 
from  interfering  returns.  Approximations  for  the  evaluation  of  the  back- 
scattering  area  are  described  and  discussed.  The  results  are  presented  and 
compared  with  previously  reported  measurements.  At  the  higher  grazing 
angles,  specular  reflection  from  normal  incidence  dominates  the  returns. 
However  at  medium  and  small  grazing  angles  data  have  a  large  plateau 
region,  with  small  variations.  A  Lambert's  rule  relationship,  with  a 
coefficient  of  -33  dB,  can  approximate  results  for  grazing  angles  as  low 
as  about  20°. 

INTRODUCTION 

Low-frequency  bottom  backscattering  strength  values  have  been  pre¬ 
viously  measured  and  reported  by  several  authors  [1  to  5].  These  results 
refer  to  experiments  with  broad-band  sources  (explosives)  and  omnidirec¬ 
tional  receivers.  However  the  data  presented  herein  were  collected  using 
narrowband  sources  and  directional  receivers.  By  taking  advantage,  on 
reception,  of  the  directional  properties  of  a  towed  linear  array,  it  was 
possible  to  accurately  discriminate  returns  at  various  grazing  angles. 
Indeed,  it  is  well  known  that  beams  formed  by  a  linear  array  have  conical 
symmetry  around  the  array's  axis.  This  property  was  exploited  to  discri¬ 
minate  returns  in  the  vertical  plane  and  to  estimate  arrivals  from 
progressively  longer  times  and  thus  smaller  grazing  angles. 

Measurements  were  made  in  a  deep  water  area  in  the  Balearic  Abyssal 
Plain.  While  the  results  are  pertinent  only  to  a  similar  environment  the 
evaluation  criteria  can  be  generally  applied  to  horizontal  towed  arrays. 

The  next  section  provides  a  description  of  the  experiment  and  is 
followed  by  a  discussion  on  the  measurement  criteria  that  emphasizes  the 
evaluation  of  the  backscattering  area.  Finally,  the  results  are  presented 
and  discussed. 
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1  DESCRIPTION  OF  THE  EXPERIMENT 

The  experiment  was  conducted  in  an  area  of  the  Balearic  Abyssal  Plain 
which  is  typical  of  Mediterranean-type  abyssal  plains.  The  water  depth  at 
the  site  is  2800  o.  Cores,  6  m  long  that  were  previously  collected  in  the 
same  area,  show  layers  of  sand  interspersed  with  thin  layers  of  clayey  and 
silty  deposits.  An  isothermal  water  column  was  present,  giving  rise  to 
typical  Mediterranean  winter  propagation  conditions  characterized  by  a 
totally  upward  refracting  profile.  Wind  speed  was  less  than  7  lot  and  sea 
state  was  1  to  2,  with  a  moderate  swell. 

The  experimental  set-up  for  the  collection  of  monostatic  rever¬ 
beration  used  the  SACLANTCEN  R/V  Marla  Paolina  to  tow  a  sound  source  and  a 
linear  array  receiver  at  constant  speed,  course  and  depth.  During  the 
run,  the  source  generated  pulses  at  regular  intervals,  and  the  bottom 
echoes  were  received  by  the  towed  linear  array  find  subsequently  processed 
on  board.  Both  the  source  and  array  depth  were  100  m. 

Figure  1  shows  a  simplified  block  diagram  of  the  data  acquisition  and 
processing  system.  The  array  has  32  hydrophones  with  half  wavelength 
spacing.  Hamming  shading  is  applied  to  the  hydrophone  data.  The  filtered 
and  digitized  returns  are  beamformed  via  a  time-domain,  programmable 
multiple  beamformer,  and  subsequently  match-filtered  in  the  array  pro¬ 
cessor.  Figure  2  depicts  the  beamformer  file  that  was  used  during  the 
experiment  and  indicates  both  a  plan  and  side  view  of  the  beam-pointing 
directions,  as  determined  by  conical  symmetry  and  assuming  a  constant 
sound  speed  profile. 

The  transmitted  waveform  was  a  linear  frequency  modulated  (LFM)  pulse 
with  a  one  second  duration  and  10  Hz  swept  bandwidth  at  a  centre  frequency 
of  365  Hz.  The  received  signals  were  processed  through  a  quadrature 
replica  correlator  with  Hanning  shading;  the  final  equivalent  slant  range 
resolution  is  approximately  120  m. 


2  BOTTOM  SCATTERING  STRENGTH  EVALUATION  METHOD 
2.1  The  Data  Base 

Returns  from  a  number  of  pings  were  ensemble-averaged  to  determine  a 
mean  reverberation  envelope.  Reverberation  at  initial  ranges  is  dominated 
by  bottom  backscattering  that  affects  returns  in  two  distinct  ways: 

-  The  first  is  referred  to  as  the  "fathometer  effect",  and  is  caused 
by  the  energy  that  is  radiated  by  the  omnidirectional  source  in  the 
vertical  direction  and  that  propagates  vertically  in  the  water 
column,  repetitively  bouncing  off  the  jttom  and  the  surface.  This 
is  evidenced  in  the  processed  beam  output  time  series  as  a  sequence 
of  peaks  of  the  received  power,  which  are  regularly  spaced  at  a 
distance  equal  to  twice  the  water  depth. 

-  The  second  effect  is  the  progressive  decay  of  the  received  power 
which  is  caused  by  backscattering  froa  the  bottom  at  progressively 
longer  ranges,  and  thus  smaller  grazing  angles. 

Both  effects  can  be  seen  in  Fig.  3  where  the  ping-to-ping  (ensemble) 
average  of  the  returns  received  at  the  broadside  beam  over  a  limited  set 
of  ranges  is  shown  on  a  logarithmic  scale  in  an  A-scan  format.  Also  shown 
is  the  measured  ambient  noise  level  in  the  processed  band  (noise  floor). 
The  lower  half  shows  the  standard  deviation  of  the  averaging  procedure  of 
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Fig.  1  Block  diagram  of  SACLANT's  data  acquisition  and 
processing  system  for  reverberation  measurements. 
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Average  of  reverberation  returns  et  the  38 •  from  broaside  beam, 
ambient  noise  floor,  and  standard  deviation. 

(water  depth:  2800  mi  centre  frequency:  365  Hz) . 
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the  figure  which  has  a  mean  value  of  5.5  dB,  as  expected  from  the 
averaging  of  Rayleigh  distributed  log  normalized  data.  The  fathometer 
peaks  and  the  regular  decay  of  bottom  reverberation  vs  range  are  evident. 
However,  broadside  beam  data  cannot  be  used  to  accurately  determine  bottom 
backscattering  strength  because  of  the  high  sidelobes  associated  with  the 
fathometer  returns.  To  overcome  that,  the  data  from  the  off-broadside 
beams  have  been  used.  In  this  case,  the  fathometer  effect  is  attenuated 
because  vertically  propagating  energy  is  picked  up  through  the  angular 
sidelobes  of  the  various  beam  patterns.  Also,  the  beam  power  time  series 
shows  a  local  maximum  corresponding  to  the  onset  of  the  bottom  returns  at 
the  grazing  angle  related  to  the  main  beam  pointing  direction.  This 
correspondence  confirms  the  hypothesis  of  the  predominance  of  bottom 
reverberation. 

This  is  illustrated  in  Fig.  h  which  show  the  time  series  associated 
with  the  beam  pointing  at  38°  from  broadside,  which  corresponds  to  a  52° 
grazing  angle.  An  arrow  marks  the  main  beam  bottom  return.  The  advan¬ 
tages  of  evaluating  backscattering  strength  in  this  manner  are  that:  a)  a 
(local)  maximum  of  the  signal-to-noise  ratio  i3  used;  l)  the  influence  of 
the  fathometer  effect  can  be  considered  negligible;  c)  the  propagation 
loss  can  easily  be  estimated  by  assuming  spherical  spreading. 

Bottom  backscattering  strength  vs  grazing  angles  can  then  be  deter¬ 
mined  by  direct  application  of  the  sonar  equation,  using  the  power  output 
time  series  of  the  various  beams. 


2.2  Evaluation  of  the  Backscattering  Area 

In  order  to  determine  scattering  strength  by  means  of  the  sonar 
equation,  an  estimate  of  the  backscattering  area  (Aj,)  must  be  obtained. 
The  evaluation  of  Aj,  is  determined  by  three  factors:  propagation  mecha¬ 
nism,  the  characteristics  of  the  waveform  and  the  beam  shape  at 
reception. 

Figure  5  provides  a  plan  and  side  view  of  the  geometry  of  the 
experiment  which  is  applicable  to  the  evaluation  of  A^  and  defines  the 
pertinent  axes  and  parameters.  The  y-axls  is  orientated  along  the  arra> 
axis  direction.  In  the  figure, 

A^  is  backscattering  area  at  medium  to  small  "razing  angles 

A^rt  is  backscattering  area  at  high  grazing  angles 

dw  i3  water  depth 

h  is  water  depth  less  source  depth 

y  is  the  grazing  angle  (angle  'etween  bottom  and  propagation 

path) 

Yeq  is  the  "equivalent”  grazing  angle,  at  small  grazing  angles 
ar  is  the  annular  width  related  to  the  signal  range  resolution 

dj,  is  the  annular  width  related  to  the  mainbeam  Intersection  and 

dcr  is  the  "cross-range"  dimension  of  Aj,. 
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Fig.  5  Geometry  C or  the  backscattering  area. 


The  energy  radiated  by  the  source  propagates  spherically  and  inter¬ 
sects  the  (flat)  bottom  along  circles,  whereas  the  various  beams,  which 
define  cones  around  the  array  axis,  intersect  the  bottom  along  hyperbolic 
curves  B  and  B'  in  Fig.  5.  However,  such  representation  is  only  appli¬ 
cable  to  high  and  medium  grazing  angles  (  >  20°)  where  ray  bending  caused 
by  the  true  sound  speed  profile  are  negligible.  In  this  case  the  grazing 
angle  is  also  equal  to  the  complementary  angle  of  the  ray  at  the 
source/receiver.  For  smaller  grazing  angles,  say  below  20*,  it  has  been 
assumed  that  the  geometry  remains  approximately  valid,  provided  that  the 
actual  grazing  angle  is  replaced  by  an  ’’equivalent”  angle  obtained  l if  the 
straight-line  approximation  to  the  actual  Intercept  point  (IP  in  the 
figure)  as  determined  by  the  actual  sound  speed  profile.  For  this  experi¬ 
ment,  the  ray  at  13.86*  at  the  source  grazed  the  bottom  at  0*  with  an 
equivalent  grazing  angle  of  6.93®  and  slant  range  of  22370  m. 

For  an  omnidirectional  receiver,  power  comes  from  a  circular  annulus 
with  a  width  (dr)  expressed  by: 

dr  »  dg  /  cos(ir)  (1) 
where  d8  “  slant-range  signal  resolution  (here,  120  a). 

Since  the  receiver  is  directional,  only  that  part  of  the  annulus 
which  is  Included  in  the  malnbeam  contributes  to  backscattering,  under  the 
assumption  that  contributions  from  the  bean’s  sldelobes  are  negligible. 
The  intersection  of  the  malnbeam  vilit  the  bottom  (hyperboll'  curves  B  and 
B’  in  Fig.  5)  Is  expressed  by: 


2 


2  2 
x  +h 


2 


i 


2 


(2) 
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where  Yi  and  Y2  are  the  grazing  angles  corresponding  to  the  3  dB  bean- 
width. 

This  area  has  a  mlrJnnim  radial  width,  db,  along  the  trace  of  the 
array  axis  on  the  bottom,  expressed  by 

db  =  h  *  bw(90°)  /  (sin  (y))3  (3) 

where  bw{90<>)  =  array  beaowldth  at  broadside,  with  the  assumption  that  the 
beamwidth  for  any  beam  pointing  direction  is  given  by 

byOi)  =  byOO*)  /  sin(y). 

The  backscattering  area  can  then  be  estimated  through  the  rec¬ 
tangular  approximation: 

Ab  -  Dr  *  Dcr  (4) 

where  Dr  tbe  radial  dimension,  and 
^cr  is  the  cross-range  dimension. 

Two  cases  may  thus  be  distinguished: 

•  dr  >  db.  This  is  valid  at  high  grazing  angles  (AbH  in  Fig.  5). 
In  this  case  the  width  of  the  transmission  annulus  is  larger  than 
the  width  of  the  hyperbolic  annulus  generated  by  the  malnbeam 
intersection;  therefore,  in  Eq.  4  it  is: 

Dr  «  db  (5) 

•  and,  for  the  cross-range  dimension,  simple  geometrical  con¬ 
siderations  on  returns  included  between  R  and  R+d„,  where  R  is  the 
slant  range  in  the  malnbeam  give: 


«  dr  <  db.  This  is  valid  at  medium  to  small  grazing  angles.  The 
backscattering  area  (A^  in  Fig.  5)  is  that  part  of  the  circular 
annulus  that  is  delimited  by  the  intersection  with  the  beam  hyper¬ 
bola.  For  a  given  range,  the  area  can  then  be  expressed  by  the 
rectangular  approximation  (4),  where  the  radial  dimension  is 

0r  -  dr  ,  (7) 

and  the  cross-range  dimension  Dcr  is  the  arc  of  circle  generated  by 
the  intersection  of  the  lowermost  hyperbola  and  the  circle  related 
to  that  range.  Then: 


where  Y,  “  higher  grazing  angle,  which  corresponds  to  the  3  dB 
°  beamwidth  boundary 

Yr  “  grazing  angle  corresponding  to  the  given  range. 
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At  small  grazing  angles,  nowever,  ray  bending  Increases  the  distance 
at  which  the  rays  intercept  the  bottom  and  the  beam  intersections  with  the 
seabed  are  not  expressed  as  hyperbolas.  However  the  evaluation  of  is 
still  based  on  formulas  (7)  and  (8),  in  which  equivalent  grazing  angles 
that  correspond  to  the  linear  approximation  to  the  actual  rays  are  used. 

The  relationship  between  the  actual  and  the  equivalent  grazing  angle 
ccn  easily  be  established  since  the  experiment  was  carried  out  in  isother¬ 
mal  water  and  thus  with  a  constant  velocity  gradient.  In  this  case  the 
rays  are  arcs  of  circles  whose  centers  lay  at  a  fixed  distance 

Ro  *  CREF  I  ?  W 

where 

CrEF  is  c*le  reference  sound  velocity  (1500  m/s),  and 
C  is  the  velocity  gradient  (0.017  s-1) 

It  is  then: 

cos  Y  -  (1  +  — )  cos  YD4v  (10) 

Ro  “Y 

and 


Yeq  '  J  <*RAY  +  (11) 

where 

Y  is  the  angle  of  the  grazing  ray 

YRAY  *s  the  aRRle  t^'e  (•onostatic)  ray  at  the  source/array  and 

Yeq  is  the  equivalent  grazing  angle. 

The  approximations  used  in  the  evaluation  of  A;,  at  small  grazing 
angles  can  be  considered  satisfactory  because  dr  is  much  smaller  than  the 
other  dimensions. 

The  behaviour  of  the  bottom  backscattering  area  for  beams  off  broad¬ 
side,  as  a  function  of  the  true  grazing  angle,  is  Indicated  in  Fig.  6 
along  with  the  two-way  propagation  loss  at  the  pertinent  ranges.  A],  is 
expressed  in  dB  relative  to  one  square  aeter  and  is  evaluated  with  the 
parameter  values  applicable  to  the  experiment. 


Fig.  6 

Bottom  backscattering  area 
and  two-way  propagation  loss 
as  a  function  of  the  true 
grazing  angle. 
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rig.  7 

Backscattering  strength 
coefficient  s.,  versus 

D 

grazing  angle,  au  a 
centre  frequency  of  365  Hz. 
The  dotted  lines  show  the 
limits  of  reported  values 
for  low  frequency  bottom 
backscattering  strength . 


3  RESULTS  AND  CONCLUSIONS 


By  applying  the  sonar  equation  (6]  and  using  the  appropriate  curves 
for  the  backscattering  area  and  the  propagation  loss,  a  point  estimation 
of  the  backscattering  s' rength  coefficient  versus  grazing  angles  is 
obtained  (see  solid  line  In  Fig.  7). 

Two  cases  must  be  distinguished: 

•  At  high  grazing  angles,  because  of  the  signal  slant-range  resolu¬ 
tion  returns  from  73°  to  90"  cannot  be  separated  and  are  merged 
into  the  fathometer  peak  and  sidelobes.  Nevertheless,  the  fatho¬ 
meter  peak  itself  can  be  regarded  as  high  grazing  angle 
backscatter.  Therefore,  an  estimation  of  Sj,  at  near  incidence  can 
be  obtained  by  measuring  the  power  at  the  peak  of  the  fathometer 
return.  However,  since  the  first-bounce  fathometer  return  satura¬ 
tes  the  hydrophones,  it  is  impossible  to  determine  Sj,  reliably. 
Instead  the  second-bounce  fathometer  return  was  used  with  the 
addition  of  the  estimated  bottom  lose  incurred  at  the  first  bounce 
(6  dB).  The  applicable  backscattering  area  here  is  that  of  the 
broadside  beam  at  90°.  The  result  is  given  in  Fig.  7  by  the  value 
at  the  90"  grazing  angle. 

•  The  evaluation  of  Sb  for  medium  to  small  grazing  angles  is 
outlined  in  Sect.  2.2.  Both  rear  and  forward  beam  pointing  direc¬ 
tions  were  used;  however  forward  beams  have  a  higher  noise  level 
due  to  the  towship  contributions  and  >re  not  suitable  for  estima¬ 
tes  at  small  grazing  angles.  The  'imate  is  quite  accurate  for 
angles  above  30®  because  the  bacKscattering  area  has  small 
variations  and,  more  Importantly,  because  the  measurements  are 
made  before  the  onset  of  the  second  fathometer  return  and  are 
thereby  free  of  multiple  bounce  mattering.  At  smaller  grazing 
angles,  the  estimate  of  Sb  becomes  more  difficul '  because  the  out¬ 
put  ratio  of  signal-to-ambient  noise  is  lc  -er,  11  variations  of 
the  array  tilt  cause  large  variations  of  t  '  .v.ercept  area  and 
contributions  from  multiple-bounce  scatt-. .mg  are  more  likely. 
However  measurements  can  still  be  made  and  an  approximate,  upper 
bound  for  Sb  can  be  obtained  for  those  smaller  grazing  angles. 

The  overall  behaviour  of  the  bottom  backscattering  strength  indicated 
in  Fig.  7  shows  that  Sb  has  high  values  in  that  region  where  backscat- 
tering  is  dominated  by  specular  reflection.  The  range  of  angles  of  high 
backscattering  is  determined  by  the  system  geometry  and  the  type  of 
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waveform  such  that  the  coherent  (specular)  component  dominates  incoherent 
(diffuse)  scattering. 

For  decreasing  angles,  S(j  quickly  decays  by  about  30  dB.  Then,  for 
medium  and  small  grazing  angles,  it  has  a  plateau  region  where  the 
variations  are  small.  Values  of  -33  to  -37  dB  are  measured  for  medium 
angles  (60°  to  30°)  and  of  -AO  to  -42  dB  for  smaller  angles  (around  10°). 

A  comparison  with  Lambert's  rule  for  diffuse  backscattering  shows 
that  a  value  of  -33  dB  for  the  Lambert  coefficient  is  reasonable  for 
medium  grazing  angles.  The  limits  of  reported  data  for  low-frequency  bot¬ 
tom  backscattering  strength  range  from  -20  to  -40  dB,  and  are  shown  in 
Fig.  7  by  the  dotted  lines  (see,  for  instance,  Ref.  6,  page  246).  These 
include  results  at  various  frequencies  and  for  different  types  of  bottom. 
As  can  be  seen,  the  range  of  previously  reported  data  spans  the  present 
measurements,  although  similar  types  of  seafloor  are  usually  assumed  to 
have  slightly  higher  Sj,  values. 

In  conclusion,  a  method  has  been  presented  to  evaluate  low  frequency 
bottom  backscattering  strength  which  takes  advantage  of  the  multiple  beam¬ 
forming  capability  of  a  towed  linear  array.  The  intrinsic  system  direc¬ 
tivity  confines  backscattering  to  small,  homogeneous  seafloor  patches  and 
suppresses  Interference  from  other  directions. 
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ABSTRACT 

Variability  of  surficlal  sediment  geoacoustic  properties  was  deter¬ 
mined  from  cores  collected  at  eight  shallow-water  continental  shelf 
regions  in  the  U.S.,  Italy  and  Australia.  Highly  porous  muds  found  in  low 
energy  environments  exhibited  the  lowest  range  of  values  in  physical  and 
acoustic  properties;  mixtures  of  sand  and  shell  found  in  higher  energy 
environments  exhibited  the  highest  range  of  values.  Compressional  wave 
attenuation  consistently  exhibited  the  highest  variability  followed  by 
mean  grain  size,  porosity  and  compressional  wave  velocity.  Vertical 
variability  was  generally  greater  than  horizontal  variability  for  all  pro¬ 
perties  measured.  Sediment  geoacoustic  properties  of  most  coastal  marine 
sediments  are  controlled  by  the  Interaction  of  biological  and  hydrodynamic 
processes.  Biological  processes  tend  to  dominate  in  finer  sediments, 
whereas  hydrodynamic  processes  control  sediment  geoacoustic  properties  in 
sandy  substrates.  Understanding  these  processes  in  various  environments 
not  only  explains  the  spatial  distribution  of  sediment  geoacoustic  proper¬ 
ties  but  leads  to  Improvement  of  predictive  geoacoustic  models. 


INTRODUCTION 

Knowledge  of  the  spatial  variability  of  surficlal  sediment  geoacoustic 
properties  is  important  to  the  prediction  of  acoustic  scattering  from  the 
sediment-water  interface  and  to  the  prediction  of  propagation  of  acoustic 
energy  through  the  sediment1.  It  Is  within  this  surficlal  zone  (upper  50 
cm  of  sediment)  that  the  most  active  and  rapid  diagenic  changes  in  sedi¬ 
ment  properties  are  found.  Gradients  and  variability  of  geoacoustic  pro¬ 
perties  result  from  dewatering  caused  by  overburden  pressure,  as  well  as 
active  chemical,  sedimentological,  biological  and  hydrodynamic  processes 
which  mix  and  alter  sediment  properties2’  . 

For  high  frequency  (  >10  kHz)  acoustic  applications,  geoac.,ustic  pro¬ 
perties  of  the  upper  tens  of  centimeters  must  be  known,  whereas  for  low 
frequency  applications  surficlal  geoacoustic  properties  provide  the  ini¬ 
tial  conditions  used  for  prediction  of  depth  gradients  of  sediment  physi¬ 
cal  properties1.  Accurate  values,  Including  variability,  of  surficlal 
sediment  geoacoustic  properties  are  therefore  required  for  geoacoustic. 
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models  covering  the  wide  range  of  frequencies  of  interest  to  those 
studying  underwater  acoustics,  marine  sedimentology,  geophysics  and  marine 
geotechnique. 

In  this  paper  I  compare  data  on  the  spatial  variability  of  sediment 
geoacoustic  properties  among  eight  shallow-water  continental  shelf  regions 
in  the  U.S,  Italy  and  Australia.  I  also  discuss  the  relative  importance  of 
biological,  hydrodynamic  and  sedimentological  processes  in  determining  spa¬ 
tial  distribution  of  these  geoacoustic  properties.  Results  are  then  com¬ 
pared  to  a  similar  study  of  geoacoustic  properties  of  three  sedimentary 
provinces  in  the  Venezuela  Basin  (3500  to  50^0  m  water  depth). 


MATERIALS  AND  METHODS 

Replicate  sediment  samples  were  collected  from  eight  sites  along  the 
continental  shelves  of  the  United  States,  Italy  and  Australia  (Table  1). 
Sediment  types  ranged  from  clayey-silts  found  in  Long  Island  Sound  and  the 
Arafura  Sea  (off  northern  stralia)  to  coarse  sands  off  San  Diego, 
California.  All  sediments  vt  e  collected  with  6.1  cm  (inside  diameter) 
cylindrical  cores  either  in  situ  by  scuba  divers  or  from  relatively 
undisturbed  0.25  m2  USNEL  box  core  samples  collected  from  deeper  sites 
(Washington  coast  and  Arafura  Sea).  Collection,  measurement  and  handling 
procedures  were  designed  to  minimize  sampling  disturbance  and  maintain  an 
intact  sediment  water  interface  with  overlying  sea  water.  Sediment 
compressional  wave  velocity  and  attenuation  were  measured  after  sediments 
equilibrated  with  laboratory  temperature,  usually  one-half  to  two  days 
after  collection.  Temperature  and  salinity  of  the  overlying  water  were 
measured  with  a  YSI  model  43TD  temperature  probe  and  an  AD  Goldberg 
temperature-compensated  salinity  refractometer.  Compressional  wave  velo¬ 
city  and  attenuation  were  measured  at  1  cm  intervals  on  sediments  in  129 
subcores  using  a  pulse  technique. 


Table  1.  Collection  sites,  including  mean  values  of  porosity,  mean 
grain  size,  compressional  wave  velocity  (V -ratio) 
and  attenuation  (dB/m  at  400  kHz)  (* unpublished  data). 


Site 

Date 

water 

Depth 

(•) 

Substrate 

type 

Porosity 

(X) 

Crsln  Site 
<♦> 

Vp-ratio 

Attenuation 

(dl/m) 

Keferencca 

Long  Island  Sound 

xm  so 

m 

FOAH 

10 

clayey-silt 

73.2 

0.99 

NWC 

16 

clayey-silt 

77.2 

" 

San  Diego, CA 

IV  -  V  82 

i 

line  sand 

18 

very  fine  sand 

- 

1.10 

188 

■ 

coarse  sand 

18 

coarse  sand 

” 

1.0 

1.13 

116 

HontauV  Point, NY 

V  82 

35 

fine  sand 

36.6 

2.1 

1..4 

88 

Qulnault,  Washington 

IV  83 

49 

fine  sand 

41.2 

2.9 

1.11 

160 

Charleston  S.C. 

VI  83 

18 

medium  sand 

37.9 

1.6 

1.12 

292 

La  Spezla,  Italy 

XI  83 

MT 

17 

gravel 

- 

-1.6 

1.20 

- 

ST 

sllty-clay 

69.8 

0.98 

104 

rv 

stlty-clay 

67.6 

0.98 

64 

AV 

■tfl 

very  fine  sand 

43.3 

■  n 

1.10 

136 

Arafura  Sea 

V  84 

clayey-sand 

69.2 

0.99 

336 

a 

Panama  'Uty,  Florida 

IX  84 

V-  ■•! 

10 

fine,  medium  sand 

39.3 

1.8 

1.14 

216 

experimental 

m 

fine  sand 

39.0 

2.6 

1.13 

228 
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Time  delay  measurements  were  made  through  sediments  and  a  distilled 
water  reference  with  a  Hewlett-Packard  1743A  dual  time  Interval 
Oscilloscope.  Signals  were  generated  by  driving  a  Underwater  Systems, 
Inc.  (Model  USI-103)  transducer-receiver  head  with  a  400  kHz,  20  volts  p- 
p  sine  wave  triggered  for  25ps  duration  every  10  ms  with  a  Tektronic  FG 
SOI  Pulse  Generator  and  FG  504  Function  Generator.  Differences  In  time 
delay  between  distilled  water  and  sediment  samples  were  used  to  calculate 
sediment  compressions!  wave  velocity  (Vp).  Compresslonal  wave  velocity 
was  expressed  as  the  dimensionless  ratio  of  measured  sediment  -elocity 
divided  by  the  velocity  in  the  overlying  water  that  was  calculated  for  the 
same  temperature,  salinity  and  depth11.  This  ratio  is  independent  of 
sediment  temperature,  salinity  and  depth  and  therefore  ideal  for  com¬ 
parison  to  other  geoacoustic  properties.  Attenuation  measurements  were 
calculated  as  20  log  of  the  ratio  of  received  voltage  through  distilled 
water  to  received  voltage  through  sediment12.  Values  of  attenuation  were 
extrapolated  to  a  1  m  pathlength  and  expressed  as  dB/m. 

After  completion  of  the  acoustic  measurements,  core  samples  were 
extruded  and  sectioned  at  2  cm  intervals  for  determination  of  sediment 
porosity  and  grain  size  distribution.  Porosities  were  determined  by 
weight  loss  of  sediment  dried  at  105*0  for  24  hours.  Values  were  not 
corrected  for  pore  water  salinity.  Salt-free  porosity  may  be  obtained  by 
multiplying  values  by  1.012.  Sediment  grain  size  distribution  was  deter¬ 
mined  on  disagregated  samples  by  dry  sieving  for  sand-sized  particles  and 
with  a  Micromeritics  Sedigraph  and/or  pipette  for  silt  and  clay-sized  par¬ 
ticles.  Grain  size  statistics  were  determined  using  the  graphic  formula 
of  Folk  and  Ward13. 


RESULTS 

Sediment  geoacoustic  data  were  measured  op  129  cores  that  were 
collected  for  this  study  at  eight  different  sites.  I  have  restricted  the 
data  presentation  to  a  summary  of  sediment  geoacoustic  properties  from 
each  site  together  with  selected  graphical  representations  of  the  spatial 
variability  of  sediment  geoacoustic  properties.  Detailed  data  presen¬ 
tation  can  be  found  elsewhere3-1®. 

Sediment  samples  were  collected  from  two  sites  in  Long  Island  Sound. 
Both  sites  were  characterized  by  high  porosity,  fine  grained,  low  velocity 
sediments.  At  the  deeper  (18  m)  NWC  site  sediments,  were  characterized  by 
a  uniform  distribution  of  geoacoustic  properties  (Table  2)  on  scales  of 
centimeters  to  meters  with  no  apparent  depth  gradients  in  the  upper  20  cm. 
Richardson  et  al.  attributed  this  low  variability  to  sediment  mixing  by 
macrofaunal  animals  that  feed  on  bottom  deposits  (bioturbation).  At  the 
shallower  FOAM  site  (10  m)  sediments  exhibited  a  much  higher  variability 
in  sediment  geoacoustic  properties  (Table  2).  This  higher  variability  was 
a  result  of  storm  induced  eroslonal  and  depositions!  events  creating  lami¬ 
nations  witiun  the  upper  35  cm  of  sediment.  These  laminations,  each  with 
different  values  of  geoacoustic  properties,  were  preserved  because  sedi¬ 
ment  mixing  by  macrofaunal  animals  at  the  FOAM  site  rarel/  extended  below 
the  upper  few  centimeters  of  sediment3. 

The  shallow  water  samples  collected  off  San  Diego  came  from  two 
distinct  sediment  types:  a  fine  sand  that  had  lower  compresslonal  wave 
velocity  and  higher  attenuation  than  the  coarser  sand  over  which  it  was 
migrating.  Of  particular  Interest  were  the  significant  positive  gradients 
in  compresslonal  wave  velocity  and  attenuation  with  depth  (Fig.  1)  for 
both  substrate  types  without  a  concordant  change  in  mean  grain  size. 
These  gradients  probably  resulted  from  Increased  packing  and  compaction  of 
sands  with  depth  in  the  sediment. 
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Fig.  1  Regressions  of  compressional  wave  velocity  ratio  (Vp-ratio)  and 
attenuation  (dB/m  at  400  kHz)  with  depth  for  fine  (+)  and  coarse 
(a)  sand  substrates  found  off  San  Diego,  California. 

The  fine  sand  sediments  collected  in  35  m  of  water  25  km  east  of 
Montauk  Point,  Long  Island,  New  York  exhibited  relatively  low  variability 
in  sediment  geoacoustic  properties,  especially  attenuation.  This  low 
variability  was  attributed  to  mixing  of  sediments  by  an  abundant  popula¬ 
tion  of  sand  dollars6.  The  positive  gradient  of  compressional  wave  velo¬ 
city  with  depth  with  little  change  in  mean  grain  size  was  similar  to 
gradients  found  in  sediments  from  San  Diego.  Porosi':y  generally  decreased 
with  depth  supporting  the  argument  chat,  for  sand,  positive  gradients  in 
compressional  wave  velocities  can  result  from  increased  compaction. 
Surficial  sediment  samples  collected  within  a  0.5  km?  area  around  the 
prime  collection  site  exhibited  considerable  variability  in  mean  grain 
size  (CV»108)  because  the  study  site  was  located  at  the  base  of  a  drowned 
barrier  spit.  At  this  site  coarse  grained  lag  deposits  formed  during  the 
last  glacial  advance  were  covered  intermittently  by  fine-grained  sediments 
which  were  in  equilibrium  with  current  hydrodynamic  conditions. 
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Fine  sand  surficial  sedinents  collected  in  50  ra  water  depth  on  the  con¬ 
tinental  shelf  off  Washington  State  exhibited  little  variation  in  values  of 
geoacoustic  properties  from  the  surface  to  12  cm  depth  (Fig.  2).  Below 
12  cm,  the  sediments  had  higher  percentages  of  silt-sized  particles  that 
resulted  in  higher  values  of  porosity,  mean  grain  size,  and  attenuation  as 
well  as  lower  values  of  compress  tonal  wave  velocity.  Variability  of 
geoacoustic  properties  was  also  highest  below  12  cm  depth.  A  high  abun¬ 
dance  of  molluscan  shells  below  12  cm  also  contributed  to  higher  variabi¬ 
lity  of  geoacoustic  properties,  especially  attenuation. 

2 

Scuba  divers  observed  two  surficial  substrate  types  in  a  one  km  study 
area  off  the  Charleston,  S.C.  coast.  Fine  to  medium  sand  (mean  grain  size 
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Fig .  2  Depth  distribution  of  values  of  compressional  wave  attenuation 
(dB/m  at  400  kHz)  and  porosity  (X)  for  sediments  collected  in 
the  Quinault  acoustic  tracking  range  -  Washington . 


1.8  to  2.2$)  with  porosity  values  of  37  to  42 X  were  found  in  troughs  of 
shallow  sand  waves.  The  coarser  (0.9  to  1.7$)  sediments  on  crests  con¬ 
tained  abundant  shell  material  and  generally  had  lower  values  of  porosity 
(30  to  40Z).  Below  6  to  10  cm,  shell  material  was  abundant  in  all  cores. 
Coefficients  of  variation  for  geoacoustic  properties  of  the  Charleston 
samples  were  higher  than  at  most  other  locations  (Table  2).  Most  of  the 
variability  at  the  Charleston  site  was  associated  with  samples  that  had 
greater  than  IX  gravel-sized  particles.  If  these  samples  are  excluded, 
coefficients  of  variation  are  reduced  to  0.77.Z  for  compressional  wave  velo¬ 
city;  23.2515  for  attenuation;  5.19J5  for  grain-size;  and  2.20Z  for 
porosity. 

Sediment  samples  were  collected  from  four  locations  in  the  vicinity  of 
La  Spezla,  Italy  as  part  of  an  international  program  to  compare  various 
acoustical,  geotechnical  and  geophysical  measurement  techniques9. 
Locations  were  chosen  to  include  a  wide  spectrum  of  sediment  types 
(Table  1).  The  Monasteroli  site  (MT)  was  located  in  a  high  energy  environ¬ 
ment  next  to  the  Ligurian  coast.  Sediments  were  composed  of  gravel-sized 
particles  (75Z  by  weight)  making  it  difficult  to  collect  adequate  samples. 
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The  coarser  sediments  (mean  grain-size  -1.57$)  probably  originated  from 
erosion  of  the  nearby  steep  coastline9 .  Frequent  storms  winnowed  finer 
material  from  the  gravel  sediments  leaving  only  0.173!  silt-and  clay-sized 
particles.  Compressional  wave  velocity  ratio  (mean  1.20)  and  its  coef¬ 
ficient  variation  (3.452)  were  higher  than  at  any  other  location. 
Porosity,  mean  grain  size  and  attenuation  were  not  accurately  measured  to 
permit  the  calculation  of  coefficients  of  variation.  Sediments  from  the  PV 
site  located  at  the  eastern  entrance  of  the  Portovenere  channel  were  pro¬ 
tected  from  physical  disturb-^ce,  e.g.  storms.  Sediments  were  composed  of 
silt  and  clay  mixed  with  up  to  52  gravel-sized  shell  particles. 
Compressional  wave  velocity  and  mean  grain  size  values  varied  little 
(Table  2).  However,  the  high  coefficient  of  variation  for  porosity 
reflected  10  to  152  higher  values  of  porosity  found  in  the  uppermost  4  cm 
of  sediment.  Both  biological  (bioturbation)  and  hydrodynamic  processes 
probably  maintained  this  higher  surficial  porosity.  The  Santa  Teresa  site 
(ST)  was  located  in  a  small  protected  harbor  on  the  eastern  side  of  the 
Gulf  of  La  Spczia.  Samples  consisted  of  very  soft  muds.  Compaction  of  the 
upper  6  cm  produced  a  negative  gradient  of  porosity  with  depth  without 
apparent  change  in  sediment  acoustic  properties.  The  lack  of  vertical  gra¬ 
dients  in  compressional  wave  velocity  resulted  in  a  low  coefficient  of 
variation.  The  Venere  Azzurra  site  (AV)  was  located  on  a  shallow  (6  m), 
flat,  hard-packed  sandy  bottom  off  Lerici  in  the  Gulf  of  La  Spezia.  These 
well  sorted  and  hard  packed  sediments  had  the  lowest  coefficient  of 
variation  for  values  of  attenuation  of  any  site  listed  in  Table  2. 

Twelve  sediment  cores  were  collected  from  a  one  km2  area  in  the  central 
Arafura  Sea.  The  clayey- sand  sediments  contained  an  abundant  quantity  of 
both  whole  and  broken  gravel-sized  molluscan  shells  (2  to  232).  The  lack  of 
orientation  of  shells  suggested  considerable  biological  mixing  had 
occurred.  Variability  of  mean  grain  size  was  a  function  of  the  patchy 
distribution  of  shells  whereas  variability  of  porosity  was  related  to  ver¬ 
tical  gradients  caused  by  bioturbation  and  by  strong  currents  which 
resuspend  sediments.  Very  high  values  of  attenuation  measured  at  400  kHz 
(mean  336  dB/m)  were  a  direct  result  of  scattering  of  the  high  frequency 
signal  from  shell  material.  Even  at  125  kHz  values  of  attenuation  were 
high  (mean  65  dB/m)  and  highly  variable  (45.72).  Compressional  wave  velo¬ 
city  ratio  appeared  to  be  much  less  affected  by  shell  material,  and  was 
controlled  primarily  by  the  clayey-sand  matrix  (0.989  at  125  kHz  vs  0.988 
at  400  kHz). 

As  part  of  a  pre-site  survey,  27  sediment  cores  were  collected  from  a 
variety  of  sandy  sites  (18  to  33  m  water  depth)  off  Panana  City,  Florida. 
The  upper  2  cm  of  sediments  ranged  from  coarse  to  fine  sand  (0.42$  to 
2.75$)  with  a  relatively  narrow  range  of  porosities  (36.3  to  43.92).  The 
coefficient  of  variation  for  attenuation  and  compressional  wave  velocity 
were  approximately  the  same  as  other  sandy  sites  (see  Table  2).  As  part  of 
the  acoustic  experiment  15  cores  were  collected  from  a  one  kn£  area  that 
had  uniform  surficial  sediment  properties10.  As  expected  the  coefficient 
of  variation  for  all  geoacoustic  properties  was  lower  than  that  calculated 
from  the  pre-site  survey  (Table  2,  Fig.  3). 


DISCUSSION 

The  data  in  Table  2  demonstrate  that  shallow-water  geoacoustic  proper¬ 
ties  of  sediment  can  be  quite  variable  on  scales  of  a  kilometer  or  less. 
Not  evident  from  Table  2  are  the  spatial  sources  of  this  variation  (either 
vertical  or  horizontal),  the  scales  (cm  to  km)  or  the  processes  that  create 
this  variation.  Understanding  these  relationships  can  increase  the  predic¬ 
tability  of  sediment  geoacoustic  properties  in  time  and  space. 
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a)  ATTENUATION  (dB/m) 


b)  ATTENUATION  (dB/m) 
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Fig.  3  Comparison  of  the  variability  and  depth  distribution  of  attenuation 
(dB/m  at  400  kHz)  for  a  one  km 7  experimental  site  (a)  and  a  larger 
presite  survey  area  (b)  off  Panama  City,  Florida. 


Several  general  observations  can  be  made  from  the  data  presented  in 
Table  2.  Congressional  wave  attenuation  consistently  had  the  highest 
values  of  variation,  and  velocity  the  lowest.  These  results  were  expected 
given  the  previously  reported  range  of  values  for  compresslonal  wave  velo¬ 
city  and  attenuation12*  l4*  15 .  The  coefficient  of  variation  for  mean 
grain  size  was  generally  higher  than  for  porosity.  The  lowest  coefficients 
of  variation  for  mean  grain  size  and  congressional  wave  velocity  were  found 
at  muddy  sites,  whereas  attenuation  and  porosity  exhibited  the  lowest  coef¬ 
ficients  of  variation  in  fine  hard  packed  aands.  The  highest  variability 
of  attenuation  was  found  at  sites  with  an  high  percentage  of  shell 
material. 

Depth  gradients  in  geoscoustic  properties  at  several  sites  (FOAM,  San 
Diego,  Montauk  Point,  Washington  and  Charleston)  accounted  for  a  con¬ 
siderable  percentage  of  variation  of  sediment  geoacoustic  properties.  At 
the  San  Diego  site  compresslonal  wave  velocity  and  attenuation  signifi¬ 
cantly  correlated  with  depth  (F-test,  >0.001)  for  both  fine  and  coarse  sand 
substrates  (see  Fig.  1).  The  coefficients  of  variation  for  these  proper¬ 
ties  after  correction  for  depth  gradients  were  11.59  to  11.79  for  atte¬ 
nuation  and  0.72  to  0.79  for  velocity.  Depth  gradients  in  geoacoustic 
properties  at  the  San  Diego  and  Montauk  Point  sites  were  the  result  of 
increas'd  packing  of  sand-sized  particles  with  depth5*6.  At  the  FOAM  site 
sediment  laminations  produced  by  storm  Induced  eroslonal  and  deposltlonal 
events  resulted  in  the  higher  values  of  mean  grain  size  and  porosity  and 
lowest  values  of  compresslonal  wave  velocity  between  5  to  15  cm  depth  in 
the  cores3*4.  Increased  percentages  of  silt  with  depth  at  the  Washington 
continental  shelf  site  resulted  in  higher  values  of  porosity  and  lower 
values  of  compresslonal  wave  velocity  below  12  cm  depth  in  the  cores7.  At 
both  the  Washington  and  Charleston  sites  the  higher  abundance  of  shell 
material  with  depth  resulted  in  higher  values  of  attenuation  with  depth7'6. 
Knowledge  of  the  depth  distribution  of  values  of  sediment  geoacoustic  pro¬ 
perties  can  greatly  increase  predictability  of  these  properties. 
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The  magnitude  of  horizontal  variation  in  geoacoustic  properties  varied 
considerably  among  sites.  At  the  FOAM  site  in  Long  Island  Sound  horizontal 
patchy  distribution  of  animals  resulted  in  considerable  variation  in 
geoacoustic  properties  on  scales  of  10  to  100  meters4.  At  Ran  Mego,  sedi¬ 
ments  were  classified  as  fine  and  coarse  sands  based  on  a  mosiac  created 
from  overlapping  tracks  of  side-scan  sonar  records5.  Knowledge  of  the 
distribution  of  these  two  substrate  types  greatly  reduced  coefficients  of 
variation  of  geoacoustic  properties  used  in  correlations  with  the  values  of 
acoustic  bottom  backscatter1® .  Side-scan  mosaics  from  Che  Charleston,  S.C. 
and  Panama  City  sites  also  Indicated  considerable  spatial  variability  in 
sediment  types.  The  improved  predictability  of  sediment  geoacoustic  pro¬ 
perties  for  the  Panama  City  site  is  clearly  demonstrated  in  Table  2.  At 
the  Charleston  site  both  horizontal  and  vertical  variations  in  sediment 
geoacoustic  properties  must  be  taken  into  account  to  improve  the  predic¬ 
tability  of  sediment  properties.  Large  scale  horizontal  variation,  such  as 
found  in  data  collected  from  the  Gulf  of  La  Spezia,  can  usually  be 
accounted  for  by  utilizing  detailed  sediment  distribution  naps  which  are 
available  for  many  continental  shelf  regions. 

Table  3.  Dominant  factors  affecting  geoacoustic  properties  of 
sediment  (xx  =  very  important;  x  =  important) . 


Experiment 

Hydrodynaeilc 

Processes 

Biological 

Processes 

Relict 

(Inherited) 

Features 

Long  Island  Sound 

X 

XX 

San  Diego  CA 

XX 

X 

Montauk  Point 

X 

X 

XX 

Quinault,  Wash 

XX 

Charleston  S.C. 

XX 

X 

X 

La  Spezia,  Italy 

XX 

XX 

Arafura  Sea 

X 

XX 

X 

Panana  City,  Fl. 

XX 

X 

Sediment  geoacoustic  properties  of  many  coastal  marine  sediments  are 
controlled  by  the  interaction  of  biological  and  hydrodynamic  processes3. 
Table  3  presents  the  relative  importance  of  these  processes  as  well  as  that 
of  relict  (inherited)  features  in  determining  values  and  variability  of 
sediment  geoacoustic  properties.  Bioturbation  (mixing  of  sediment  by 
deposit-feeding  animals)  tends  to  reduce  horizontal  variation  in  sediment 
physical  properties  especially  in  fine  grained  sediments  (NWC,  Arafura 

Sea).  The  vertical  gradients  created  by  bioturbation  are  usually 

restricted  to  the  upper  10  cm  and  are  predictable  on  a  seasonal  basis  given 
knowledge  of  biological  processes  present.  A  high  abundance  of  tube 

dwelling  species  tends  to  stabilize  the  sediment  surface,  thus  preserving 
the  spatial  variability  of  geoacoustic  properties  (FOAM).  In  fine  sediment 
where  the  dominance  of  hydrodynamic  and  biological  processes  alternate 
horizontal  and  vertical  variability  is  great  and  predictability  poor 
(FOAM). 

In  sandy  substrates,  surface  deposit  feeding  species,  such  as  sand 

dollars,  tend  to  reduce  horizontal  variability5.  Feeding  also  reduces 
packing  in  the  upper  few  centimeters  increasing  porosity  and  decreasing 
corapressional  wave  velocity  (Montauk  Point,  Panama  City).  At  Montauk  Point 
sand  dollars  were  so  abundant  that  they  acted  as  surface  point  scatterers 
of  acoustic  energy17.  The  presence  of  surface  and  burled  shell  material 
contributes  to  considerable  fine  scale  variation  of  geoacoustic  properties 
(Quinault,  Charleston,  Arafura  Sea  and  Panama  Cir  >ites).  The  high  values 
and  variability  of  corapressional  wave  attenuate .  i  associated  with  shell 
material  are  due  to  scattering  from  the  shells  as  opposed  to  the  intrinsic 
absorption  described  by  Hamilton12.  The  distribution  of  shells  in  the 
sediment  is  controlled  by  both  biological  and  hydrodynamic  processes. 
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The  random  distribution  and  lack  of  preferred  orientation  of  shells  at  the 
NWC  and  Arafura  Sea  sites  suggests  biological  mixing  processes  were  domi- 
.  ant,  whereas  hydrodynamic  mixing  processes  control  the  distribution  of 
shells  at  sandy  sites  (Charleston,  Qulnault  and  Panama  City). 

Large  scale  variability  Is  generally  controlled  by  hydrodynamic  pro¬ 
cesses  acting  on  given  rel<ct  features.  At  the  San  Diego  site  a  recently 
deposited  fine-grained  sand  was  found  to  be  migrating  over  a  coarser- 
grained  offshore  Pleistocene  deposit6.  At  the  Montauk  Point  site  a  light 
colored  fine  grained  sand  of  modern  origin  discontinously  covered  a  reddish 
granular  sediment  which  was  a  lag  deposit  formed  by  the  erosion  of  a 
a.  owned  Pleistocene  barrier  spit6.  At  both  locations  the  sand  of  modern 
origin  was  in  dynamic  equilibrium  with  recent  storm  events,  whereas  the 
coarser  Pleistocene  deposits  were  in  equilibrium  with  winter  storm  events. 

Hamilton  and  Bachman14  presented  data  on  geoacoustic  properties  from 
340  sediment  samples  collected  on  the  continental  shelf  and  slope.  They 
classified  sediments  into  nine  groups  based  on  grain  size.  The  coefficient 
of  variation  for  porosity  (range  5.0  to  13. 8Z;  mean  10. 3Z)  and 
compresslonal  wave  velocity  ratio  (0.7  to  3.2Z;  mean  2.3Z)  were  much 
greater  within  those  nine  sediment  classes  then  for  the  individual  sites 
included  in  this  study  (Table  2).  Direct  measurements  of  geoacoustic  pro¬ 
perties  are  therefore  preferred  to  predictions  based  on  sediment  type 
despite  considerable  fine-scale  variability  of  geoacoustic  properties 
demonstrated  in  this  paper. 

Briggs  et  al18,  using  the  same  techniques  employed  here,  characterized 
sediment  geoacoustic  property  variability  for  three  sedimentary  provinces 
in  the  Venezuela  Basin  (3450  to  5050  m  water  depth)  (see  Table  2).  Most  of 
the  variability  was  associated  with  vertical  gradients  of  geoacoustic  pro¬ 
perties  (cm)  as  opposed  to  horizontal  variability  (on  scales  of  cm  to 
km)18.  As  would  be  expected  in  the  deep  sea  coefficients  of  variation  for 
the  carbonate  and  hemi pelagic  sites  were  lower  than  at  most  shallow-water 
sites,  wheress  the  turbidite  site  had  very  'nigh  coefficients  of  variation. 
This  variability  resulted  from  the  presence  of  alternating  layers  of  pela¬ 
gic  and  terrigeneous  sediments*  This  study  indicates  the  importance  of 
considering  'he  variability  of  geoacoustic  properties  for  geoacoustic 
modelling  of  both  deep-sea  and  shallow  water  sediments. 
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EVALUATION  OF  EXPERIMENTAL  TECHNIQUES  FOR  DETERMINING 
THE  PLANE  WAVE  REFLECTION  COEFFICIENT  AT  THE  SEA  FLOOR 


Henrik  Schmidt  and  Finn  B.  Jensen 
SACLANT  ASW  Research  Centre 
1-19026  La  Spezla,  Italy 


ABSTRACT 

One  of  the  most  commonly  used  techniques  for  determination  of  the 
geoacoustical  properties  of  the  sea  bed  is  the  measurement  of  plane  wave 
reflection  coefficients  at  the  ocean  bottom.  An  incident  field  is  produced 
by  either  an  explosive  source  or  a  beam  generating  device  and  the 
reflected  field  is  then  detected  by  means  of  an  array  of  hydrophones.  The 
associated  angles  of  Incidence  have  traditionally  been  determined  by 
simple  geometrical  considerations.  However  such  simple  interpretations 
give  results  that  in  some  cases  depend  on  the  actual  experimental 
geometry.  For  example  critical  angles  can  appear  far  away  from  the  correct 
values  or  not  be  present  at  all.  Here  an  exact  numerical  model  Is  used  to 
examine  the  different  experimental  techniques.  The  observed  discrepancies 
are  explained,  both  for  the  point  and  beam  source  experiments.  In  addi¬ 
tion,  guidelines  are  given  for  interpretation  of  results  obtained  by  the 
different  experimental  techniques. 


INTRODUCTION 

The  importance  of  the  sea  bed  geoacoustic  properties  for  shallow 
water  sound  propagation  is  well  established,  and  reliable  transmission 
loss  predictions  obtained  by  means  of  numerical  propagation  models  there¬ 
fore  require  accurate  knowledge  of  the  sea  bed  properties.  The  traditional 
raytrace  propagation  models  require  a  plane  wave  reflection  coefficient  at 
the  bottom  in  order  to  account  for  the  bottom  loss.  More  recent  wave 
theory  models  like  those  based  on  normal  modes,  the  full  wavefield  fast 
field  programs  (FFP)  and  the  parabolic  equation  models  (PE),  require  a 
more  detailed  knowledge  of  the  wave  speeds,  attenuations  and  densities  in 
the  sea  bed.  These  parameters  could  in  principle  be  obtained  from  samples, 
but  due  to  the  fact  that  the  low-frequency  acoustic  waves  penetrate 
deeply,  very  deep  and  expensive  boreholes  would  be  required.  Further,  and 
often  more  important,  the  de-pressurisation  and  change  of  temperature, 
unavoidable  In  the  core  sampling  process,  tend  to  deteriorate  the  mechani¬ 
cal  and  chemical  bondings  in  the  sediment  material,  and  thus  heavily 
influence  the  properties,  shear  in  particular. 

The  geoacoustic  properties  therefore  primarily  have  to  be  determined 
from  in-situ  propagation  experiments.  One  of  the  most  common  experimental 
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techniques  Is  the  determination  of  the  plane  wave  reflection  coefficient 
at  the  sea  floor.  This  approach  has  two  advantages:  the  results  can  be 
used  directly  by  the  raytrace  models,  and  the  plane  wave  reflection 
coefficient  is  needed  as  an  input  parameter  in  most  inverse  schemes  for 
the  determination  of  bottom  properties  fl). 

Several  experimental  techniques  have  been  devised  with  the  objective 
of  determining  the  reflection  coefficients  directly.  They  are  all  based 
on  the  detection,  by  hydrophones,  of  the  bottom  reflected  part  of  the 
field  produced  by  a  sound  source  placed  in  the  water  column.  The  source 
has  been  either  an  omnidirectional  explosive  source  or  a  device  producing 
a  narrow  beam  of  sound  directed  towards  the  sea  floor  at  a  variable  angle 
of  incidence.  However,  several  authors  have  shown  that  the  results 
obtained  are  usually  not  directly  interpretable  as  plane  wave  reflection 
coefficients. 

Here  we  will  use  the  full  wavefield  SAFARI  model  [2,3],  to 
demonstrate  how  these  discrepancies  arise,  and  to  show  how  numerical 
models  are  not  restricted  to  interpretational  purposes,  a  traditional 
application  in  underwater  acoustics,  but  can  also  be  used  for  the  design 
and  planning  of  experimental  setups. 


EXPLOSIVE  SOURCE  EXPfiRIKENTS 

The  most  common  experimental  technique  for  determination  of  plane 
wave  reflection  coefficients  uses  an  explosive  omnidirectional  source  to 
generate  a  transient  field.  The  hydrophone  array  used  as  a  detector  may 
be  either  horizontal  or  vertical,  moving  or  fixed  (Fig.  1). 

By  assuming  that  the  source  and  receivers  are  so  far  apart  that  the 
bottom-interfering  eigenrays  can  be  considered  plane  waves  when  hitting 
the  bottom,  a  very  simple  interpretation  technique  has  been  used.  First 
the  nominal  specular  reflection  angle  at  the  bottom  is  determined  for  each 
receiver  by  means  of  simple  raytracing.  Then  the  received  signals  are 
split  into  a  direct  part  and  a  bottom  reflected  part  by  inspection.  After 
correction  for  different  travel  paths,  the  reflection  coefficient  is  found 
simply  by  dividing  the  frequency  spectrum  of  the  reflected  signal  by  that 
of  the  direct  signal. 


Fig.  1  Experimental  set  up  for  bottom  reflection-loss  measurements  by 
means  of  an  explosive  source  and  a  vertical  hydrophone  an-j. 
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The  obvious  advantage  of  the  explosive  source  experiments  is  the 
possibility  of  separating  the  different  arrivals  in  the  time  domain,  and 
hence  to  eliminate  unwanted  surface  multiples.  Further,  the  experiments 
are  cost-effective  because  they  do  not  require  any  specialized  equipment. 
As  pointed  cut  by  several  authors,  however,  the  simple  interpretation 
technique  outlined  above  is  only  rarely  applicable.  Stickler  [A]  and 
Santanlello  et  al  [5]  demonstrated  that  the  simple  interpretation  tech¬ 
nique  will  give  wrong  reflection  coefficients  when  the  bottom  is  upward 
refracting  or  has  deeper,  reflecting  interfaces,  because  of  the  inter- 
*e~ence  between  the  different  arrivals.  Non-physical  effects  like  negative 
•-  om  loss  and  source/receiver  r jsition-dependent  results  arise.  The 

.cal  angle  shift  pointed  out  by  Sticklet  [4]  is  due  to  the  same 
e  ct.  Even  in  the  case  of  an  isovelocity  bottom  with  a  sound  speed 
nigher  than  that  of  th'  ocean,  the  headwave  formation  will  give  rise  to 
the  same  phenomenon.  The  simple  Interpretation  technique  can  be  applied 
only  in  the  rather  unusual  case  of  a  purely  downward  refracting  or  homoge¬ 
neous  bottom  with  a  sound  speed  less  than  that  of  water. 

We  will  here  illustrate  the  limitations  of  the  simple  interpretation 

principle  by  simulating  an  explosive  source  experiment  in  a  very  simple 

ocean  environment  by  means  of  .he  SAFARI  full-wavefield  model.  The  sound 
speed  profile  is  shown  in  Fig.  2.  The  water  is  characterized  by  an  upward 
refracting  profile  close  to  the  bottom,  and  sure«ce  mu’tiples  of  no  pre¬ 
sent  interest  are  avoided  by  replacing  the  ocean  aoovi  -800  m  by  an  isove¬ 
locity  halfspace.  The  fluid  bottom  is  upward  reflecting  to  30  m  below  the 

seabed.  Below  this  depth  it  is  represented  by  an  isovelocity  halfspace. 
In  order  to  more  clearly  illustrate  the  phenomenon,  the  bottom  is  con¬ 
sidered  lossless.  Thus  the  reflection  loss  is  identically  zero  for 
grazing  angles  less  than  the  critical  14.6°  onto  the  water-bottom  inter¬ 
face. 


An  explosive  source  is  assumed  to  be  placed  400  m  above  the  bottom 
and  the  radiated  pressure  pulse  has  a  duration  of  10  m3  and  a  centre 


Fig.  2  Sound  speed  profile  fc;  simulation  of  experiment . 
Depth  0  m  refers  to  the  sea  bed. 
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Fig.  3  Transmission  loss  contours  at  100  Hz  for  a  point  source 
400  m  above  the  bottom.  Hydrophone  array  used  for  pulse 
calculations  is  indicated  by  o . 

frequency  of  100  Hz.  Figure  3  shows  the  transmission  loss  contours  in 
depth  and  range  at  the  centre  frequency  (b'.ack  indicating  highest 
intensity).  The  lloyd-mirror  pattern  due  to  the  interference  between  the 
direct  and  the  bottom  reflected  fields  is  evident  and  illustrates  the 
complexity  of  the  sound  field  even  in  this  very  simple  case. 

A  vertical  array  of  7  hydrophones  with  100  m  spacing  is  placed  3  Urn 
from  the  source  as  indicated  in  Fig.  3.  The  synthetic  hydrophone  signals 
for  this  array  are  shown  :.n  Fig.  4,  with  each  trace  being  identified  by  the 
nominal  sp^ .  ilar  angle  of  the  reflected  signal  at  the  ocean  bottom,  as 
separated  f-om  the  reflected  parts,  which,  however,  do  not  clearly 


Direct  arrival 


Reflected  arrival 


ojj:  -»uu  m 

Range:  3  km 


Time  (seconds) 

Fig.  4  Synthetic  hydrophone  signals  for  vertical  array  at  3  km 
range.  Each  trace  is  identified  by  its  corresponding 
nominal  specular  reflection  angle. 
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Frequency  (Hz) 


Fig.  5  Reflection  coefficients  determined  by  deconvolution 
of  synthetic  signals:  a)  9.3°,  b)  11.4° 


indicate  the  critical  angle  of  14.6°.  Further,  the  trace  corresponding  to 
9.3°  indicates  an  apparent  negative  bottom  loss.  These  properties  are 
also  indicated  by  the  reflection  coefficients  obtained  by  the  simple 
deconvolution  principle  for  the  grazing  angles  9.3°  (Fig.  5a)  and  11.4° 
(Fig.  5b)  respectively.  Both  angles  are  less  than  critical;  therefore  the 
reflection  loss  should  be  zero  in  both  cases.  The  "ringing"  at  low  and 
high  frequencies  is  due  to  the  very  low  energy  content  of  the  source  pulse 
at  these  frequencies,  but  even  in  the  central  frequency  interval  errors  of 
several  dB  are  obtained. 

Figure  6  outlines  the  different  travel  paths  yielding  errors  in  the 
reflection  coefficients  obtained  by  the  simple  technique.  The  possibility 
of  multiple  paths,  due  to  headwaves  (1),  upward  refracting  profiles  (2)  or 
reflecting  interfaces  (3),  yields  results  which  are  dependent  on  whether 
the  different  arrivals  are  interfering  constructively  or  destructively, 
which  is  again  dependent  on  the  source-receiver  positions.  As  is  also 
clear  from  Fig.  6,  the  multiple  arrivals  do  not  correspond  to  the  same 
angle  as  the  nominal  specular  reflection  angle  at  the  water-bottom 


SACLANTCEN  CP-37 


8-5 


SCHMIDT  &  JENSEN:  Plane  wave  reflection  coefficients 


Fig.  6  Additional  arrivals  obstructing  simple  deconvolution. 
1)  Headwaves,  2)  Upward  refracted  arrivals,  3)  Deep 
reflection . 


interface.  The  field  detected  by  each  hydrophone  is  therefore  not  a  single 
plane  wave  component,  as  is  assumed  when  the  simple  interpretation  tech¬ 
nique  is  used,  but  a  complicated  interference  between  several  different 
components . 

This  phenomenon  directly  leads  to  the  conclusion  that  in  order  to 
obtain  the  plane  vsve  reflection  coefficient  from  explosive  source  experi¬ 
ments,  some  kind  of  beamforming,  or  plane  wave  de-.omposition,  has  to  be 
applied  to  the  received  signals.  In  principle  this  could  be  done  by  using 
towed  arrays.  However,  at  the  low  grazing  angles,  often  of  main  interest, 
this  technique  requires  very  long  arrays  because  the  bottom  refracted 
signals  travel  a  long  distance  before  re-entering  the  water  column, and  the 
whole  reflected  field  has  to  be  covered  by  the  array.  Therefore,  a 
synthetic  aperture  technique,  as  proposed  by  Frisk  et  al  [6],  is  often 
more  convenient  than  the  use  of  a  towed  array. 


BEAM  EXPERIMENTS 

In  contrast  to  the  omnidirectional  explosive  source  described  above, 
a  source  that  generates  an  approximate  plane  wavefield  would  provide  the 
opportunity  to  directly  measure  the  plane  wave  reflection  coefficients. 
This  is  the  philosophy  behind  the  use  of  beamed  sources,  as  shown  schema¬ 
tically  in  Fig.  7.  A  beam  is  directed  towards  the  water-sediment  inter¬ 
face  at  a  nominal  angle  of  incidence,  and  the  reflected  beam  is  measured 
in  the  specular  direction  by  means  of  a  hydrophone. 

It  is  well  known,  however,  that  very  wide  beams  are  required  in  order 
to  simulate  plane  wave  behaviour.  This  is  due  to  the  fact  that  a  beam  of 
finite  width  has  a  finite  ®ngular  spectrum,  whereas  for  a  plane  wave  the 
angular  spectrum  is  infinitely  narrow  [3].  For  practical  reasons  this 
technique  has  therefore  not  been  a  realistic  alternative  to  the  explosive 
source  experiments  until  the  development  of  the  parametric  transducer  [7]. 
In  principle,  this  transducer  generates  a  virtual  endfire-array  ''hat  is 
the  source  of  a  highly  directional  beam  which  can  be  used  for  measurement 
of  plane  wave  reflection  coefficients  [8],  In  practice,  however,  the  vir¬ 
tual  array  will  have  a  finite  length,  and  thus  results  in  the  beam  having 
a  finite  angular  spectrum. 
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We  will  here  use  the  SAFARI  model  to  demonstrate  how  a  velocity  pro¬ 
file  in  the  ocean  close  to  the  bottom  influences  the  reflection  of  a 
realistic,  narrow  beam  at  low  grazing  angles.  The  environment  shown  in 
Fig.  2  is  also  used  for  this  study.  A  vertical  source  array  of  61  ele¬ 
ments  at  half- wavelength  spacing,  placed  400  m  above  the  bottom  is  used 
to  generate  a  narrow  beam,  6  wavelengths  wide  at  the  sea  bed,  measured 
between  the  3-dB  down  points.  The  array  is  phased  to  yield  a  nominal 
grazing  angle  of  2.5°,  58  and  10°  at  the  sea  bed. 

The  resulting  fields  are  shown  in  Fig.  8  (black  indicating  highest 
intensity).  The  contour  interval  is  2  dB,  but  the  actual  dB  values  are 
arbitrary.  Obviously  the  beams  are  not  specularly  reflected.  At  the  two 
smallest  grazing  angles,  Figs.  8a  and  8b,  a  clear  beam  splitting  occurs, 
and  the  reflected  beam  with  highest  amplitude  has  a  much  smaller  grazing 
angle  than  the  incident  beam.  A  comparison  with  Fig.  3  shows  that  the 
directions  of  the  split  beams  correspond  to  those  of  the  lowermost 
lobes  of  the  Lloyd-mirror  pattern  in  the  point  source  field.  This  indica¬ 
tes  that  the  angular  spectrum  of  the  beam  is  apparently  so  wide  that  the 
interference  between  the  direct,  but  upward  refracted,  purely  waterborne 
arrivals  and  those  reflected  off  the  bottom  is  impcrtant.  This  assump¬ 
tion  is  supported  by  the  calculated  angular  spectrum  of  the  5°  beam  at  the 
bottom,  shown  in  Fig.  9.  The  beam  is  seen  to  contain  significant  energy 
at  grazing  angles  in  the  interval  of  1°  through  12°.  Figure  8c  indicates 
that  the  splitting  effect  due  to  the  interference  decreases  for  higher 
grazing  angles,  as  expected,  but  although  the  energy  is  here  concentrated 
in  a  sing:*  beam,  a  significant  widening  of  the  beam  cross-section  has 
occured.  Even  at  this  relatively  high  grazing  angle,  a  single  hydrophone 
in  the  specular  direction  would  not  yield  the  right  value  of  the  reflec¬ 
tion  coefficient;  instead  it  would  indicate  too  high  a  reflection  loss. 

Even  in  the  ideal  case  of  an  isovelocity  water  column  over  a  homoge¬ 
nous  bottom,  it  was  shown  by  Muir  et  al  [8],  that  narrow  parametric  beams 
could  penetrate  into  the  bottom  at  grazing  angles  less  than  critical. 
This  problem  has  been  treated  theoretically  by  several  authors.  Tjotta 
and  Tjotta  [9]  presented  a  wide  bean  approximation  showing  the  effect 
qualitatively.  Schmidt  and  Jensen  [3)  used  the  SAFARI  model  to  simulate 
the  experiments  in  which  extremely  narrow  beams,  not  covered  by  the 
approximate  theory,  were  used.  The  reported'  beam  cross-section  at  the 
bottom  interface  was  simulated  by  means  of  l.  focusing  linear  array. 
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Fig.  9  Angular  spectrum  of  beam  at  the  ocean  bottom 
for  nominal  grazing  angle  of  5°. 

The  results  demonstrated  that  subcrltlcal  penetration  is  a  simple  con¬ 
sequence  of  the  basic  physical  principle  that  a  narrow  beam  has  a  wide 
angular  spectrum.  Thus  the  narrow  beam  may  contain  energy  that  propagates 
at  grazing  angles  greater  than  critical  although  the  nominal  angle  of 
incidence  is  subcritical.  The  results  further  showed  that,  for  the  same 
reason,  the  reflected  beam  was  not  specular,  but  shifted  towards  a  smaller 
grazing  angle.  In  conclusion,  the  parametric  beams  in  practice  are  narrow 
and  not  highly  directional  as  often  stated,  and  hence  they  are  a  poor 
approximation  to  a  plane  wave. 

To  summarize,  the  use  of  beam  sources  does  not  yield  the  possibility 
of  measuring  the  plane  wave  reflection  coefficients  directly  by  means  of  a 
single  hydrophone.  As  was  the  case  for  the  explosive  sources,  plane  wave 
decomposition  has  to  be  applied  to  the  reflected  field.  However,  the 
advantage  of  using  beam  sources  is  that  smaller  arrays  or  synthetic  aper¬ 
tures  can  be  used  because  only  a  limited  part  of  the  angular  spectrum  is 
activated  by  the  beam. 


CONCLUSION 

A  full  vavefield  numerical  model  has  been  used  to  demonstrate  that 
neither  <  :plosive  source  techniques  nor  narrow  beam  techniques  yield  the 
plane  wave  reflection  coefficients  directly.  It  has  been  demonstrated, 
that  the  narrow  beams  are  not  highly  directional,  and  thus  do  not  behave 
like  plane  waves.  Therefore,  both  experimental  techniques  require  the 
application  of  either  plane  wave  decomposition  or  beamforming  to  the 
reflected  field  in  order  to  correctly  determine  the  plane  wave  reflection 
coefficients. 
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INTERFACE  WAVE  STUDIES  ON  THE  LIGURIAN  SHELF  USING  AN  OBS  ARRAY: 

EXPERIMENTAL  RESULTS  AND  PROPAGATION  MODELS 

Michael  Snoek,  Guido  Guidl,  Enzo  Michelozzi 
SACLANT  ASW  Research  Centre 
Viale  San  Bartolomeo  400 
1-19026  La  Spezla,  Italy 


ABSTRACT 

Seismic  waves  travelling  In  tie  water/sediment  or  subbottora 
sediment/sediment  interface  have  been  the  subject  of  considerable  interest 
In  recent  years.  Experiments  have  confirmed  the  existence  of  Interface 
waves  in  different  geological  environments.  Some  progress  has  been  made 
In  understanding  the  propagation  and  attenuation  characteristics  of  these 
waves.  However,  the  generating  mechanisms  are  poorly  understood.  In  par¬ 
ticular,  what  are  the  roles  of  geometry  and  bottom  characteristics  in 
coupling  waterborne  sound  into  the  sea  bottom?  What  is  the  acoustic-to- 
seismic  energy  conversion  process?  Using  explosives  as  sound  sources,  an 
experiment  was  performed  in  which  the  propagation  along  different  inter¬ 
face  paths  was  measured  by  an  array  of  three  tri-axial  ocean  bottom 
seismometers.  To  obtain  realistic  input  data  for  the  propagation  models, 
continuous  seismic  profiles  were  run  and  cores  were  taken  along  the  propa¬ 
gation  path.  Results  of  these  experiments  and  a  comparison  of  measured 
and  modelled  data  are  presented. 


INTRODUCTION 

Results  from  SACLANTCEN  experiments  and  other  published  data  show 
that  the  propagation  and  attenuation  characteristics  for  seismic  interface 
waves  are  still  poorly  understood  and  are  ofttiraes  ambiguous  (All  A 
Schmalfeldt  [1]  and  Jensen  &  Schmidt  12]). 

Rauch  [3,4]  has  outlined  the  theoretical  background  for  the  mathema¬ 
tical  treatment  of  interface  waves.  Jensen  &  Schmidt  [2]  have  also  pre¬ 
sented  a  modelling  approach  to  reconstruct  shear  speed  and  shear 
attenuation  from  interface  waves  for  the  upper  relevant  sediment  layers. 

Unclear,  however,  is  how  and  under  which  conditions  acoustic  energy 
converts  into  seismic  energy.  It  has  been  frequently  observed  that  in 
changing  the  bearing  of  a  profile  or  the  frequency  spectrum  of  the  source 
propagation,  interface  wave  characteristics  can  cnange  dramatically.  The 
changes  can  range  from  signal  enhancement  to  disappearance.  Before  a 
reasonable  prediction  of  attenuation  and  propagation  characteristics  can 
be  made,  the  driving  mechanism  for  interface  wave  conversion  and  propaga¬ 
tion  have  to  be  carefully  studied. 
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To  describe  this  conversion,  one  must  determine  the  sensitivity  of 
the  process  to  differences  in  the  physical  properties  and  geometry  of  the 
surface  layers,  to  geological  structure,  to  the  frequency  content  of  the 
seismic  source  and  to  the  bearing  and  position  of  the  source  with  respect 
to  the  receiving  sensors.  The  interrelations  of  these  parameters  affect 
the  complexity  of  the  conversion  process.  The  interpretation  of  the  data 
has  to  be  supported  by  a  modelling  phase  based  on  input  data  which  should 
be  as  complete  as  possible.  This  requires  a  multi-parameter  survey 
involving  methods  and  systems  schematically  shown  in  Fig.  1. 


SEISMICS  SEAFLOOR  PROFIUNG 


Fig.  1 

Direct  and  indirect 
methods  of  investi¬ 
gations  for  a 
detailed  ocean 
bottom/subbottom 
survey 


INTERFACE  WAVE  STUDIES 

We  have  conducted  a  set  of  experiments  on  the  Ligurian  shelf  off  La 
Spezia,  Italy,  to  study  the  interface  waves  along  different  propagation 
paths  (Fig.  2)  in  a  region  where  we  could  expect  silty  clay  and  clayey 
silt  layers  over  a  more  sandy  3ubbottom. 

The  instrumentation  and  the  setup  of  SACLANTCEN  interface  wave  experi¬ 
ments  have  been  described  in  detail  by  lauch  &  Schmalfeldt  [5)  and  All  & 
Schmalfeldt  [1].  Wp  operated  with  two  ships  and  deployed  three  tri-axial 
ocean  bottom  seismometers  in  a  triangular  configuration,  with  distances  of 
537  m,  352  m  and  315  m  between  them.  All  distances  were  relative  to  a 
marker  buoy.  Positions  were  controlled  by  Loran  C,  satellite  navigation, 
radar,  and  terrestrial  navigation. 

Charges  were  fired  electrically  on  the  sea  bottom,  their  sizes  ranging 
from  180  g  (standard  size)  to  1080  g.  The  data  were  available  a  short 
time  after  each  shot,  thus  allowing  for  correction  of  charge  size  if 
necessary.  After  the  shot  instant,  the  distances  to  the  receivers  were 
controlled  by  using  the  unprocessed  recordings;  all  data  were  bandpass- 
filtered,  2-14  Hz,  and  stacked  in  time-distance  plots.  The  response  of 
the  four  sensors  of  seismometer  No.  2,  for  different  runs  B  and  C,  are 
displayed  in  Figs.  3  and  4.  In  Fig.  5  we  see  the  response  of  the  z- 
component  (the  vertical  siismometer)  for  profiles  of  different  azimuth. 
We  can  trace  interface  wave  phases  up  to  a  range  of  4  km;  beyond  that  the 
8ignal-to-noi8C  ratio  becomes  too  small.  However,  even  for  the  close 
ranges  different  modes  cannot  to  identified  easily.  This  is  a  further 
indication  of  a  rather  complex  subbottom  structure.  We  marked  events  of 
the  interface  wave  train  and  correlated  them  »dth  group  velocities,  as 
shown  in  Table  l.  The  values  obtained  are  in  the  expected  range  and  typi¬ 
cal  for  the  environment  encountered.  They  show  a  certain  scatter, 
amounting  to  about  6Z  for  the  slow  end  and  35Z  for  the  higher  end  of  the 
interface  wave  train  for  profiles  A-B  and  C-D.  On  profile  F,  however, 
velocities  are  considerably  higher  (see  Table  1),  remarkable  because  pro¬ 
files  B  and  F  vary  in  direction  by  only  9°.  The  data  also  contains  p-vave 
phases  travelling  as  headwaves  or  refracted  waves  in  deeper  interfaces 
with  velocities  of  1.8  to  2.4  km/s.  These  data  will  be  included  when 
modelling  the  full  wave  field. 
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We  calculated  the  interface  wavelet  energy  generated  by  standard 
charges  for  different  runs  and  varying  ranges.  Having  only  a  few  profi¬ 
les,  we  obtained  information  on  the  distribution  of  interface  wave  energy 
for  selected  azimuthal  angles,  rather  than  a  continuous  azimuthal  distri¬ 
bution  of  energy.  We  compared  two  pronounced  frequency  bands  of  4-5  Hz 
and  8-9  Hz  for  three  different  ranges  (Fig.  6).  There  are  considerable 
differences  in  energy  levels,  ranging  as  high  as  18  dB.  At  this  stage, 
however,  the  results  apparently  are  contradictory,  with  inconsistency  in 
range  and  frequency. 


Fig.  2  Area  of  operation:  Ligurian  Sea  off  the  coast  of  Italy 


Whereas  we  can  observe  a  trend  in  the  behaviour  for  the  8-9  Hz  fre¬ 
quency  band  -  lowest  attenuation  in  direction  of  profile  D  and  highest 
attenuation  on  profile  B  (data  were  not  available  at  all  ranges  for  pro¬ 
file  F)  -  no  general  trends  can  be  observed  i or  the  4-5  Hz  frequency  band. 
As  with  the  group  velocity  analysis,  the  most  pronounced  effects  can  be 
observed  on  profile  F,  where  we  find  extremal  valued  for  the  attenuation, 
hignest  for  the  8-9  Hz  frequency  band  and,  at  the  same  range,  the  lowest 
for  the  4-5  Hz  frequency  band.  This  analysis  was  based  on  only  three 
ranges  and  two  frequency  groups;  further  analysis  should  reveal  a  more 
coherent  picture  of  spatial  energy  distribution  when  extended  to  full 
r^nge  and  all  frequencies. 
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Fig .  5  Stacks  of  seismograms  for  profiles  A,  B,  C,  D. 
Only  the  Z-component  is  displayed. 
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BEARING  OF  PROFILES  RELATIVE  TO  OBS 
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Fig.  6 

The  interface  mare  energy 
as  a  function  of  profile- 
bearing ,  frequency  and 
range.  All  data  are 
relative  to  Seismometer 
No.  2  and  calculated  for 
standard  bottom  charges 
(ISO  g). 
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A  further  approach  in  analyzing  the  interface  wave  trains  was  to 
apply  the  multiple  filter  technique  or  dispersion  analysis  suggested  by 
Dzlewonskl  [8J.  This  technique  provides  a  tool  to  separate  modes  of  dif¬ 
ferent  order  in  a  seismogram.  However,  this  procedure  has  limitations  in 
this  case  because  the  interface  wave  trains  display  only  a  poor  modal 
structure. 

n’e  will  present  herein  a  characteristic  case,  which  will  also  be  used 
in  modelling  the  data  by  applying  an  inverse  method.  The  dispersion  ana¬ 
lysis  was  run  and  displayed  in  a  Gabor  matrix  of  the  data  from  profile  C. 
Prior  to  the  analysis  a  conceptual  rotation  of  the  geophone  orientation 
was  performed,  i.e.  an  operation  was  applied  to  obtain  maximum  r-sponse 
of  the  radial  sensor  in  the  direction  of  the  profile.  We  determine  the 
radial  (R)  and  transverse  (T)  components  from  the  actual  X  and  Y  com¬ 
ponents  of  the  seismometer  for  a  given  source  bearing  6  by: 

R  =  x  cosfl  +  y  sin8 

T  =  -x  sinfl  +  y  cos8 

Hence,  only  Gabor  matrices  for  the  radial  and  vertical  components  are 
displayed,  (Fig.  7).  The  most  striking  feature  is  the  high  velocity  gra¬ 
dient  effective  over  a  very  narrow  frequency  band  of  2  Hz.  We  find  a 
reduction  in  group  velocity  from  about  260  m/s  to  74  ra/s  over  a  frequency 
range  from  3.5  to  5.5  Hz.  In  addition  to  this  we  find  some  further 
features  of  interest  on  the  radial  geophone  data,  marked  as  local  maxima, 
which  at  this  stage  are  taken  as  indicators  of  the  complexity  of  the  sub- 
bottom  structure. 

Earlier  experiments  by  Essen  et  al  [6],  Rauch  &  Schmalfeldt  [5]  and 
Schmalfeldt  &  Rauch  [7]  have  *  clearly  demonstrated  a  correlation  between 
the  physical  properties  of  the  seafloor  and  subbottom  and  the  propagation 
characteristics  of  interface  waves.  Therefore  in  this  area,  with  its  pro¬ 
nounced  lateral  variations  in  sedimentation  rates  and  processes,  the 
effects  depend  on  the  bearing  (direction)  of  the  profiles.  Further,  there 
is  always  a  certain  amount  of  scatter  involved  in  arrival  times  of  seismic 
waves  travelling  near  the  surface  of  the  crust,  or,  in  this  case,  in  the 
top  sediment  layers.  Because  of  this,  one  cannot  attribute  with  certainty 
the  observed  effects  to  anisotropy  or  to  lateral  heterogeneities. 


>IM.O 


S.o  8.7 
FREQUENCY  (HZ! 


FREQUENCY  (HZ) 


Fig.  7  Gabor  a ntrix  for  shot  Ho.  2  on  profile  r  fsee  Fig.  2)  for 
range  0.552  km.  Displayed  are  the  vertical  and  radial 
components . 
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Prior  to  the  actual  interface  wave  experiment  a  set  of  high  resolution 
seismic  UNIBOOM  profiles  were  run.  The  experiment  had  three  goals:  to 
determine  experimentally  relevant  coring  locations,  to  enable  later  corre¬ 
lation  of  geoacoustical  parameters  and  geological  features  and  to  map  the 
topography  of  the  top  layers. 


Two  typical  profiles  are  shown:  profile  A-B  with  a  bearing  of  320° 
(Fig.  8),  and  profile  C-D  with  a  bearing  of  25° ,  representing  a  sloping 


Fig.  8  Reflection  seismic  UNIBOOM  profile  A-B 


seabottom  (Fig.  9).  Profile  A-B  has  a  uniform  water  depth  (it  follows  the 
100  m  contour;  all  other  profiles  crossed  contour  lines). 

The  major  structural  feature  recorded  on  profile  A-B  is  a  slightly 
tilted  reflecting  surface  that  outcrops  on  the  slope.  This  reflector, 
which  is  interpreted  as  the  bedrock  surface  during  the  last  glacial 
period,  shows  fracturing  as  well  as  erosion.  It  is  overlain  by  sediment 
wedges  that  thin  out  towards  the  slope  and  consist  of  some  weaker  reflec¬ 
tors,  a  dark  band  of  sediments  with  some  surface  roughness  and  a  laminated 
top  layer.  In  general  we  can  distinguish  two  different  regions:  the 
western  part  with  a  more  regular  and  uniform  sedimentation  and  an  eastern 
'  part  that  displays  strong  lateral  variations  in  sediment  accumulation. 

*  This  coincides  with  results  from  the  seismic  runs,  which  indicate  that  at 

i  a  range  greater  than  2.5  km  eastwards  fror:  seismometer  No.  2  either  poor 

f  or  no  interface  wave  trains  were  recorded.  On  profile  C-D  the  most 

I  interesting  features  are  cross  bedding  and  the  slumps,  indicators  of 

recent  sediment  dynamics.  Marks  of  earlier  gravitational  sliding  as  well 
as  the  stratification  of  the  upper  sediments  can  be  clearly  identified. 

\ 

{  Eleven  cores  were  taken  with  a  recovery  range  from  2.0  to  3.95  m  and 

If  an  average  recovery  of  3.80  m.  Three  of  the  cores  were  taken  along  pro¬ 

file  A-B:  One  was  taken  at  the  site  of  the  seismometer,  one  on  the  slope, 
and  one  in  the  region  of  anomalous  sediment  accumulation  (see  Fig.  2). 
•  Eight  cores  were  taken  along  the  other  profiles  to  allow  for  spatial 

interpretation  of  sediment  distribution.  The  No.  7  core,  closest  to  the 
array,  has  been  analyzed  and  the  others  have  been  inspected. 
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The  core  analysis  shows  that  the  lower  end  (3.70-3.85  m)  has  a  sandy 
composition  (46. 8%  sand,  27.6%  silt,  25.5%  clay)  that  may  correlate  with 
sand  at  3.95  m  at  a  nearby  core  site.  We  find  some  dark  laminations  that 
correspond  to  regions  of  higher  p-wave  velocities  (see  Fig.  10).  From  the 
available  data  we  can  conclude  that  there  exists  c  sandy  interface  at 
about  4  m  of  depth.  We  can  correlate  this  e'*ent  to  the  reflector  S  marked 
in  Fig.  7. 
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The  continuous  seismic  reflection  profiles  indicate  strong  lateral 
heterogeneities.  It  therefore  can  be  assumed  that  the  observed  dependen¬ 
cies  of  interface  velocities  and  energy  on  tlte  direction  of  propagation 
are  caused  primarily  by  lateral  inhomogeneities.  However,  this  does  not 
rule  out  the  effects  of  anisotropy.  The  further  analysis  of  data  from  all 
three  ocean  bottom  seismometers  (OES)  should  allow  more  precise  discrimi¬ 
nation  of  the  various  effects  involved. 


MODELLING 


Fig.  10  a)  Site  of  the  seismometer  array  and  core  No.  7  with 
indicated  analyzed  column, 
b)  P-wave  analysis  of  core  No.  7 


The  inverse  method  for  calculating  the  dispersion  of  interface  waves 
trains  needs  realistic  input  data,  such  as  measured  velocity-depth  profi¬ 
les  for  P  and  S  waves.  We  have  at  this  stage  the  following  kind  of  data 
and  information: 

-  Sound  speed  profiles  in  the  water  column. 

-  A  velocity  depth  profile  for  P  waves  in  the  core. 

-  The  actual  thickness  of  the  uppermost  soft  layer. 

-  Evidence  for  a  sandy  substratum. 

-  Group  velocities  from  seismogram  stacks. 

-  Geoacoustic  properties  from  the  core  analysis. 

Using  these  data  we  applied  modelling  codes  from  which  the  group  and 
phase  velocities  for  different  modes  were  calculated  as  a  function  of  fre¬ 
quency  (Fig.  11). 


Fig.  11 

Calculated  dispersion  curve 
from  input  model  based  on 
experimental  results . 
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The  first  or  fundamental  modo  haa  bean  successfully  modelled.  The 
reaulta  fit  the  data  reaaonably  well.  However,  mode a  of  higher  order 
could  not  be  matched  to  the  data,  although  aoveral  mode  la  and  varlationa 

were  tried.  The  final  model  need  waa  a  aandy  half apace  aubbottom  with  a 

shear  velocity  of  Vfl-300  m/a,  overlain  by  a  4  m  thick  layer  of  soft  sedl- 
mente  with  a  ahear  velocity  profile  of  78  m/e  on  top  and  96  a /a  on  the 
bottom.  A  ahear  velocity  of  300  m/a  waa  uaed  for  the  halfapace,  alnce 

thla  haa  proved  to  be  a  realistic  value  for  other  measurements  In  this 
region . 

Not  all  features  seen  In  the  Gabor  matrix  display  can  be  modelled. 

This  Is  because  the  local  maxima  revealed  In  the  dispersion  analysis  of 
the  field  data  may  be  caused  by  structural  effects  of  the  subbottom  rather 
than  higher  order  modes.  Also,  the  input  model  needs  conaiderable  adjust¬ 
ment  to  the  real  world  in  order  to  allow  the  modelling  of  higher  order 
modes.  Nevertheless,  if  theoretical  reaulta  agree  reaaonably  with  field 
data  for  at  least  the  fundamental  mode,  then  the  geoacoustic  properties  of 
the  uppermost  layer,  that  is  the  dispersion  behaviour  of  Interface  waves, 
can  be  predicted  to  a  reasonable  degree  of  accuracy.  With  the  sensors 
sitting  on  the  top  layer  of  the  seafloor,  in  general  only  conversion  pro¬ 
cesses  at  the  seafloor  or  near  seafloor  interfaces  can  be  recorded  because 
Interface  wave  energy  decreases  exponentially  perpendicular  to  the  propa¬ 
gation  path. 

CONCLUSION 

By  evaluating  all  of  the  measured  data,  we  can  begin  to  resolve  the 
pertinent  problems  that  are  related  to  the  effects  and  the  role  of  the 
geological  environment  on  Interface  wave  propagation.  The  preliminary 
data  analysis  has  shown  that  a  multiparameter  survey  combined  with  a 
powerful  analysis  and  modelling  phase  can  provide  a  deeper  understanding 
of  the  relevant  conversion  processes.  This  will  lead,  to  a  greater 
accuracy  for  predictions  of  interface  wave  characteristics. 
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LOW-FREQUENCY  ANOMALIES  IN  THE  REFLECTION  BEHAVIOUR  OF  MARINE  SEDIMENTS 


Markus  von  Haumeder 
SACLANT  ASW  Research  Centre 
Vlale  San  Bartolomeo  400 
1-19026  La  Spezla,  Italy 


ABSTRACT 

It  is  frequently  found  that  the  reflectivity  of  the  sea  bottom 
deviates  considerably  from  the  values  that  are  expected  by  considering  its 
impedance.  A  type  of  sea  bottom  in  which  anomalies  in  the  reflectivity 
are  particularly  strong  is  a  coarse  marine  sediment  consisting  of  two  com¬ 
ponents:  particles  that  form  a  stiff  frame  and  a  liquid  that  fills  the 
interstices.  Such  sediments  can  be  considered  porous  materials.  A  theory 
that  analyses  the  response  of  a  porous  material  to  an  elastic  wave  is  the 
Biot  theory.  It  treats  the  two  components  as  two  interpenetrating  elastic 
contlnua  and  predicts  some  remarkable  anomalies  in  the  acoustic  behaviour 
of  porous  materials.  One  important  prediction  of  the  Biot  theory  is  that 
acoustic  energy  of  a  certain  frequency  band  can  penetrate  the  sediment 
even  in  the  range  of  total  reflection.  The  porous  material  behaves  here 
like  a  frequency  filter  with  respect  to  broadband  acoustic  energy.  This 
energy  transfer  into  the  porous  material  is  associated  with  a  relaxation 
process.  We  give  here  a  derivation  of  the  pertinent  relaxation  times  and 
demonstrate  how  the  relaxation  time  can  be  evaluated  from  the  reflectivity 
anomaly.  The  relaxation  time  is  expressed  in  terms  of  parameters  descrip¬ 
tive  of  the  internal  flow.  Use  of  this  interrelation  permits  us  to  esti¬ 
mate  the  flow-permeability  of  a  sediment  on  the  basis  of  a  bottom 
reflection  loss  measurement.  The  method  is  applied  to  a  real  marine  sedi¬ 
ment. 


INTRODUCTION 

A  large  fraction  of  the  naturally  occurring  marine  sediments  can  be 
considered  porous  materials.  Being  Immersed  in  a  fluid,  such  a  porous 
material  consists  of  two  clearly  separated  components:  a  porous  skeleton 
or  frame  and  a  fluid  that  fills  the  pore  space.  Each  of  the  two  com¬ 
ponents  forms  a  large,  interconnected  cluster  and  the  two  clusters  fit 
exactly  into  each  other.  This  situation  can  be  viewed  as  two  inter¬ 
penetrating  elastic  contlnua.  A  pressure  wave  impinging  on  such  a  com¬ 
posite  material  results  not  only  in  a  deformation  of  the  single  components 
but  also  in  a  relative  motion  of  the  two  contlnua  with  respect  to  each 
other.  The  fluid  is  free  to  move  through  the  frame  from  areas  of  high 
local  pressure  to  areas  of  low  local  pressure. 
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The  theory  that  describes  the  propagation  of  elastic  waves  through  a 
porous,  fluid-filled  material  was  developed  by  Maurice  Blot  [1-3].  The 
Biot  theory  Is  a  microscopic  theory  that  takes  the  particular  geometry  of 
the  porous  skeleton  Into  consideration.  It  Introduces  two  particle- 
displacement  vectors,  one  for  the  solid  frame  and  one  for  the  pote  fluid, 
and  expresses  the  Internal  flow  in  terms  of  a  difference  vector.  The  com¬ 
posite  nature  of  the  porous  material  is  reflected  in  the  theory  by  a  pair 
of  wave  equations  that  are  coupled  by  a  flow  term.  Thus,  the  Biot  theory 
is  capable  of  describing  the  underlying  microscopic  processes  in  some 
detai l . 


According  to  the  Biot  theory,  porous  materials  exhibit  some  anomalies 
in  their  response  to  a  pressure  wave.  One  of  these  anomalies  pertains  to 
the  reflection  process  at  the  interface  between  a  liquid  and  a  porous 
material.  In  an  angle-of-incidence  range  close  to  grazing  angle  where 
total  reflection  is  expected  for  a  solid  of  equivalent  acoustic  impedance, 
a  fluid-filled  porous  material  absorbs  energy  of  a  certain  frequency  band. 
The  porous  material  acts  here  as  a  band-pass  filter  with  respect  to  broad¬ 
band  acoustic  energy,  and  the  fraction  of  energy  that  is  transferred  into 
the  interior  of  the  porous  material  is  used  to  drive  the  internal  flow. 
For  many  frequencies  this  internal  flow  will  follow  the  excitation  by  the 
pressure  wave  with  a  delay.  Such  a  delay  results  in  energy  dissipation 
due  to  relaxation.  Relaxation  can  be  described  best  in  terms  of  a  charac¬ 
teristic  response  time,  the  so-called  relaxation  time,T. 

This  paper  gives  a  derivation  of  the  relaxation  times  for  a  porous, 
fluid-filled  material  on  the  basis  of  the  Biot  theory.  The  relaxation 
times  are  expressed  entirely  in  terms  of  flow  parameters,  which  gives  evi¬ 
dence  that  the  underlying  microscopic  process  is  a  flow  process. 

Further  we  present  a  measurement  of  the  bottom-reflection  loss  that 
shows  the  reflectivity  anomaly  and  permits  the  relaxation  time  to  be 
determined  experimentally.  This  measurement  is  used  in  combination  with 
the  given  formula  to  calculate  the  permeability/porosity  ratio  of  the  set 
bottom.  The  procedure  used  can  be  considered  a  new  technique  for  remote 
bottom-permeability  determination. 


1  THEORY 


The  Biot  theory  treats  an  isotropic,  homogeneous,  porous  material  and 
its  pore-filling  fluid  as  two  interpenetrating  elastic  contlnua  [1-5]. 
Each  of  the  two  contlnua  shows  its  own  characteristic  response  to  a 
pressure  wave,  and  the  composite  of  the  two  materials  requires  the  use  of 
two  particle  displacement  vectors,  one  for  the  skeletal  frame,  u,  and  one 
for  the  pore  fluid,  U.  The  difference  between  the  two  displacements  Is  a 
measure  of  the  internal  fluid  flow,  which  is  expressed  in  terms  of  the 


vector  quantity, w : 

w  =  0  [u  -  u]. 


(1) 


The  quantity  w  represents  the  average  fluid  displacement  relative  to  the 
frame  and  f)  is  the  porosity  defined  as  the  fraction  of  volume  occupied  by 
the  pore  fluid.  By  expressing  the  vectors  u  and  w  in  terms  of  scalar 
potentials,  <t>  ,  and  vector  potentials,  'If  , 

u  =  v4>s‘  curl  "  =  V4>,  ♦  curl 

the  Biot  equations  of  motion  [3]  can  be  seperated  [5] 


SACLANTCEN  CP-37 


10-2 


(2) 

into  a  pair  of 


VON  HAUMEDER:  Reflection  behaviour  of  marine  sediments 

equations  Involving  only  dllatational  modes  and  another  pair  of  equations 
for  the  rotational  mode.  The  subscript  s  refers  to  'solid'  and  the 
subscript  f  to  'fluid'.  We  start  our  consideration  with  the  simpler  set 
of  wave  equations  for  the  shear-wave: 


H  =  p  -  P,  .  Vk  'i',  =  Pt%-  % 


(3) 


where  p  is  the  shear  modulus,  p  and  p(  are  the  densities  of  the  fluid  and 
of  the  fluid-filled  porous  material,  respectively,  q  stands  for  the  fluid 
viscosity,  and  k  designates  the  flow  permeability.  The  composite  nature 
of  the  porous,  fluid-filled  material  is  reflected  in  the  theoretical 
description  by  a  set  of  two  coupled  wave  equations.  The  term  on  the  left- 
hand  side  of  the  second  equation  is  a  flow  term  that  has  to  be  compared  to 
the  lag  term  in  relaxation  theory.  This  term  contains  the  particular  phy¬ 
sics  of  porous  materials.  Its  prefactor  p/k  is  the  viscous  resistance  to 
fluid  flow.  For  high  frequencies  the  fluid  viscosity,  jj  ,  is  not  constant 
and  the  viscous  resistance  to  fluid  flow  must  be  made  frequency  dependent. 
Considering  the  friction  forces  within  pore  channels  of  simplified 
geometry.  Blot  derived  [1]  a  complex  correction  factor,  F(x)  ■  Fr(«)  + 
'F|(*)»  for  oscillatory  motion.  The  factor  7j/k  may  be  replaced  by 
Fij/k  in  the  frequency  domain,  where 


F(«) 


1  K  T|«| 

<  l-Tj«|/i  K 


berffc)  ♦  i  bei  [k] 
ber[x]  .  i  bei  [k] 


(4) 


The  functions  ber(x)  and  bei(x)  are  the  real  and  imaginary  parts  of 
the  Kelvin  function  with  ber’(  x  ),  bel'(x)  being  their  derivatives.  The 
argument,  <s  ,  of  the  correction  factor,  F, 


k  =  a  yj  pt  to/ Tj ' 


(5) 


is  non  dimensional  and  depends  on  the  pore-size  parameter  a,  the  fluid 
density  p,  ,  the  viscosity,  n  ,  and  the  angular  frequency, w •  Since  we  are 
going  to  make  use  of  this  correction  factor  F(»t)  In  a  quantitative  way,  we 
show  the  variation  of  its  real  and  imaginary  part  with  frequency  in 
Fig.  1. 


rig.  1 

Real  p art  (Fr)  and  imaginary 
part  (Ff )  of  the  frequency 
correction  function  for 
three-dimensional ,  oscillatory 
flew. 


For  very  low  frequencies  F(<)  approaches  unity,  resulting  In  the  same 
equation  as  when  Poiseullle  flow  Is  assumed.  For  high  frequencies,  and 
hence  large  values  of  x  ,  the  curves  become  straight  lines  with  a  slope  of 
1/4  fl*  0.177.  The  derivation  of  the  complex  correction  factor  Is  valid 
only  for  frequencies  where  the  wave  length  is  large  compared  with  the 
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pore  size.  For  sands,  this  puts  the  upper  Unit  on  frequencies  at  about 
10s  to  106  Hz. 

Incorporating  the  complex  correction  factor  into  Eq.  3  by  replacing 
q  by  q  F(x),  we  derive  now  an  explicit  relationship  between  attenuation 
and  frequency. 


By  Fourier  transforming  Eq.  3  into  the  frequency  domain  and 
assuming  plane-wave  type  potentials,  the  two  equations  can  be  combined  to 
give  one  homogeneous  wave  equation  of  the  form: 


V2* 


u*.-  _  fl2o> Vm 

V2  Ptosz/S  -  ionjF/k 


0 


which  permits  solutions  of  the  type 


(6) 


=  A  exp  i(u>t  «-kx)  ♦  B  exp  i[o>t  -kx)  . 


(7) 


Equation  6  can  be  cast  into  the  form  of  a  lossy  Helmholtz  equation  of  the 

•*•2.  2 

general  type  (a2  +  Tc2 )  *  0,  with  k  =  (k  -io )  .  In  this  fora  it  appears 
as 


-  iojqF/k 


}]*- 


(8) 


The  complex  k-vector,  k,  is  to  be  identified  with  the  expression  in 
the  brackets.  Separating  the  real  and  iaaglnary  parts  of  the  wave  vector 
renders  two  equations  from  which  we  derive  the  attenuation  coefficient: 


“  ~  \/7r  v^l^^7FrSKajV:*  '^1 ; %]- kHT}  (9) 


Pi  r  _  P  _k_ 

P  j3  V* 


(10) 


in  which  r  has  the  dimension  of  a  time  and  can  be  identified  as  a  relaxation 
time.  This  typical  response  time  of  the  porous  material  is  entirely  expressed 
in  terms  of  parameters  that  are  descriptive  of  the  internal  flow  of  the  pore 
fluid  through  the  frame,  i.e.  permeability,  k,  fluid  density,  viscosity, 
and  the  real  part  of  the  frequency  correction  factor,  Fr  .  K  is  a  dimen¬ 
sionless  factor  that  is  determined  by  the  density  ratio  of  the  fluid  and 
the  porous  material  (  p-0pf*l>-B)Pr  ).  This  factor  Introduces  an  altered 
frequency  into  the  attenuation  formula,  which  largely  affects  the 
reflection  coefficient. 


If  we  now  assume  the  resistance  to  fluid  flow  to  be  frequency  inde¬ 
pendent  by  setting  the  correction  factor,  F(*)»l,  we  obtain  for  the  atte¬ 
nuation  the  simpler  formula: 


Equation  9,  which  is  valid  for  all  frequencies,  takes  the  form 
of  Eq.  II  idien  »— 0,Fj/Fr  —  0.  For  the  case  of  a  pure  liquid  (/9«1),  K 

equals  one  and  Eq.  11  takes  the  fora  of  Che  well-known  formula  for  the 
relaxation  attenuation  in  a  viscous  fluid: 


a 


_ i _ 

I  ♦  u*rJ 


1  ♦ 


i  * 


(12) 
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This  shows  that  the  derived  formula  is  capable  of  producing  the 
proper  results  in  the  limit  3  =  1*  We  now  want  to  test  in  how  far  Eq.  11 
reproduces  the  overall  variation  of  the  attenuation  as  obtained  by  the 
Biot  computations. 


For  this  purpose  the  attenuation  according  to  Eq.  11  is  plotted 
versus  the  frequency  for  a  porosity  of  0.5  (curve  *b’  of  Fig.  2).  In  this 
case  of  a  frequency-independent  flow  resistance  (Poiseuille  flow)  the 
attenuation  tends  towards  a  constant  value  at  higher  frequencies.  Curve 
’a*,  however,  represents  a  plot  of  Eq.  9  in  which  the  flow  resistance  is 
made  frequency  dependent.  The  important  effect  of  a  frequency-dependent 
flow  resistance  is  apparent,  mainly  at  high  frequencies.  These  curves  are 
in  very  good  agreement  with  corresponding  curves  shown  by  Stoll  15,6), 
which  were  obtained  as  a  result  of  a  Biot  model  computation. 


FREQUENCY  (Hz) 


The  agreement  gives  evidence  that  the 
derived  formula  reflects  the  attenuation 
in  porous  materials  as  considered  by  the 
Biot  theory.  In  the  following  chapters 
we  make  use  of  the  bottom-reflection  loss 
anomaly  in  order  to  determine  t  experi¬ 
mentally. 


Fig.  2 

Frequency  dependence  of  the  attenuation  for 
a  frequency  dependent  riov  resistance  (a) 
and  a  constant  flat  resistance  (b) . 


2  BIOT  PREDICTIONS  OF  BOTTOM-REFLECTION  LOSS 

The  frequency  dependence  of  the  bottom-reflection  loss  exhibits  cae 
of  the  most  significant  anomalies  in  porous  materials.  In  an  angle-of- 
incidence  range  where  total  reflection  would  normally  be  expected,  a 
liquid/porous-composite  interface  acts  as  a  band-pass  filter  with  respect 
to  broadband  energy.  Figure  3  shows  the  Biot  predictions  for  the  fre¬ 
quency  dependence  of  the  reflection  coefficient  for  different  grazing 
angles.  At  high  grazing  angles  the  energy  is  reflected  almost  indepen¬ 
dently  of  frequency.  At  grazing  angles  below  about  10*,  however,  energy 
of  a  certain  frequency  band  penetrates  the  porous  material,  whereas  the 

energy  of  higher  and  of  lower  frequen¬ 
cies  is  reflected.  The  line  shape  of 
the  minimum  resembles  that  of  an 
Inverted  relaxation  curve.  In  order  to 
verify  that  the  internal  f  low- 


Fig.  } 

Frequency  dependents  of  the  reflection 
coefficient  for  different  grazing  angles 
Wc  ¥3*  it1  10 1  10*  as  predicted  by  the  Plot  theory - 

FftCOUCMCY  (Mi) 

relaxation  creates  such  a  minimum  in  the  reflectivity  we  display  the 
relaxation  and  reflectivity  curves  together  in  Fig.  4.  A  maximum  of  the 
relaxation  manifests  Itself  most  clearly  in  the  logarithmic  decrement  of 
the  viscous  losses,  which  are  plotted  in  the  lower  part  of  Fig.  4  for  dif¬ 
ferent  permeability/ pore-size  combinations.  The  appearance  of  the  maximum 
in  the  relaxation  Is  closely  correlated  to  the  appearance  of  the  minimum 
in  the  reflection  loss. 
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By  varying  the  permeability,  the  curve  is  strongly  shifted 
within  the  frequency  spectrum;  for  comparison  we  display  on  top  of  Fig.  4 
the  corresponding  minima  in  the  bottom-reflection  loss.  It  turns  out  that 
the  corresponding  minima  and  maxima  do  not  show  up  at  the  same  frequency. 
Instead  they  are  offset  by  a  constant  ratio:  the  frequency  at  which  the 
minimum  in  the  bottom-reflection  1-^ss  occurs  is  0.53  times  the  frequency 
at  which  the  maximum  of  the  logarithmic  decrement  is  observed.  This  is  a 
direct  result  of  the  Biot  theory  and  will  be  used  later  for  the  evaluation 
of  T. 

The  location  of  the  maximum  of  the  logarithmic  decrement  in  the  fre¬ 
quency  spectrum  permits  T  to  be  determined  via  the  peak-condition 

ft 


10'  10*  10*  10* 
FREQUENCY  (H*) 


Fig.  4 

Comparison  of  the  viscous  loss 
per  wave-length  (log. decrement  I 
with  corresponding  reflection 
loss  for  four  different  perm ea- 
bilitg(k) /pore-size(a)  combi¬ 
nations: 
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10 
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IV  : 

• 
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He  are  now  in  the  position  to  use  the  minimum  in  the  reflection  loss 
to  find  the  maximum  of  the  corresponding  logarithmic  decrement,  and  from 
the  frequency  where  this  maximum  occurs,  to  Jind  the  relaxation  time.  In 
principle,  the  minimum  in  the  bottom-reflection  loss  can  be  measured, 
although  it  is  not  clear  whether  the  Biot  effect  is  large  enough  to  be 
measurable  in  a  real  marine  sedlaer.* 


3  MEASUREMENT  OF  BOTTOM- REFLECTION  LOSS 

In  1978  the  SACLANT  ASH  Research  Centre  had  used  a  vertical  array  of 
twelve  hydrophones  to  measure  the  bottom-reflection  loss  in  certain  „eep 
basins  of  the  central  Mediterranean  at  frequencies  from  50  H*  to  7  kHz 
f7J.  To  measure  the  reflection  losses  over  this  wide  frequency  band, 
explosive  sources  of  500  g  TNT  were  used.  The  experimental  set-up  Is 
sketched  in  Fig.  j.  The  receiving  ship  suspended  a  750-m  long  vertical 
hydrophone  string. 
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RECEIVER 

ARRAY 


Fig.  5 

The  experimental  set-up  for 
botto-r-re  flection  loss  *  X* 
surMMti  17}. 


while  the  source  ship  moved  on  a  predetermined  fixed  course,  launching 
charges  set  to  explode  at  550  a  depth.  The  direct  and  reflected  acoustic 
signals  were  recorded  In  digital  form  by  the  receiving  shtp  together  with 
a  radio  pulse  transmitted  from  Che  source  ship  at  the  moment  when  It 
received  the  direct  acoustic  pulse  from  tt>e  charge.  Some  of  the  data 
collected  in  those  studies  will  be  used  here. 
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Because  the  sea  bottom  often  exhibits  surface  roughness  and/or 
layering,  the  reflection  anomaly  studied  in  this  paper  is  frequently 
masked  by  distortion  of  the  acoustic  signals.  To  isolate  the  reflection 
features  related  to  the  Biot  effects  from  those  related  to  the  bottom 
morphology  we  choose  a  measuring  site  that  shows  a  negligible  degree  of 
geomorphclogical  discontinuities.  This  choice  was  made  by  examining 
underwater  photographs  and  profiles  of  relative  sound  speed  (0.93),  rela¬ 
tive  density  (1.5),  and  porosity  (0.5),  deduced  from  core  samples.  The 
area  selected  for  the  measurements  used  in  this  paper  was  on  the  Messina 
abyssal  plain  in  comparatively  deep  water  (3700  m).  The  morphology  of  the 
sea  floor  over  the  abyssal  plain  is  extremely  flat,  with  a  very  soft  top 
layer  of  fine  clay.  The  area  is  homogeneous  enough  to  permit  the  measure¬ 
ment  of  the  Biot  reflectivity  anomaly. 

Accordingly  the  signals  themselves  are  undistorted,  almost  textbook¬ 
like,  which  easily  permits  Fourier  analysis.  As  one  result  of  this  analy¬ 
sis,  we  obtain  the  Impulse  response  for  different  angles  of  incidence,  as 
shown  in  Fig.  6. 


TIME  (ms) 


Fig.  6 

Defiled  displsy  of  the  first  20  ** 
of  impulse  responses  for  different 
g rezing  engles  [ 7 ). 


An  analysis  of  the  signal  pattern  indicates  a  very  soft  sea -bottom. 
This  bottom  type,  with  high  permeability  and  high  porosity,  implies  that 
the  sound  speed  in  the  bottom  is  less  than  that  in  sea  water,  which  repre¬ 
sents  an  intromission-angle  situation  for  the  frequencies  discussed  here. 
Accordingly,  the  first  reflection  remains  small  for  all  grazing  angles. 
However,  the  first  layer  is  terminated  at  a  depth  of  14  m  by  another, 
high-velocity  layer,  as  can  be  deduced  from  the  signal  pattern  in  Fig,  6 
(extrapolation  to  vertical  incidence).  As  long  as  the  grazing  angle  is 
less  than  10*,  all  reflections/refractions  merge  here  to  create  a  single, 
large-amplitude,  phase-shifted  pulse.  It  is  mainly  at  this  14-m  Interface 
that  a  reflectivity  anomaly  similar  to  the  one  predicted  by  the  Blot 
theory  was  measured. 


4  COMPARISON  BETWEEN  MODEL  PREDICTIONS  AND  MEASUREMENTS 

By  Fourier  analyzing  the  measured  signal  we  obtain  the  bottom- 
reflection  loss  at  different  frequencies.  If  we  use  the  whole  signal  for 


SACLANTCEN  CP-37 


10-7 


VON  HAUMEDER:  Reflection  behaviour  of  marine  sediments 


the  computation,  i.e.  the  full-length  signal  of  Fig.  6  comprising  the 
contributions  of  all  layers,  we  obtain  the  oscillating  curve  depicted  in 
Fig.  7.  the  bottom-reflection  loss  is  here  shown  for  a  very  small  grazing 
angle  of  10°. 

It  is  striking  that  the  overall  tendency  of  this  curve  resembles  very 
much  the  one  predicted  by  the  Biot  theory  for  liquid/porous-composite 
interfaces,  as  shown  in  Fig.  3.  Some  features  are  different,  due  to  the 
layering  and  mixing  of  different  materials  in  a  real  sea  floor  and  due  to 
the  fact  that  the  signal  is  also  effected  by  the  refracted  arrival. 
Nevertheless,  an  appropriate  choice  of  the  input  parameters  to  the  Blot 
model  permits  the  measured  curve  to  be  fitted  in  a  way  that  reproduces  the 
overall  tendency;  this  is  particularly  true  for  the  position  of  the  mini¬ 
mum  which  is  strongly  shifted  in  frequency  by  small  changes  of  the  per¬ 
meability  (see  Fig.  4).  The  superimposed  solid  line  in  Fig.  7  represents 
the  fit  by  the  Biot  theory.  The  values  of  the  parameters  used  foi  this 
computation  are  given  in  Table  1  and  were  chosen  to  obtain  agreement  bet¬ 
ween  the  experimental  and  the  theoretical  curves  in  Fig.  7.  However,  in 
choosing  their  values,  constraints  were  imposed  in  order  to  keep  them 
within  physically  realistic  limits.  Accordingly  they  are  of  the  same 
order  of  magnitude  as  commonly  used  parameter  sets  given  by  Stoll  [81. 


Fig.  7 

Bottom-reflection  loss  versus 
frequency  for  a  20*  grazing  angle. 
The  smooth  line  represents  the 
prediction  by  the  Biot  theory. 


TABLE  1 


POROSITY 

f}- 0.5 

ISOTROPY  FACTOR 

a*  3 

PERMEABILITY 

BULK  MODULUS  OF  FRAME 

K0»  «■»'(■» 

PORE  SIZE 

a  s  1  -  to'5  m 

SHEAR  MODULUS  OF  FRAME 

fl.  2.5-to’p* 

DENSITY  OF  GRAIN  MATERIAL 

2.5 -to’  kg/mJ 

LOG  DECREMENT  FOR  SHEAR 

d,’  o.ots 

VIBRATIONS  OF  THE  FRAME 

DENSITY  OF  FLU© 

Qj*  10J  k9/W 

LOG  DECREMENT  FOR  BULK 

d'  0.015 

VISCOSITY 

1)  s 

VIBRATIONS  OF  THE  FRAME 

Regarding  the  position  of  the  minluoum  in  the  frequency  spectrum, 
which  is  critically  influenced  by  the  input  parameters,  the  fit  can  be 
considered  accurate.  From  this  fit  we  obtain  values  for  the  bottom  para" 
meters,  in  particular  for  those  determining  the  above-derived  relaxation 
time,  7  .  This  allows  us  to  apply  our  concept  of  a  relaxation  time 
directly  to  a  measurement.  On  the  one  hand  the  relaxation  time  is  known 
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by  the  frequency  at  which  the  measured  minimum  occurs  according  to  the 
peak  condition:  7  »  (i.l-F|/Ft)/?7T(i.8l),  where  the  factor  1.8  Is  determined 
by  the  shift  factor  of  Eq.  10.  On  the  other  hand  we  can  easily 
calculate  T  by  Introducing  the  fit  parameters  Into  our  derived  formula.  A 
comparison  of  the  two  values  of  T  shows  close  agreement  and  proves  the 
applicability  of  the  derived  formula  to  real  marine  sediments. 
Experimentally  we  obtain  7  »  1.5xI0-ss;  the  formula  gives  7  -  1.6xl0-5s. 
The  close  agreement  convinces  us  that  we  are  dealing  with  a  measured  Biot 
effect.  It  opens  up  the  possibility  of  determining  bottom  parameters  from 
a  measurement  of  the  bottom-reflection  loss. 


CONCLUSIONS 

The  derived  expression  for  the  relaxation  time  —  which  correlates  an 
overall  response  time  with  the  parameters  of  the  microscopic  structure- 
can  be  used  for  determining  bottom  parameters  remotely.  The  idea  is  first 
to  use  the  minimum  in  a  measured  bottom-reflection-loss  curve  to  evaluate 
the  relaxation  time  and  then  to  employ  Eq.  10  to  deduce  values  for  the 
correlated  structure  parameters.  The  relaxation  time  itself  is  expressed 
in  terms  of  the  permeabillty/poroslty  ratio  and  of  additional  parameters 
descriptive  of  the  liquid.  With  a  simple  assumption  concerning  the  rela¬ 
tion  between  flow  permeability  and  porosity  we  are  thus  in  a  position  to 
determine  the  flow  permeability  of  the  bottom  under  examination.  This 
method  of  remotely  sensing  the  permeability  may  not  only  be  useful  for 
characterizing  the  sea  bottom,  but  may  also-  have  applications  to  the 
sensing  of  liquid-carrying  layers.  By  analyzing  bottom-reflection 
losses  in  the  central  Mediterrane--’  w  have  demonstrated  the  applicability 
of  this  method  to  an  actual  sea-floor  measurement. 
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